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Chapter  1 
Introduction 


1.1  Motivation 

There  are  times  when  it  is  desirable  to  have  a  fluid  flow  in  a  highly  agitated,  tur¬ 
bulent  state.  In  such  cases,  the  early  transition  of  flow  from  laminar  to  turbulent 
is  actively  sought.  Dimples  on  a  golf  ball  and  mixing  vanes  in  jet  engine  combus¬ 
tors  and  industrial  furnaces  are  examples.  In  the  first  case,  the  turbulence  delays 
boundary-layer  separation,  thus  lowering  the  pressure  drag  and  allowing  the  ball  to 
fly  farther.  In  the  second  case,  the  turbulence  enhances  the  mixing  of  oxidizer  and 
fuel,  increasing  the  combustion  efficiency.  On  the  other  hand,  turbulence  can  also  be 
a  hindrance.  Examples  here  include  the  boundary  layers  on  aircraft  wings  and  ship 
hulls,  where  the  turbulence  increases  the  skin  friction  drag  and  lowers  fuel  efficiency. 
The  actions  of  the  turbulent  motions  are  also  transmitted  as  forces  on  the  constrain¬ 
ing  surfaces,  which  can  couple  with  the  vibrational  modes  of  the  structure  [1.1].  The 
result  can  be  accelerated  structural  fatigue  and  increased  radiated  noise,  both  inside 
and  outside  of  the  structure.  These  examples  on  the  negative  effects  provide  the 
motivation  to  control  the  turbulence  to  mitigate  or  eliminate  the  undesirable  effects. 

Efforts  to  control  turbulent  flows  have  focused  primarily  on  the  manipulation  of 
quasi-periodic,  coherent  flow  structures  [1.2],  since  it  is  generally  accepted  that  such 
structures  are  responsible  for  the  production  and  maintenance  of  boundary-layer 
turbulence  [1.3].  Most  of  the  work  to  date  has  been  in  the  area  of  passive  control, 
where  the  inputs  to  the  flow  boundary  conditions  are  static.  In  contrast,  active 
control  involves  dynamic  inputs  to  the  flow.  The  topic  has  been  recently  reviewed 
by  Bushnell  and  McGinley  [1.4],  who  point  out  that 

The  probable  key  to  increasing  the  effectiveness  of  dynamic  inputs  ...  is 
correct  phasing.  The  turbulence  dynamics,  due  to  phase  jitter  and  vari¬ 
ability  in  size/strength  of  the  various  motions,  occurs  randomly  in  both 
space  and  time.  . . .  sensor-feedback  control  loops  are  evidently  required 
for  treating  the  random-phase  problem  in  the  turbulence  dynamics. 

For  a  feedback  control  system,  it  is  desirable  that  the  sensor  be  nonintrusive. 
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This  naturally  leads  one  to  consider  flow  variables  at  the  wall,  of  which  the  fluctuat¬ 
ing  wall  pressure  is  the  most  readily  measured.  If  the  turbulence  is  to  be  controlled 
through  the  sensing  of  the  wall  pressure,  it  is  imperative  that  the  relationship  be¬ 
tween  the  two  be  thoroughly  understood.  Numerous  investigations  have  established 
a  link  between  near-wall  shear  layers  associated  with  turbulence  production  and 
wall-pressure  peak  events;  the  so-called  bidirectional  relationship  (cf.  Section  1.2.2). 
However,  these  studies  were  for  fully  developed  turbulent  boundary  layers  flowing 
over  smooth  flat  plates  with  zero  pressure  gradient,  or  for  fully  developed  turbulent 
flow  through  smooth  pipes.  Real-world  applications  will  rarely  involve  such  a  flow. 
A  host  of  disturbances  to  the  ideal  case  are  possible,  indeed,  probable.  The  present 
research,  therefore,  has  two  broad  aims:  to  clarify  the  causal  relationship  between 
the  organized  structures  and  the  wall-pressure  events  for  equilibrium  flows,  and  to 
examine  this  relationship  in  a  disturbed  flow. 


1.2  Background 


1.2.1  Physical  features  of  turbulent  boundary  layers 

The  following  overview  of  the  structural  features  of  the  turbulent  boundary  layer 
(TBL)  is  based  primarily  on  the  summary  work  of  Robinson  [1.3]  and  Kline  and 
Robinson  [1.5].  Organized  motions  or  coherent  structures  are  defined  as  spatial  re¬ 
gions  of  the  flow  over  which  at  least  one  fundamental  flow  variable  exhibits  significant 
correlation  with  itself  or  with  another  variable.  The  inner  layer  or  near-wall  region 
refers  to  the  region  y+  =  yuT/v  <  100,  where  y  is  the  distance  from  the  wall  in  the 
wall-normal  direction,  uT  is  the  shear  velocity,  and  v  is  the  kinematic  viscosity.  It 
includes  the  sublayer,  the  buffer  region,  and  at  least  part  of  the  logarithmic  region. 
The  rest  of  the  boundary  layer  is  called  the  outer  layer. 

To  understand  the  mechanisms  of  turbulence  production  and  dissipation,  one 
must  examine  the  governing  equations.  The  total  velocity  in  the  ith  direction  is 
represented  as  a  mean  and  a  fluctuation,  Ui  +  tq.  If  this  representation  is  substi¬ 
tuted  into  the  Navier-Stokes  equations  for  incompressible  flow,  the  equation  for  the 
turbulence  kinetic  energy  Kj  =  ( u\  +  u\  +  u\)/2  =  rquj/2  can  be  obtained  as: 
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—UiUj 
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dxj 


dui  dui 
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where  D/Dt  is  the  substantial  derivative,  p  is  the  density,  and  p  is  the  fluctuating 
pressure.  Overbars  represent  time  averages.  In  this  expression,  the  first  term  on 
the  right-hand  side  represents  turbulence  production  driven  by  the  mean  shear  and 
the  second  term  represents  dissipation  through  viscosity.  The  third  term  in  brackets 
denotes  three  forms  of  turbulence  transport.  The  goal  of  active  control  is  to  inhibit 
the  production  term. 
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For  positive  production  to  occur,  the  Reynolds  stresses  —  uiu]  must  be  positive, 
which  means  that  the  fluctuating  components  must  be  of  opposite  sign.  In  wall- 
bounded  flows,  the  mean  shear  is  dU/dy,  where  U\  =  U  is  the  mean  streamwise 
velocity  and  x2  =  y.  Production,  then,  occurs  if  iq  =  u  <  0,u2  =  v  >  0  or 
u  >  0,  v  <  0,  which  correspond  to  the  second  and  fourth  quadrants  of  the  u-v  plane. 
The  first  case  is  referred  to  as  a  Q2,  uv2,  or  ejection  motion.  The  second  case  is 
referred  to  as  a  Q4,  or  sweep  motion. 

The  production  and  dissipation  of  turbulence  in  wall-bounded  flows  is  carried 
out  predominantly  through  quasi-periodic,  organized  motions.  Across  the  boundary 
layer,  these  coherent  motions  exhibit  a  broad  range  of  spatial  and  temporal  scales.  In 
general,  they  consist  of  both  inner-layer  motions  and  outer-layer  structures  that  act 
across  the  complete  boundary  layer.  The  following  (unsorted)  classification,  taken 
directly  from  Robinson  [1.3,  pp.  606-607],  summarizes  the  experimentally  observed 
forms  of  coherent  motions: 

1.  Low-speed  streaks  in  the  viscous  sublayer 

2.  Ejections  of  low-speed  fluid  outward  from  the  wall,  including  lifting  low-speed 
streaks 

3.  Sweeps  of  high-speed  fluid  inward  toward  the  wall,  including  inrushes  from  the 
outer  region 

4.  Vortical  structures  of  various  forms 

5.  Sloping  near- wall  shear  layers,  exhibiting  local  concentrations  of  spanwise  vor- 
ticity  and  du/dx 

6.  Near-wall  “pockets”  visible  in  the  laboratory  as  regions  swept  clean  of  marked 
near-wall  fluid 

7.  Large  (5-scale)  motions  capped  by  three-dimensional  bulges  in  the  outer  tur¬ 
bulent/potential  interface 

8.  Shear-layer  “backs”  of  large-scale  outer-region  motions,  consisting  of  sloping 
(5-scale)  discontinuities  in  the  streamwise  velocity 

The  inner-layer  motions,  items  1-6,  are  responsible  for  the  majority  of  the  turbulence 
production.  They  are  termed  active  motions,  and  occur  intermittently  over  relatively 
small  time  and  space  scales.  In  contrast,  the  temporal  and  spatial  scales  of  the 
organized  structures  in  the  outer  region  of  the  TBL  are  much  larger  and  roughly 
scale  on  the  freestream  velocity  and  boundary-layer  thickness.  These  structures  are 
termed  passive  because  they  are  not  directly  linked  to  turbulence  production. 

The  first  three  items  of  the  classification  are  related  to  what  has  been  historically 
called  the  bursting  process.  The  term  “bursting”  originally  described  the  oscillation 
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and  violent  breakup  of  a  low-speed  streak  after  lifting,  followed  by  an  inrush  of  fluid, 
or  sweep.  The  emphasis  was  on  a  process  that  is  highly  intermittent  in  time.  An 
alternative  viewpoint  has  evolved  in  which  the  burst  is  an  ejection  caused  by  the 
passage  of  a  tilted  vortex,  which  is  localized  in  space  but  persists  for  long  times. 
Such  a  region  passing  a  stationary  probe  would  produce  sudden  velocity  excursions 
at  the  measurement  location  without  making  clear  the  passage  of  an  associated  flow 
structure.  The  disparate  meanings  of  “bursting”  throughout  the  literature  prompted 
Kline  and  Robinson  to  exclude  the  term  from  the  above  classification  [1.3]. 

Vortices  and  vortex  dynamics  (item  4)  play  a  central  role  in  the  description  of 
turbulent  flows.  At  a  fundamental  level,  the  transfer  of  turbulent  energy  from  large 
to  small  scales  occurs  through  vortex  stretching  [1.6].  In  this  process,  vortices  are 
elongated  by  the  mean  shear.  As  the  cross-sectional  area  of  the  vortex  decreases, 
conservation  of  angular  momentum  leads  to  an  increase  in  vorticity  and  kinetic 
energy.  The  stretched  vortex  induces  a  strain-rate  field,  which  stretches  vortices  of 
different  orientation.  As  the  process  continues,  kinetic  energy  is  redistributed  among 
the  coordinate  directions,  leading  to  local  isotropy  at  the  smallest  scales. 

Vortical  structures  have  been  identified  by  numerous  investigations,  and  are  the 
basis  for  most  conceptual  models  of  the  turbulence  production  dynamics.  Head  and 
Bandyopadhyay  [1.7],  using  motion-picture  flow  visualization  and  hot-wire  measure¬ 
ments,  documented  the  existence  of  large  numbers  of  vortex  structures  in  the  TBL. 
The  vortices  originate  at  the  wall  and  extend  through  a  large  part  of  the  boundary 
layer  or  beyond.  At  low  Reynolds  numbers,  these  vortices  form  horseshoe  or  loop 
shapes.  At  high  Reynolds  numbers,  they  become  elongated  into  hairpin  vortices  or 
vortex  pairs. 

Thomas  and  Bull  [1.8]  synthesized  a  model  for  the  large-scale  motions  associated 
with  the  burst-sweep  cycle  from  simultaneous  measurements  of  velocity,  wall  shear 
stress,  and  wall  pressure.  Figure  1.1  shows  plan  and  side  views  of  their  model. 
The  most  prominent  features  are  a  (5-scale  horseshoe  vortex  and  an  inclined  shear 
layer.  Recently,  Snarski  and  Lueptow  [1.9]  inferred  a  similar  structure  from  detailed 
correlations  between  the  wall  pressure  and  streamwise  velocity.  Their  model  includes 
secondary,  counter-rotating  vortices  between  the  fronts  and  backs  of  the  primary 
large-scale  structures. 

Robinson  [1.3,  1.10]  constructed  a  model  for  low-Reynolds  number  TBLs,  which 
was  deduced  from  a  review  of  both  experimental  and  direct  numerical  simulation 
data.  Inclined  quasi-streamwise  and  spanwise  vortices,  sometimes  forming  horseshoe 
or  hairpin  shapes,  form  the  core  structures.  Quasi-streamwise  vortices  dominate  the 
inner  layer,  while  archlike  vortices  are  the  most  common  structure  in  the  outer  layer. 
In  the  logarithmic,  or  overlap,  region,  both  types  of  vortex  exist  together,  often  as 
elements  of  the  same  structure.  The  two  types  of  structure  and  the  kinematical 
relationships  between  ejection/sweep  motions  are  illustrated  in  Figure  1.2. 

Even  with  a  wide  consensus  on  the  above  boundary-layer  structures,  there  are 
still  open  questions.  In  particular,  the  coupling  between  the  inner  and  outer  flow 
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Figure  1.1:  Plan  and  side  views  of  the  large-scale  structure  associated  with  the 
burst-sweep  cycle  in  a  frame  of  reference  moving  with  the  structure.  From  Thomas 
and  Bull  [1.8]. 


Figure  1.2:  Vortical  structures  and  ejection/sweep  motions  in  the  turbulent  bound¬ 
ary  layer.  From  Robinson  [1.3]. 
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motions  is  not  completely  understood.  On  the  one  hand,  the  inner-region  production 
should  intuitively  scale  on  the  wall  variables  uT  and  v.  Indeed,  some  researchers  have 
found  that  the  rate  of  bursting  near  the  wall,  usually  deduced  by  the  passage  of 
streamwise  shear  layers,  is  constant  when  scaled  on  wall  variables  [1.11].  However, 
others  have  concluded  that  the  burst  rate  is  linked  to  the  outer  flow  [1.9,  1.12]. 
Addressing  this  inner/outer  controversy,  Kline  and  Robinson  [1.5]  conclude  that 
while  the  majority  (70  percent  or  more)  of  the  turbulence  production  in  the  inner 
layer  is  self-maintaining  by  events  in  the  inner  layer  itself,  it  can  be  significantly 
affected  by  outer-layer  or  freestream  motions.  This  raises  the  following  question: 
Are  some  flow  structures  unaffected  by  outer-layer  motions,  and  hence,  universal  to 
turbulent  boundary  layers,  or  are  all  motions  modified  in  some  way?  The  answer  has 
implications  for  the  application  of  turbulence  control  in  the  presence  of  real-world 
disturbances,  since  the  control  strategy  may  have  to  change  with  the  flow  situation. 


1.2.2  Wall-pressure  fluctuations 

Wall-pressure  fluctuations  beneath  turbulent  boundary  layers  have  been  extensively 
studied.  Eckelmann  [1.13]  provides  a  recent  review  of  the  subject.  For  an  incom¬ 
pressible  flow,  the  fluctuating  pressure  is  related  to  the  fluctuating  velocity  through 
the  Poisson  equation: 


d2p 

dx\ 


dlJidu^  <92 

dxj  dxl  +  dxidxj  UlU^ 


UiUj  ) 


(1.2) 


The  right-hand  side  of  this  equation  represents  two  sources  of  wall-pressure  fluctua¬ 
tions.  The  first  term  is  linear,  and  is  due  to  interaction  of  the  mean  shear  with  the 
turbulence.  The  second  term  is  nonlinear  and  represents  the  interactions  of  the  ve¬ 
locity  fluctuations  with  themselves.  The  fluctuating  pressure  at  the  wall  is  obtained 
as  an  integration  of  the  Poisson  equation  over  the  entire  upper  half-plane.  Since  the 
totality  of  fluid  motions  within  the  boundary  layer  are  reflected  onto  the  surface, 
the  wall-pressure  fluctuations  contain  a  broad  range  of  temporal  and  spatial  scales. 

The  statistical  and  spectral  properties  of  the  fluctuating  wall  pressure  have  been 
addressed  in  detail  by  Schewe  [1.14]  and  Karangelen  [1.15].  They  found  that,  com¬ 
pared  to  a  random  signal  with  Gaussian  probability  density  function,  the  wall  pres¬ 
sure  contains  many  more  large-amplitude  excursions  from  the  mean,  and  that  neg¬ 
ative  large-amplitude  events  occur  more  often  than  positive  events.  This  latter 
finding  is  not  supported  in  the  work  of  Wilczynski  [1.16]  or  in  the  present  research, 
where  positive  and  negative  events  were  found  to  be  equally  likely  (cf.  section  3.3.3). 
Schewe  also  investigated  the  effects  of  transducer  size  on  the  measured  statistics.  He 
showed  that  the  probability  density  function  approaches  a  Gaussian  distribution  as 
the  transducer  size  increases.  This  is  because  an  increasing  number  of  independent 
stochastic  processes  occur  over  the  active  sensing  area. 
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Figure  1.3:  Schematic  of  near-wall  turbulent  motion  and  associated  wall  pressure. 
From  Astolfi  and  Forestier  [1.21]. 


Many  researchers  have  investigated  the  relationship  between  the  wall  pressure 
and  coherent  flow  structures.  The  measurements  suggest  at  least  two  groups  of 
pressure  fluctuations  stemming  from  the  inner-  and  outer-flow  structures.  Kobashi 
and  Ichijo  [1.17]  separated  the  two  signatures  by  performing  correlations  between 
velocity  and  high-  and  low-pass  filtered  pressure  data.  The  low-pass  filtered  pressure 
exhibited  high  correlation  with  the  velocity  across  the  boundary  layer.  The  high- 
pass  filtered  data  was  correlated  with  the  velocity  only  in  the  near- wall  regions. 
Thomas  and  Bull  [1.8]  performed  similar  filtering  and  found  that  the  high-frequency 
pressure  fluctuations  were  associated  with  the  burst-sweep  cycle  of  near-wall  events. 
The  above  results  are  consistent  with  the  findings  of  Farabee  and  Casarella  [1.18], 
who  showed  that  the  high-frequency  part  of  the  wall-pressure  spectrum  scales  on 
wall  variables,  while  the  low  frequency  part  scales  on  outer-flow  variables. 

Conditional  sampling  based  on  large-amplitude  wall-pressure  peaks  and  near¬ 
wall  flow  structures  has  confirmed  a  bidirectional  relationship  between  the  two.  Jo- 
hannson,  Her,  and  Haritonidis  [1.19],  conditioning  on  near-wall  shear  layers,  found 
that  negative  wall-pressure  peaks  were  associated  primarily  with  sweep-type  mo¬ 
tions.  Haritonidis,  Gresko,  and  Breuer  [1.20]  reached  a  similar  conclusion  and  also 
state  that  the  primary  connection  between  the  wall  pressure  and  the  flow  field  is 
through  the  normal  component  of  velocity,  v.  Wilczynski  [1.16]  performed  condi¬ 
tional  sampling  on  the  Reynolds  stress  and  found  that  large-amplitude  wall-pressure 
peaks  are  related  to  Q2  and  Q4  motions.  The  sequence  has  been  clarified  by  Astolfi 
and  Forestier  [1.21]  who  performed  conditional  sampling  and  correlation  on  the  wall 
pressure  and  flow  regions  of  large  du/dy.  Their  schematic  view  of  the  near- wall 
turbulence  and  associated  wall  pressure  is  shown  in  Figure  1.3.  The  results  of  this 
dissertation  are,  to  a  large  extent,  consistent  with  this  model. 

Though  the  link  between  the  wall  pressure  and  near- wall  motions  is  established, 
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the  causal  relationship  between  them  is  not  understood.  In  a  study  of  data  obtained 
from  a  large-eddy  simulation  of  turbulent  channel  flow,  Kim  [1.22]  found  that  the 
ejection  process  at  the  end  of  a  sweeping  motion  was  associated  with  a  localized 
adverse  pressure  gradient.  This  is  in  contrast  with  Thomas  and  Bull  [1.8],  who  con¬ 
cluded  that  the  adverse  pressure  gradients  associated  with  characteristic  pressure 
patterns  are  not  the  direct  cause  of  the  bursting  process.  Furthermore,  Haritoni- 
dis,  et  al.  [1.20]  state  that  “Positive  pressures,  in  general,  appear  to  be  the  result 
of  wall  ward  moving  fluid  being  retarded  by  the  presence  of  the  wall.”  This  issue 
needs  to  be  resolved  before  proper  strategies  for  active  turbulence  control  can  be 
implemented. 

1.2.3  Conditional  sampling 

In  outlining  the  history  of  research  into  turbulent  boundary-layer  structure,  Robin¬ 
son  identified  the  period  1972-1979  as  the  conditional  sampling  era  [1.3,  p.  611]. 
Conditional  sampling  has  played  a  central  role  in  confirming  and  establishing  rela¬ 
tionships  among  the  various  types  of  TBL  structure.  As  pointed  out  in  the  review 
by  Antonia  [1.23],  it  has  helped  to  close  the  gap  between  data  collected  in  an  Eu- 
lerian  frame  of  reference  and  information  in  a  Lagrangian  frame,  as  provided  by 
flow-visualization  experiments.  It  continues  to  figure  prominently  in  the  analysis  of 
virtually  all  digital  experimental  data,  including  the  present  investigation.  There¬ 
fore,  a  brief  review  is  in  order. 

Conditional  sampling  is  generally  performed  by  first  identifying  those  periods  of 
a  signal  that  are  “interesting;”  i.e.,  some  event  detection  criteria.  These  periods 
define  an  indicator  function, 


“  {  0,  otherwise  <13) 

The  times  ti  are  those  for  which  the  event  criteria  is  satisfied.  By  identifying  align¬ 
ment  points  tm  within  the  contiguous  ranges  of  t/,  conditional  or  ensemble  averages 
corresponding  to  the  event-detection  criteria  can  be  formed: 

-i  M 

</M>  =  m  £  /(*»  +  r)  (1.4) 

771=1 

In  this  expression,  M  is  the  number  of  event  realizations  and  r  is  the  delay.  The 
indicator  function  can  be  formed  from  /  or  another  function  at  the  same  or  a  differ¬ 
ent  position  in  space.  When  the  triggering  is  at  a  different  spatial  location,  random 
variations  in  transit  times  of  the  flow  structures  can  significantly  degrade  the  co¬ 
herence  of  the  ensemble  average.  This  problem  is  referred  to  as  phase  jitter,  and 
has  been  addressed  by  Blackwelder  [1.24].  Different  conditional  sampling  meth¬ 
ods  are  distinguished  from  one  another  by  the  criteria  or  algorithms  that  define 
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I.  The  most  common  methods  in  the  literature  include  the  variable-interval  time 
average  (VITA),  the  u-v  quadrant,  and  the  peak-event  detection  methods.  Direct 
comparisons  between  the  various  methods  have  recently  been  done  by  Yuan  and 
Mokhtarzadeh-Dehghan  [1.25]. 

The  VITA  method  computes  the  short-time  variance  of  the  signal  and  identifies 
an  event  as  occurring  when  this  quantity  exceeds  a  certain  threshold.  The  method  is 
usually  applied  to  the  streamwise  velocity,  but  has  been  applied  to  the  wall-normal 
velocity  and  the  Reynolds  stress  as  well.  Variations  on  VITA  include  distinguishing 
the  slope  of  the  function  as  well  as  the  absolute  level.  When  applied  to  velocity 
signals,  the  VITA  technique  is  essentially  a  shear  layer  detector.  Most  of  the  work 
establishing  the  wall-pressure/near- wall  flow  structure  relationship,  cited  in  the  pre¬ 
vious  section,  was  based  on  the  VITA  technique. 

The  u-v  quadrant  method  queries  the  sign  of  u  and  v  to  distinguish  the  four 
quadrants  of  the  u-v  plane.  This  allows  one  to  focus  directly  on  the  turbulence- 
producing  motions  (quadrants  2  and  4).  If  only  strong  events  are  of  interest,  an 
additional  threshold  value  on  the  uv  product  may  also  be  required.  Similar  to  the 
quadrant-plus-level  method,  peak-event  detection  is  perhaps  the  most  straightfor¬ 
ward  means  of  conditioning  ensemble  averages.  Here,  an  event  is  deemed  to  occur 
when  the  signal  amplitude  exceeds  a  given  threshold.  The  work  of  Schewe  [1.14]  and 
Karangelen  [1.15]  made  extensive  use  of  this  technique  in  studying  the  statistics  of 
the  wall  pressure.  Wilczynski  [1.16]  used  both  quadrant  and  peak-event  techniques 
in  proposing  that  Q2-Q4  event  pairs  are  linked  to  wall-pressure  peak  events. 

While  conditional  sampling  has  been  extremely  useful  in  illuminating  coherent 
motions  in  the  boundary  layer,  one  must  exercise  care  in  its  application.  All  of 
the  methods  require  the  input  of  one  or  more  arbitrary  parameters,  in  the  form  of 
averaging  times  or  thresholds.  The  results  obtained  are  highly  dependent  on  the 
values  chosen.  Another  difficulty,  pointed  out  by  Tiederman  [1.26],  occurs  when  a 
detection  scheme  triggers  on  different  portions  of  the  same  flow  structure.  For  exam¬ 
ple,  from  an  Eulerian  viewpoint,  a  single  burst  may  contain  several  ejections.  If  the 
multiple  ejections  are  not  grouped  into  one  event  realization,  skewed  results  will  be 
obtained.  A  more  general  event-detection  scheme  that  requires  a  minimum  number 
of  independent  parameters  is  preferred.  In  addition,  an  approach  that  recognizes 
groupings  of  events,  or  clusters ,  that  are  associated  with  the  passage  of  organized 
motion  is  required. 

1.3  Research  objectives  and  experimental  approach 

This  research  seeks  to  address  two  related  questions  that  are  central  to  the  design 
of  active  turbulence  control  systems: 

1.  What  are  the  detailed  features  of  near- wall  flow  structures  associated  with 
wall-pressure  peak  events  beneath  turbulent  boundary  layers? 
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2.  Are  these  near- wall  structures  coupled  to  and  affected  by  the  outer-layer  dy¬ 
namics? 

To  answer  these  questions,  a  series  of  experiments  was  performed  in  which  si¬ 
multaneous  measurements  of  the  fluctuating  wall-pressure  and  the  streamwise  and 
wall-normal  velocity  were  obtained.  High-quality  databases  of  time-resolved  digital, 
statistical,  and  spectral  data  were  constructed  for  each  of  two  separate  flows.  The 
first  flow  consisted  of  a  fully  developed,  turbulent  boundary  layer  over  a  smooth  wall 
with  zero  pressure  gradient  —  the  equilibrium  or  canonical  case.  The  second  flow  was 
a  nonequilibrium  or  disturbed  turbulent  boundary  layer,  in  which  the  disturbances 
were  primarily  confined  to  the  outer  regions  of  the  boundary  layer.  Specifically,  the 
disturbed  flow  was  that  formed  downstream  of  reattachment  of  the  flow  over  an 
aft-facing  step. 

Once  the  data  were  collected,  the  following  tasks  were  necessary  to  address  the 
stated  objectives: 

•  Confirm  that  wall-pressure  peak  events  are  footprints  of  near-wall  organized 
motion  in  the  equilibrium  layer  using  traditional  conditional  sampling  tech¬ 
niques. 

•  Develop  appropriate  filtering  to  extract  the  near- wall  motions  associated  with 
wall-pressure  peak  events  in  the  disturbed  layer. 

•  Formulate  a  new  conditional  sampling  approach  that  recognizes  cluster  group¬ 
ings  of  events  in  order  to  test  the  robustness  of  the  results. 

•  Document  the  spatial  extent  of  the  observed  near- wall,  turbulence-producing 
motions. 

The  two  turbulent  boundary  layers  studied  will  now  be  described. 

1.3.1  Equilibrium  turbulent  boundary  layer 

The  criteria  for  the  existence  of  a  fully  developed,  equilibrium  turbulent  boundary 
layer  are  given  by  Hussain  [1.27]: 

•  A  mean  streamwise  velocity  profile  with  conventional  shape  factor  H  «  1.4, 
Clauser  shape  factor  G  ~  6.3,  and  shear-to-freestream  velocity  ratio  ur/UP  = 
f(Reg)  »  0.04. 

•  A  mean  streamwise  velocity  profile  with  a  log-law  region  of  adequate  extent 
based  on  the  flow  Reynolds  number. 

•  A  wake  strength  that  is  appropriate  for  the  flow  Reynolds  number. 

•  A  peak  in  the  turbulence  intensity  profile  of  u/uT  =  2.5  ±  10 percent. 
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Figure  1.4:  Features  of  the  flow  over  a  backward-facing  step.  From  Bradshaw  and 
Wong  [1.30]. 


•  A  location  of  the  peak  turbulence  intensity  at  y+  =  15. 

•  A  monotonic  decrease  in  u  to  the  freestream  value. 

•  A  broadband,  continuous  u  spectrum  that  contains  an  inertial  subrange. 

In  addition,  the  wall  pressure  rms  level  should  be  appropriate  for  the  flow  Reynolds 
number  and  have  a  smooth  spectrum  with  no  apparent  discontinuities.  These  char¬ 
acteristics  are  confirmed  for  the  present  study  in  Chapter  3. 


1.3.2  Disturbed  turbulent  boundary  layer 

The  disturbed  flow  chosen  for  this  study  approximated  the  flow  downstream  of  reat¬ 
tachment  of  flow  over  a  backward-facing  step.  The  selection  was  based  on  previous 
studies  of  step  flows  conducted  in  the  same  facility  [1.28,  1.29].  The  flow  was  ob¬ 
tained  by  installing  a  ramp  on  the  equilibrium  TBL  test  surface.  The  ramp  had  an 
inclination  of  1.7  deg  and  terminated  in  a  backward-facing  step  of  height  1.5  cm. 
One  advantage  is  a  fixed  point  of  separation.  While  not  strictly  the  same  as  flow 
over  a  step,  the  ramp  flow  provided  the  desired  disturbed  TBL  for  comparison  with 
the  equilibrium  case. 

The  characteristic  features  of  flow  over  a  backward-facing  step  are  illustrated 
schematically  in  Figure  1.4,  taken  from  Bradshaw  and  Wong  [1.30].  The  following 
discussion  is  based  on  the  detailed  measurements  of  Chandrsuda  and  Bradshaw  [1.31] 
and  Farabee  and  Casarella  [1.28,  1.29].  The  flowfield  can  be  described  in  terms  of 
overlapping  zones.  In  the  first  zone,  an  initial  or  inflow  boundary  layer  flows  along 
the  top  of  the  step.  At  the  edge  of  the  step,  the  boundary  layer  separates  and  a  free 
shear  layer  forms  over  a  second  zone  of  recirculating  flow.  The  shear  layer  grows  and 
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makes  contact  with  the  wall  over  a  third  zone  of  reattachment,  which  contains  high 
levels  of  turbulence  activity.  Immediately  downstream  of  reattachment,  a  boundary 
layer  develops  which  is  highly  energized  and  has  the  characteristics  of  an  attached 
free  shear  layer.  This  developing  layer  has  a  log  law  region  even  very  close  to 
reattachment.  As  the  flow  moves  downstream,  the  boundary  layer  relaxes,  eventually 
becoming  an  equilibrium  TBL. 

The  rms  wall  pressure  peaks  just  upstream  of  reattachment,  decreasing  mono- 
tonically  thereafter.  The  magnitude  of  the  peak  is  three  to  five  times  greater  than 
for  a  flat  plate,  depending  on  the  thickness  of  the  inflow  boundary  layer.  The  wall- 
pressure  spectrum  exhibits  increased  low-frequency  components  compared  to  the 
flat  plate  case.  Similar  to  the  rms  values,  the  highest  spectral  levels  occur  at  reat¬ 
tachment,  followed  by  a  rapid  decrease  in  the  low-frequency  components  as  the  flow 
moves  downstream. 


1.4  Outline  of  dissertation 

The  following  chapters  present  the  details  of  the  experimental  measurements  and 
the  signal  processing  methods  used  to  extract  the  desired  information.  Chapter  2 
documents  the  experimental  setup  and  methods  used  to  obtain  simultaneous  wall 
pressure  and  flowfield  velocity  measurements  across  the  equilibrium  and  disturbed 
turbulent  boundary  layers.  The  long-time  averaged  properties  of  the  two  flows  are 
presented  in  Chapter  3.  It  will  be  shown  that  the  disturbed  TBL  differs  from  the 
equilibrium  TBL  by  the  addition  of  large-scale,  low-frequency  disturbances  in  the 
outer  regions  of  the  boundary  layer,  but  that  the  inner  layers  appear  to  be  substan¬ 
tially  the  same.  The  hypothesis  is  made  that  the  inner  and  outer  layer  portions  of 
turbulent  boundary  layers  are  uncoupled. 

To  test  this  hypothesis,  Chapter  4  implements  a  wavelet  filtering  scheme  to  iso¬ 
late  the  common  near- wall  features  of  both  flows.  The  analysis  of  the  data  is  further 
enhanced  using  the  idea  of  event  clusters,  introduced  in  Chapter  5.  Two  new  analysis 
methods,  Localized  Variance  and  Localized  Windowed  Peak  Detection,  will  be  used 
to  form  correlations  and  conditionally  sampled  ensemble  averages  of  the  flowfield 
structures  of  interest.  Correlations  will  be  used  to  compute  convection  velocities  of 
the  wall-pressure  events,  and  conditional  sampling  results  will  confirm  the  bidirec¬ 
tional  relationship  between  positive  wall-pressure  peaks  and  accelerating  u-velocity 
or  second-quadrant  uv  motions.  Finally,  Chapter  6  assembles  the  results  into  a  com¬ 
posite  picture  of  the  near-wall  motions,  allowing  the  fundamental  questions  of  the 
research  to  be  addressed. 
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Chapter  2 

Experimental  setup  and  methods 


2.1  Introduction 

The  experiments  were  performed  in  the  Catholic  University  Low  Noise  Flow  Facility 
(LNFF).  Two  flows  were  studied:  an  equilibrium  turbulent  boundary  layer  (ETBL) 
and  a  disturbed  turbulent  boundary  layer  (DTBL).  The  instrumentation  used  in¬ 
cluded  sensors  to  monitor  the  state  of  the  wind  tunnel,  two  microphones  for  the 
measurement  of  the  fluctuating  wall  pressure,  and  single-component  hot  wire  and 
two-component  hot  film  anemometers  for  the  measurement  of  fluctuating  streamwise 
and  wall-normal  velocities.  In  this  chapter,  the  facility,  setup,  measurement  systems 
and  calibration  techniques,  and  data-acquisition  hardware  will  be  described,  along 
with  a  description  of  the  testing  procedures  and  limitations.  A  detailed  measurement 
uncertainty  analysis  is  presented  in  Appendix  B. 

2.2  Facility 

The  LNFF  was  constructed  specifically  for  making  acoustic  measurements  of  wall 
pressure  beneath  turbulent  flows.  A  complete  description  can  be  found  in  [2.1].  A 
schematic  is  shown  in  Figure  2.1.  The  inlet  is  covered  by  honeycomb  and  multiple 
screens  for  turbulence  management.  This  is  followed  by  a  16:1  contraction  section 
designed  to  provide  a  gradual  static  pressure  drop  and  to  avoid  flow  separation 
along  the  tunnel  walls.  The  test  section  is  2.4  m  long  and  0.6  m  x  0.6  m  in  cross 
section.  The  tunnel  is  powered  by  a  25  horse  power,  variable  speed  dc  motor  driving 
a  centrifugal  blower  in  an  open-return  arrangement.  The  motor/blower  assembly  is 
located  in  a  separate  room  and  is  isolated  by  acoustic  mufflers  to  minimize  facility 
vibrations  and  acoustic  noise  levels. 

The  tunnel  speed  can  be  varied  from  approximately  3  to  30  m/s.  At  a  nomi¬ 
nal  speed  of  16  m/s,  the  freestream  turbulence  level  ?w/t/e  is  approximately  0.2 
percent.  As  documented  in  earlier  studies  [2. 1-2.3]  and  confirmed  in  the  present 
research,  the  acoustic  performance  is  excellent  above  approximately  50  Hz.  Below 
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Figure  2.1:  The  Catholic  University  Low  Noise  Flow  Facility. 

this  level,  the  wall  pressure  is  contaminated  by  blower-generated  acoustic  standing 
waves.  As  discussed  in  a  later  section,  this  noise  is  well  below  the  frequency  range 
of  interest  and  is  removed  by  conventional  digital  filtering. 

The  boundary  layers  studied  in  this  investigation  were  formed  on  one  wall  of 
the  wind  tunnel.  The  test  wall  was  constructed  of  plywood  covered  with  formica. 
A  turbulent  layer  was  ensured  through  the  use  of  two  0.81-mm  diameter  trip  wires 
epoxied  to  the  wall  0.3  and  0.9  m  upstream  of  the  test  section.  The  measurements 
were  obtained  on  the  tunnel  center  line  at  approximately  1.4  m  downstream  of  the 
test  section  entrance  («2100  trip  heights  from  the  second  trip).  At  this  location, 
the  boundary  layer  is  approximately  3  cm  thick  for  a  tunnel  speed  of  16  m/s.  The 
equilibrium  state  of  the  boundary  layer  will  be  demonstrated  in  the  next  chapter. 

The  disturbed  turbulent  boundary  layer  approximated  the  flow  over  an  aft-facing 
step.  This  was  accomplished  by  installing  a  formica  ramp  of  length  0.51  m  on  the 
tunnel  wall  ending  0.13  m  upstream  of  the  measurement  station.  The  step  height, 
h,  was  1.5  cm.  The  ramp  spanned  the  width  of  the  tunnel,  and  the  leading  edge  was 
faired  into  the  tunnel  wall  using  modeling  clay.  The  characteristics  of  the  resulting 
boundary  layers  are  presented  in  Chapter  3. 

2.3  Facility  instrumentation 

The  tunnel  conditions  were  determined  through  the  use  of  a  test  cell  Pitot-static 
probe.  The  Pitot-static  difference  was  measured  directly  with  an  Edwards  High 
Vacuum  (formerly  Datametrics)  model  570D-10T-2C2-V1X  0-10  torr  differential 
pressure  transducer  and  model  1174  electronic  manometer  and  power  supply.  The 
transducer  system  zero,  full  scale,  and  sensitivity  were  set  daily  using  internal  cal¬ 
ibration  circuitry.  Atmospheric  pressure  and  temperature  were  measured  using  a 
Welch  Scientific  fixed  cistern  mercury  barometer  and  a  mercury-in-glass  thermome- 
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ter.  The  density  was  computed  using  the  perfect  gas  equation  of  state.  Tunnel 
velocity  was  computed  using  the  measured  Pitot-static  difference  and  the  calculated 
density  through  Bernoulli’s  equation. 

As  outlined  in  Section  2.5,  the  output  of  a  hot-wire  anemometer  is  a  function 
of  the  static  temperature  of  the  fluid.  To  compensate  for  temperature  changes  over 
the  course  of  the  experiments,  the  static  temperature  of  the  tunnel  was  monitored 
with  a  chromel-alumel  (type  K)  thermocouple  and  Omega  type  DP11  thermocouple 
reference  junction/monitor.  The  thermocouple  was  inserted  into  the  flow  through 
the  test  wall  approximately  0.86  m  downstream  of  the  measurement  station. 

2.4  Wall  pressure  sensors 

Wall-pressure  measurements  were  made  using  two  microphones,  identified  as  pi  and 
P2-  Microphone  p\  was  mounted  at  the  primary  measurement  station  beneath  the 
hot-wire  probe,  while  microphone  p2  was  mounted  1.27  cm  upstream  of  px.  The 
microphones  were  Briiel  and  Kjaer  (B&K)  model  4138  |  in.  condenser  microphones. 
Each  microphone  was  connected  to  a  B&K  type  2633  preamplifier  using  a  type 
UA160  |-in.  to  |-in.  adapter.  Both  microphones  were  powered  by  a  B&K  type 
2807  two-channel  power  supply.  The  signals  were  then  amplified  an  additional  40  dB 
(disturbed  TBL)  or  50  dB  (equilibrium  TBL)  using  an  Ithaco  model  451  amplifier 
with  built-in  high-pass  filter  set  at  a  cutoff  frequency  of  1  Hz.  The  phase  delay  from 
this  filter  was  less  than  1  deg  for  frequencies  greater  than  100  Hz.  Both  the  B&K 
and  the  Ithaco  amplifiers  were  of  the  inverting  type;  their  series  connection  kept  the 
output  signal  in  phase  with  the  fluctuating  pressures.  Figure  2.2  shows  a  schematic 
of  the  microphone  system  in  relation  to  the  hot-wire  and  data-acquisition  systems. 
Only  one  microphone  is  illustrated  for  clarity. 

The  microphone/preamplifier  assemblies  were  inserted  into  plastic  holders,  which 
were  in  turn  mounted  in  a  track  attached  to  the  back  of  the  test  wall.  This  allowed 
the  microphone  to  be  positioned  at  the  desired  measurement  location.  The  micro¬ 
phone  was  pushed  into  the  hole  until  the  pinhole  cap  extended  past  the  test  surface. 
Using  the  sharp  edge  of  a  razor  blade,  the  cap  was  then  pushed  back  flush  with  the 
wall.  Modeling  clay  was  then  carefully  applied  to  the  seam  between  the  wall  and 
the  cap,  and  feathered  with  the  razor  blade.  This  ensured  a  smooth  surface  with 
minimal  disturbances. 

2.4.1  Microphone  calibration 

The  microphones  were  calibrated  against  a  B&K  type  4220  pistonphone.  The  pis- 
tonphone  produced  a  known  sound-pressure  level  (nominally  150  dB  re  1  pPa)  at 
a  frequency  of  250  Hz.  Individual  calibration  curves  provided  by  the  manufacturer 
demonstrate  that  each  microphone  has  a  flat  frequency  response  within  ±1  dB  out 
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to  50  KHz.  The  microphone  sensitivity  was  then  calculated  as 

8  =  SPL/Erms  (2.1) 

where  s  is  the  microphone  sensitivity  in  Pa/V,  SPL  is  the  pistonphone  sound  pressure 
level,  and  E ^  is  the  rms  output  voltage  measured  at  the  data  acquisition  system. 
For  the  disturbed  TBL  experiments  with  the  Ithaco  amplifiers  set  to  40  dB,  a  typical 
microphone  sensitivity  of  12.7  Pa/V  was  obtained. 

2.4.2  Microphone  resolution 

To  minimize  spatial  averaging  of  the  fluctuating  wall  pressure,  the  microphone’s 
standard  protective  screen  was  replaced  with  a  pinhole  cap.  The  cap  had  a  0.79-mm 
diameter  hole  drilled  in  the  center.  This  resulted  in  nondimensional  sensing  diame¬ 
ters  d+  =  duT/v  of  32  and  23  for  the  equilibrium  and  disturbed  turbulent  boundary 
layers,  respectively.  The  pinhole  cap  and  microphone  assembly  form  a  Helmholz 
resonator,  which  limits  the  useful  frequency  range  of  the  microphone  data.  Previous 
researchers  [2.1,  2.4]  have  established  that  the  Helmholz  resonance  frequency  for  the 
current  sensor  configuration  is  approximately  25  KHz.  This  was  confirmed  in  the 
present  research  as  a  small  spike  in  the  wall-pressure  spectrum.  As  discussed  in  the 
Chapter  3,  the  Helmholz  frequency  is  well  beyond  both  the  energy-containing  range 
of  the  wall  pressure  and  the  digital  sampling  Nyquist  frequency  of  16.3  KHz. 

2.5  Hot-wire  anemometers 

Velocity  measurements  were  obtained  using  a  constant-temperature  anemometry 
(CTA)  system  comprised  of  the  following  components: 

•  TSI  Inc.  type  1261  A-T1.5  miniature  boundary-layer  probe  or  type  1249  A-10 
miniature  “X”  probe 

•  TSI  1155  probe  support  on  a  motorized  traverse 

•  Dantec  (formerly  DISA)  type  55M10  CTA  standard  bridges  in  type  55M01 
main  units 

•  Dantec  type  55M25  linearizers 

•  TSI  type  1015C  correlator 

The  anemometer  system  is  illustrated  schematically  in  Figure  2.2. 

The  type  1261  A-T1.5  probe  consisted  of  a  single  tungsten  wire  of  length  1.52  mm, 
3.8  gm  in  diameter,  with  an  active  sensing  length  of  1.27  mm.  The  wire  was  oriented 
normal  to  the  streamwise  velocity  in  the  z-direction.  Curved  wire  supports  allowed 
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the  wire  to  be  positioned  close  to  a  solid  boundary.  The  type  1249  A- 10  “X”  probe 
consisted  of  two  ceramic  “wires”  on  which  were  deposited  a  platinum  film.  The 
wires  were  of  length  1.27  mm  and  diameter  25  pm,  with  active  film  sensing  lengths 
of  0.51  mm.  The  wires  were  oriented  in  the  x-y  plane,  at  ±45  deg  to  the  x-axis.  The 
two  wires  were  separated  in  the  z-direction  by  1.0  mm.  The  probe  incorporated  a 
90-deg  elbow  for  making  measurements  close  to  the  wall. 

2.5.1  Hot-wire  calibration 

For  digital  measurements,  the  bridge  outputs  were  input  directly  into  the  digital  data 
acquisition  system,  which  is  described  in  Section  2.6.1.  The  single-wire  voltages  were 
converted  to  velocities  in  a  conventional  way  [2.5,  pp.  15-17],  based  on  a  form  of 
King’s  law: 

E2  =  A  +  BUn  (2.2) 

where  E  is  the  anemometer  output  voltage,  U  is  the  velocity  in  m/s,  and  A,  B,  and 
n  are  constants.  The  exponent,  n,  was  allowed  to  vary  from  0.1  to  0.7  in  increments 
of  0.01,  and  a  least-squares  method  used  to  determine  A  and  B.  The  value  of  n 
that  resulted  in  the  smallest  standard  error  fixed  the  calibration.  Eight  voltage- 
velocity  pairs  were  used  for  the  calibration.  The  tunnel  velocity  was  determined 
from  the  Pitot-static  probe.  Typical  calibration  data  and  the  resulting  fit  are  shown 
in  Figure  2.3. 

The  X-film  probe  was  calibrated  according  to  the  method  of  Camelio  [2.6]  as  pre¬ 
sented  in  Wilczynski  [2.4].  The  two  velocity  components,  U  and  V,  were  represented 
by  fourth-order,  two-dimensional  polynomials  in  the  two  wire  voltages,  E\  and  E2\ 
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Calibration  velocities  U  and  V  were  obtained  by  rotating  the  probe  to  seven  different 
angles,  </>,  with  respect  to  the  tunnel  velocity,  Q,  so  that  U  =  Q  cos  <j>  and  V  = 
Qsin(f>.  Seven  different  angles  were  used:  <j>  =  —36,  —18,  —9,  0,  9,  18,  and  36  deg. 
At  each  angle,  data  were  recorded  at  eight  tunnel  speeds,  resulting  in  a  total  of  56 
calibration  points.  The  coefficients  Gy  and  6y  were  then  computed  using  a  least- 
squares  regression  method.  Typical  data  and  the  resulting  calibration  surfaces  are 
shown  in  Figures  2.4  and  2.5. 

More  common  methods  of  calibrating  X-probes  found  in  the  literature  involve 
fitting  (Q,  (f>)  to  (Ei,E2),  then  solving  for  U  and  V.  Various  curve  fitting  and  table 
look-up  implementations  can  be  found.  The  present  method,  however,  was  felt  to 
be  a  more  direct  approach,  without  sacrificing  accuracy. 
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The  calibration  methods  outlined  above  assume  that  the  flow  being  measured 
is  isothermal  and  that  the  temperature  remains  constant.  In  the  course  of  the 
experiments,  changes  of  up  to  3  deg  C  were  observed  between  calibration  and  mea¬ 
surement.  To  account  for  these  changes,  a  temperature  compensation  procedure  was 
implemented  as  given  in  [2.5,  pp.  194-198]: 

E2 It W  =  (2.5) 

Here,  E  is  the  wire  voltage,  T  is  fluid  temperature  at  the  time  of  the  measurement, 
Tw  is  the  wire  operating  temperature,  and  Tref  is  the  reference  temperature.  The 
reference  temperature  was  chosen  to  be  the  mean  temperature  observed  over  the 
course  of  collecting  the  calibration  data.  The  correction  was  applied  to  the  voltages 
prior  to  using  them  to  compute  calibration  coefficients  or  velocities.  This  correction 
scheme  is  valid  for  temperature  changes  of  a  few  degrees  or  less.  For  larger  variations, 
more  involved  corrections  must  be  made. 

Spectral  measurements  (Section  2.6.2)  require  a  voltage  signal  directly  propor¬ 
tional  to  the  velocity.  For  the  single-wire  data,  the  anemometer  output  was  fed 
through  the  55M25  linearizer  before  being  input  to  the  spectrum  analyzer.  The 
linearization  was  adjusted  at  flow  speeds  near  the  minimum  and  slightly  above  the 
maximum  velocity  measured.  For  the  X-probe,  the  wire  voltages  were  linearized 
corresponding  to  the  flow  components  ±45  deg  to  the  x-direction.  The  linearized 
signals  were  then  fed  into  the  correlator,  which  added  and  subtracted  them  to  pro- 
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vide  signals  proportional  to  U  and  V.  The  X-probe  case  is  illustrated  in  Figure  2.2. 

Solid  boundaries  in  close  proximity  to  a  hot-wire  probe  can  act  as  a  heat  sink, 
conducting  heat  away  from  the  probe.  This  effect  will  appear  as  an  increase  in 
velocity,  since  more  voltage  is  required  to  replace  the  lost  heat.  A  survey  of  the 
literature  on  this  effect  can  be  found  in  [2.5].  For  nonconducting  walls,  such  as 
the  formica  covered  plywood  of  the  LNFF,  corrections  need  not  be  made  for  mea¬ 
surements  where  y+  =  yuTjv  >  2.5,  and  then  only  for  very  low  flow  speeds.  All 
measurements  in  this  investigation  were  obtained  at  a  sufficient  distance  from  the 
wall  and  at  a  high  enough  flow  speed  such  that  corrections  for  wall  conduction  effects 
were  not  necessary. 

2.5.2  Hot-wire  resolution 

The  spatial  resolution  of  the  velocity  measurements  was  limited  by  the  active  sens¬ 
ing  length  of  the  single-wire  probe  and  by  the  separation  distance  of  the  X-probe. 
Expressed  in  wall  units  of  l+  —  lur/u,  the  resolutions  for  the  equilibrium  TBL  were 
50  for  the  single  wire  and  41  for  the  X-probe.  For  the  disturbed  TBL,  the  resolutions 
were  26  and  29  wall  units  for  the  single  and  X-probes,  respectively. 

2.5.3  Probe  displacement  measurements 

The  hot-wire  probe  support  was  mounted  on  a  three  degree-of-freedom  (3-DOF) 
traverse  inside  of  the  wind  tunnel.  The  traverse  could  be  controlled  manually  as 
well  as  remotely  via  stepper  motors.  In  practice,  the  probe  was  positioned  manually 
over  microphone  pi,  at  approximately  1-mm  downstream  of  the  pinhole. 

The  probe  was  initially  brought  close  to  the  wall  while  being  observed  through 
an  Edmund  Scientific  6— 8x  close-focus  telescope.  The  telescope  was  equipped  with 
a  combination  radius/angle  reticle  and  was  mounted  on  a  3-DOF,  micrometer- 
positioned  platform  manufactured  by  the  Line  Tool  Company.  The  platform  was 
placed  on  the  top  of  the  tunnel  such  that  the  telescope  was  located  at  the  same 
x-location  as  the  hot-wire  probe,  at  y  «  0,  and  sighted  along  the  positive  z-axis. 
From  this  shallow  viewing  angle,  with  the  telescope  focused  on  the  wire  support 
prongs,  the  reflection  of  the  probe  was  clearly  visible  in  the  test  wall.  The  probe 
could  therefore  be  brought  as  close  to  the  wall  as  possible  without  actually  making 
contact. 

Once  positioned,  the  initial  y-displacement  could  be  determined  with  the  tele¬ 
scope  platform  positioning  micrometers  by  taking  readings  with  the  reticle  cross 
hairs  on  the  wire  and  then  on  the  wall.  This  was  repeated  twice  and  the  results 
averaged.  Dust  made  the  wall  location  clear  midway  between  the  probe  image  and 
its  reflection.  The  accuracy  of  the  micrometer  readings  was  verified  by  sighting 
the  telescope  on  a  precision  scale  fixed  to  the  tunnel  wall  with  modeling  clay.  The 
position  of  the  hot  wire  at  subsequent  measurement  locations  was  monitored  with 
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an  Ono  Sokki  DG-357  digital  gauge  attached  to  the  traverse  assembly.  Both  the 
telescope  micrometers  and  the  potentiometer  system  had  a  resolution  of  26  gm. 

2.6  Data  acquisition  and  processing  systems 

Data  acquired  in  this  research  fell  into  two  categories:  digital  and  spectral.  This 
section  will  detail  the  two  systems,  which  are  illustrated  in  Figure  2.2. 

2.6.1  Digital  data  system 

Digital  data  were  acquired  on  a  system  based  on  the  VMEbus  Extensions  for  Instru¬ 
mentation  (VXIbus  or  VXI)  specification.  The  VXI  system  uses  a  mainframe  chassis 
with  slots  to  hold  modular  instruments  on  plug-in  boards.  The  specific  system  was 
manufactured  by  the  Hewlett  Packard  Company.  An  HP  E1401A  C-size  mainframe 
provided  the  VXI  backplane,  cooling,  power,  and  current  requirements.  Digitization 
was  performed  by  an  HP  E1431A  25.6  KHz  Eight  Channel  VXI  Input  Module.  This 
module  was  chosen  for  its  high  sampling  rates,  built-in  signal  conditioning,  and  high 
resolution,  16-bit,  analog-to-digital  converter. 

The  module  featured  simultaneous  sampling  on  up  to  eight  differential  input 
channels.  Sampling  rates  were  selectable  in  powers  of  two  from  1  to  65536  Hz,  with 
alias  protection  provided  internally  by  analog  and  digital  antialias  filters.  The  inter¬ 
nal  digital  filters  used  2.56  x  oversampling,  resulting  in  a  maximum  alias-protected 
input  bandwidth  of  25.6  KHz.  System  gains  allowed  input  voltage  ranges  of  5  mV 
to  10  V  peak  and  were  selectable  on  an  individual  channel  basis.  In  addition,  buck¬ 
ing  amplifiers  provided  dc  offset  compensation  individually  for  each  channel.  The 
bucking  amplifier  subtracted  a  constant  voltage  from  the  input  signal,  allowing  the 
remainder  to  be  amplified  and  digitized,  thus  maximizing  the  resolution  of  the  low- 
amplitude  fluctuations.  This  feature  was  exploited  in  the  velocity  measurements, 
where  a  typical  input  consisted  of  mV-level  fluctuations  superimposed  on  a  4-5  V 
signal.  Conventional  high-pass  filtering  is  not  acceptable  in  this  situation  due  to 
the  nonlinear  voltage-velocity  relationship  of  the  constant-temperature  anemometer 
system. 

The  system  was  controlled  by  an  HP  model  735  workstation  via  Multi-system 
extension  Interface  bus  (MXIbus),  using  an  HP  E1482B  VXI-to-MXI  interface  mod¬ 
ule  in  the  mainframe  and  an  HP  E1439I  EISA/ISA-to-MXI  interface  card  in  the 
workstation.  The  MXIbus  is  a  general-purpose,  32-bit  multimaster  system  bus  on 
a  cable.  It  uses  a  hardware  memory-mapped  communication  scheme  to  minimize 
software  overhead  for  very  high  speed.  Data  were  acquired  from  the  E1431A  direct 
to  the  HP  735  main  memory.  The  maximum  theoretical  data  rate  for  MXIbus  is 
20  Mbytes/s  [2.7].  The  use  of  MXI  for  system  control  allowed  realtime  acquisition  of 
all  eight  channels  of  the  E1431A  sampled  at  32768  samples/s  without  the  need  for 
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local  memory  modules  in  the  VXI  mainframe.  A  maximum  of  four  channels  was  used 
simultaneously  to  acquire  data  from  the  two  microphones  and  the  two-component 
hot  wire.  The  tunnel  temperature  and  the  hot-wire  probe  position  were  acquired 
and  recorded  manually. 

The  HP  735  was  a  general  purpose  Unix  workstation.  It  was  equipped  with 
64  Mbytes  of  RAM,  an  HP  PA-RISC  central  processing  unit  clocked  at  90  MHz, 
and  approximately  6  Gbytes  of  hard  disc  storage.  The  software  to  control  the  VXI 
system  was  custom  developed  in-house  [2.8].  It  provided  complete  control  of  the 
E1431A  setup  and  acted  as  the  external  trigger.  The  data  were  received  into  the 
HP  735  memory  and  then  written  to  disc.  One  file  was  created  per  active  channel. 
The  raw  data  were  stored  in  binary  format  and  included  the  channel  sample  rate, 
digital  count  to  voltage  conversion  factor,  and  dc  offset  in  the  record. 

2.6.2  Spectral  data  system 

Spectral  data  were  acquired  with  an  Hewlett  Packard  HP  3562A  Dynamic  Signal 
Analyzer.  The  HP  3562A  is  a  two-channel,  fast  Fourier  transform  analyzer  with 
internal  antialias  filters  and  a  100-KHz  input  bandwidth.  The  analyzer  was  used  to 
compute  auto  spectra,  cross  spectra,  phase  response,  and  coherence  functions. 

A  feature  of  the  HP  3562A  utilized  in  this  study  was  its  ability  to  compute  spec¬ 
tral  measurements  in  a  logarithmic  resolution  mode.  In  this  mode,  the  frequency 
span  of  interest  was  represented  by  80  discrete  frequencies  per  decade,  logarithmi¬ 
cally  spaced.  When  plotted  on  a  log-frequency  axis,  spectra  computed  in  this  mode 
had  equal  resolution  in  the  lowest  and  the  highest  decades.  This  proved  to  be  useful, 
since  the  turbulent  flow  fields  of  this  study  had  significant  spectral  features  spanning 
a  frequency  range  of  10  Hz  to  approximately  20  KHz. 

The  analyzer  was  controlled  by  an  IBM-compatible  486  personal  computer  equip¬ 
ped  with  an  HP  82335  HP-IB  (IEEE  488  or  GPIB)  interface  card.  Custom  software 
was  developed  to  initialize  the  analyzer,  trigger  the  measurements,  and  write  the 
results  to  external  disc  drives  on  the  GPIB  bus.  As  illustrated  in  Figure  2.2,  the 
data  were  then  available  to  the  PC  via  the  GPIB  bus  and  to  the  HP  735  via  Ethernet 
(IEEE  802.3)  local  area  network.  The  measurements  were  made  in  the  logarithmic 
resolution  mode  over  a  frequency  span  of  10  Hz  to  100  KHz.  The  final  results  were 
the  ensemble  average  of  fifty  individual  spectra,  using  50  percent  overlap  and  the 
rms  voltage  as  the  basic  amplitude  unit. 

2.7  Test  conduct  and  data  processing 

This  section  describes  the  experimental  procedures  for  collecting  the  raw  data  and 
the  first-order  processing  done  to  create  data  bases  needed  for  further  analyses. 
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2.7.1  Test  conduct 

Data  were  collected  for  two  different  flows:  an  ETBL  and  a  DTBL.  For  each  of 
these  flows,  a  single-wire  survey  and  several  X-wire  surveys  were  performed.  For 
each  single  wire  survey,  the  streamwise  velocity  U,  and  the  wall  pressure  at  two 
x-locations  p\  and  p2,  were  measured  at  each  of  22  or  23  y-locations  within  and 
outside  of  the  boundary  layer.  For  each  X-wire  survey,  the  wall-normal  component 
of  velocity  V  was  measured  as  well.  Several  X-wire  surveys  were  required  to  obtain 
the  desired  cross  spectral  data  pairs:  U-V,  U-pi,  and  V-p\. 

The  initial  task  was  the  calibration  of  the  microphones.  Each  microphone  was 
coupled  to  the  pistonphone  using  a  plastic  coupler,  and  its  output  recorded  on  the 
VXI  system  in  the  same  channel  to  be  used  for  the  survey.  The  rms  voltage  com¬ 
puted  from  the  calibration  data  was  checked  against  a  voltmeter  for  consistency. 
The  microphones  were  then  installed  in  the  test  surface.  Next,  the  Pitot-static  pres¬ 
sure  transducer  was  adjusted  for  zero,  full  scale,  and  sensitivity  according  to  the 
manufacturer’s  procedures. 

At  this  point,  the  single-  or  X-wire  probe  was  installed  in  the  traverse  and  the 
anemometer  bridge  was  set  up  for  a  70  percent  overheat  ratio.  The  tunnel  was  then 
set  to  a  speed  approximately  equal  to  Ue,  and  the  bridge  balanced  for  frequency 
response  using  the  internal  square  wave  circuitry.  The  useful  frequency  ranges  for 
both  types  of  probe  were  found  to  be  at  least  47.6  KHz.  Analog  linearization  of  the 
anemometer  was  then  checked  at  two  tunnel  speeds:  ~  0.2 Ue  and  1.0 Ue.  Following 
linearization,  the  hot-wire  calibration  was  performed.  If  the  flow-disturbance  ramp 
was  installed,  the  hot-wire  probe  was  traversed  well  upstream  of  the  ramp  where  the 
flow  angles  were  undisturbed.  For  the  X-wire  probe,  the  final  calibration  angle  was 
always  0  deg  so  that  the  probe  was  not  disturbed  again  prior  to  the  survey. 

With  the  preparations  complete,  the  data  collection  was  begun.  The  probe  was 
positioned  close  to  the  wall  as  described  in  Section  2.5.3,  and  the  initial  displacement 
recorded.  The  tunnel  was  then  started.  Flow  speeds  were  set  according  to  motor 
speed,  thus  some  variation  in  the  actual  velocities  between  surveys  was  observed. 
The  spectral  analyzer  was  then  triggered  to  ensemble  average  50  individual  spectra. 
Concurrent  with  the  spectral  measurements,  the  digital  system  was  triggered  to 
record  two  data  sets: 

1.  Statistics  data,  consisting  of  140  s  of  the  anemometer  channel(s)  at  500  sam¬ 
ples/s  per  channel 

2.  Time  series  data,  consisting  of  10  s  of  all  channels  at  32768  samples/s  per 
channel. 

Both  the  digital  and  the  spectral  systems  required  approximately  2.5  min  to  com¬ 
plete  their  measurements.  Once  both  systems  were  finished,  the  traverse  was  moved 
outward  and  data  collected  at  the  new  location.  This  measure-move-measure  pro¬ 
cess  was  repeated  for  22  or  23  locations  throughout  the  boundary  layer.  When  the 
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Figure  2.6:  Pressure  measurement  locations  relative  to  the  aft-facing  step. 

survey  was  complete,  the  probe  was  returned  to  the  near-wall  position  for  collection 
of  the  remaining  cross-spectral  combinations. 

The  U  power  spectra  and  U-pi  cross  spectra  obtained  with  the  single- wire  probe 
were  compared  with  those  obtained  with  the  X-probe  for  consistency.  Differences 
in  the  power  spectra  were  less  than  1  dB  over  the  entire  100-KHz  frequency  span, 
with  a  similar  agreement  in  the  coherence  data.  Therefore,  the  U-pi  cross  spectra 
were  obtained  with  the  single-wire  probe,  since  this  probe  could  be  positioned  much 
closer  to  the  wall  than  the  X-probe.  Cross-spectral  data  were  also  obtained  for  the 
P1-P2  input  pair  with  the  hot-wire  probe  at  its  closest  and  farthest  measurement 
locations  from  the  wall. 

Figure  2.6  shows  the  locations  relative  to  the  aft-facing  step  where  pressure  mea¬ 
surements  were  made.  The  stations  labeled  p\  and  P2  were  the  primary  stations, 
where  pressure  measurements  in  conjunction  with  hot-wire  anemometry  were  si¬ 
multaneously  obtained.  Wall-pressure  measurements  were  obtained  at  these  other 
stations  to  locate  the  flow  reattachment  point  and  to  separate  the  spectral  signature 
of  the  free  shear  layer  from  that  of  the  developing  TBL.  These  wall  pressure  results 
upstream  of  the  primary  location  are  presented  in  Appendix  A. 
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2.7.2  Data  reduction  and  filtering 

In  addition  to  controlling  the  VXI  data  system,  the  HP  workstation  performed  engi¬ 
neering  unit  conversions  and  all  subsequent  analyses  of  the  digitized  data.  All  com¬ 
putations  and  graphics  were  performed  using  MATLAB™  scientific  software  running 
under  the  HP-UX  operating  system.  The  raw  binary  statistics  and  time  series  data 
were  read  into  MATLAB  and  converted  to  Pa  or  m/s  according  to  Sections  2.4.1 
and  2.5.1.  The  data  were  then  digitally  filtered  to  remove  unwanted  electrical  and 
acoustical  noise.  This  was  done  using  a  5th-order  high-pass  Butterworth  filter  with 
a  cutoff  frequency  of  100  Hz.  The  filter  was  applied  both  forwards  and  backwards  to 
eliminate  phase  lags.  To  account  for  end  effects,  500  points  were  deleted  from  both 
ends  of  the  filtered  record.  The  data  were  then  stored  in  MATLAB  binary  format 
for  further  processing. 

The  MATLAB  file  format  proved  convenient  for  two  reasons.  First,  the  binary 
format  is  efficient.  One  channel  of  filtered  data  occupied  2.5  Mbytes  of  disk  space. 
One  complete  survey  therefore  required  2.5  MBytes  x  4  channels  x  22  stations  = 
220  MBytes  of  storage,  excluding  raw  data.  Second,  the  file  header  contains  a  ma¬ 
chine  signature,  allowing  the  program  to  convert  data  written  by  a  different  platform. 
This  proved  convenient  in  transmitting  results  between  the  HP  735  and  the  486  PC. 
The  format  is  completely  documented  so  that  it  can  be  read  without  the  MATLAB 
software  using  FORTRAN  or  C  routines. 

As  detailed  in  Section  2.6.2,  the  spectral  data  were  acquired  by  the  spectrum 
analyzer  and  stored  on  an  external  GPIB  disk  drive.  The  486  PC  retrieved  the  data 
and  converted  it  from  HP  3562A  binary  format  into  MATLAB  format  using  Hewlett 
Packard  Standard  Data  Format  utility  software.  The  data  could  then  be  processed 
on  the  PC  running  MATLAB  under  Microsoft  Windows  or  transferred  to  the  HP  735 
via  Ethernet. 
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Chapter  3 

Properties  of  equilibrium  and  disturbed  flows 


3.1  Introduction 

This  chapter  presents  the  first-order  analysis  of  the  equilibrium  and  disturbed  flows. 
The  mean  features  of  each  flow  are  presented,  followed  by  single-point  statistics 
and  spectra.  Peak-event  detection  results  are  then  presented  on  wall  pressure  and 
Reynolds  stress.  For  each  perspective,  the  equilibrium  and  disturbed  cases  are  com¬ 
pared  and  contrasted  for  insight  into  the  physics  of  turbulence. 

3.2  Mean  features 

The  mean  features  of  the  two  flows  studied  are  presented  in  Table  3.1,  which  also 
summarizes  the  measurement  resolutions.  For  each  flow,  the  shear  velocity,  uT, 
was  chosen  as  the  value  which  best  fit  the  profile  to  the  log-law  relation  given  by 
Clauser  [3.1]: 

U/uT  =  2.441ni/+  +  4.9  (3.1) 

where  y+  —  yurfv.  This  assumes  that  the  standard  logarithmic  law  applies  near  the 
surface  of  the  disturbed  flow,  which  was  shown  to  be  valid  in  [3.2].  The  boundary 
layer  thickness,  6,  was  defined  as  the  location  where  U  —  0.99 Ue.  The  displacement 
thickness,  5*,  and  the  momentum  thickness,  9,  were  evaluated  using  trapezoidal  rule 
integrations  of  the  mean  velocity  profiles.  Calculated  skin  friction  coefficients,  Cf , 
were  obtained  from  the  Ludwieg-Tillman  formula: 

Cf  =  0.246R^-°'268  1  0~°'678H  (3.2) 

where  the  momentum  thickness  Reynolds  number  Reg  =  9Ue/v  and  H  =  5* / 9  is  the 
conventional  shape  factor.  “Measured”  friction  coefficients  were  obtained  from  the 
definition  of  the  shear  velocity  as: 


cf 


(3.3) 
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Table  3.1:  Nomina 


boundary-layer  properties  and  measurement  resolutions. 


Equilibrium  TBL 

Disturbed  TBL 

Single  Wire 

X-wire 

Single  wire 

X-wire 

Ue 

14.79  m/s 

16.13  m/s 

13.74  m/s 

16.18  m/s 

ur 

0.60  m/s 

.  0.64  m/s 

0.32  m/s 

0.45  m/s 

%iT  jJJe 

0.040 

0.040 

0.023 

0.028 

8 

2.72  cm 

2.95  cm 

4.05  cm 

4.30  cm 

8* 

0.43  cm 

0.46  cm 

0.97  cm 

0.89  cm 

8+ 

1064 

1202 

836 

1255 

Re& 

3045 

3358 

4599 

5397 

H 

1.39 

1.42 

1.85 

1.74 

G 

6.93 

7.41 

19.72 

15.38 

cf 

0.00329 

0.00304 

0.001431 

0.001627 

(calculated) 

Cf 

0.00324 

0.00318 

0.00108 

0.00153 

(measured) 

/,  =  v/uT 

25.6  pm 

24.6  pm 

48.4  pm 

34.2  pm 

*,  =  v/u2 

42.9  ps 

38.2  ps 

151.3  ps 

76.5  ps 

/  wire 

1.27  mm 

1.0  mm 

1.27  mm 

1.0  mm 

(50/,) 

(41/,) 

(26/,) 

(29/,) 

dmic 

0.79  mm 

0.79  mm 

0.79  mm 

0.79  mm 

(31/,) 

(32/,) 

(16/,) 

(23/,) 

A  ts 

30.5  ps 

30.5  ps 

30.5  ps 

30.5  ps 

(0.71*,) 

(0.80*,) 

(0.20*,) 

(0.40*,) 

Figure  3.1  shows  the  mean  streamwise  velocity  profiles  of  the  two  flows.  The 
equilibrium  layer  is  seen  to  exhibit  the  classic  full  profile,  while  the  disturbed  layer 
has  a  significantly  retarded  section  over  most  of  the  layer.  Alternatively,  the  dis¬ 
turbed  layer  can  be  viewed  as  having  an  accelerated  or  energized  outer  layer.  This  is 
emphasized  by  Figure  3.2,  where  the  profiles  are  normalized  by  wall  variables.  Both 
the  single  wire  and  the  X-wire  surveys  are  shown  in  this  figure,  with  the  single  wire 
surveys  approaching  much  closer  to  the  wall.  The  low  values  of  skin  friction  result 
in  highly  exaggerated  wake  regions  for  the  disturbed  flows.  The  difference  in  wake 
strengths  between  the  disturbed  flow  single  and  X-wire  runs  is  due  to  the  slightly 
different  freestream  conditions  resulting  in  differences  in  the  location  of  shear  layer 
reattachment.  This  also  accounts  for  the  dip  below  the  law  of  the  wall  seen  in  the 
disturbed  flow  X-wire  survey.  The  higher  value  of  Ue  for  this  run  places  reattachment 
closer  to  the  measurement  station  than  for  the  single-wire  run.  Therefore,  the  flow 
is  exhibiting  the  constant  mixing  length  character  of  the  free  shear  layer,  having  not 
yet  relaxed  to  a  mixing  length  proportional  to  distance  from  the  wall  as  assumed 
in  the  law  of  the  wall  [3.2,  3.3].  Figure  3.3  shows  the  mean  wall-normal  velocity 


32 


Figure  3.1:  Mean  streamwise  velocity  profiles  from  single  wire  surveys 
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Figure  3.3:  Mean  wall-normal  velocity  profiles. 

profiles.  As  expected,  the  disturbed  case  is  exhibiting  significant  inflow  over  most  of 
the  layer  until  constrained  by  the  wall.  This  inflow  manifests  itself  as  higher  values 
of  fluctuating  wall  pressure,  as  will  be  seen  in  Section  3.3.3. 

3.3  Single-point  statistics  and  spectra 

In  addition  to  the  mean  values  presented  above,  higher-order  statistics  were  com¬ 
puted  from  the  statistics  data  sets  for  both  components  of  velocity,  the  uv  product, 
and  the  wall  pressures.  The  higher-order  statistics  were  the  root-mean-square  or 
rms,  the  skewness,  and  the  kurtosis  or  flatness.  These  quantities  are  (the  square 
root  of)  the  second,  the  third,  and  the  fourth  central  moments,  respectively,  of  a 
random  variable  or  signal.  They  describe  the  probability  density  function;  i.e.,  how 
the  fluctuations  are  distributed  around  an  average  value.  Using  u  for  the  zero-mean, 
random  variable,  the  statistics  are  given  by: 


o 

X 


Equilibrium 

Disturbed 


O  o 


O 


X 

X 

X 


'Xx 


rms 


skewness 
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1  N  /  v.  \4 

kurtosis  =  (3.6) 

iV  i~ l  \Ujrms' 

where  N  is  the  number  of  samples.  The  rms  value  has  the  dimensions  of  the  variable 
of  interest  and  is  used  to  normalize  the  skewness  and  kurtosis.  It  is  a  measure  of 
the  variation  from  the  mean  value.  The  skewness  relates  to  the  symmetry  of  the 
probability  density  function;  i.e.,  it  is  an  indication  of  whether  a  variable  favors 
positive  or  negative  excursions  from  the  mean.  The  kurtosis  is  indicative  of  how 
“spiky”  the  signal  is.  A  signal  with  a  high  kurtosis  takes  on  values  far  away  from 
the  mean  relatively  more  often  than  a  signal  with  a  low  kurtosis. 

Auto  or  power  spectra  and  cross  spectra  were  also  obtained  for  u,  v,  pi,  and  p2. 
The  u  power  spectrum  is  given  by  $u(o>)  =  Tu{uj)T*{uj),  where 

Tu{u)  =  ±-J  u{t)e-^dt  (3.7) 

is  the  Fourier  transform  of  the  signal  u,  T*  is  its  complex  conjugate,  u  is  the  fre¬ 
quency,  and  t  is  the  time.  Spectra  are  decompositions  of  a  signal  into  waves  of 
different  periods  or  wavelengths.  The  value  of  the  spectrum  at  a  given  frequency  is 
the  mean  energy  in  that  wave.  As  outlined  in  Tennekes  and  Lumley  [3.4],  the  spec¬ 
trum  is  indicative  of  the  way  in  which  waves,  or  eddies,  of  different  sizes  exchange 
energy  with  each  other. 

The  spectra  as  presented  above  are  frequency  spectra  where  a  function  of  time 
is  mapped  into  the  frequency  domain.  If  the  signals  were  measured  as  a  func¬ 
tion  of  space  at  a  fixed  time,  the  mapping  would  be  from  the  spatial  domain  into 
the  wavenumber  domain.  The  transformation  between  frequency  and  wavenum¬ 
ber  spectra  can  be  made  by  assuming  Taylor’s  hypothesis,  k\  =  co/Uc,  which  gives 
$(ku  y)  =  Uc$(uj,y).  Here,  is  the  streamwise  component  of  the  wavenumber 
vector  and  Uc  is  the  convection  velocity,  taken  to  be  constant  and  equal  to  the  local 
mean  velocity,  U (y).  The  wavenumber  has  dimensions  of  length-1. 

Compared  to  mean  measurements,  higher-order  statistics  are  particularly  sensi¬ 
tive  to  probe  resolution.  Blackwelder  and  Haritonidis  [3.5],  using  hot  wires  to  detect 
the  burst  frequency  in  an  equilibrium  TBL,  found  that  only  sensors  having  a  spatial 
scale  less  than  20 lu  were  free  from  spatial  averaging  effects.  Bandyopadhyay  and 
Gad-el-Hak  [3.6]  have  recently  reviewed  the  literature  on  turbulent  boundary-layer 
flows,  concentrating  on  Reynolds  number  effects.  They  argue  that  the  velocity  probe 
resolution  expressed  in  wall  units  “Z+  should  not  be  much  larger  than  5,  where  Z  is 
the  probe  length.”  As  given  in  Table  3.1,  the  hot-wire  probe  resolutions  for  this 
study  were  in  the  range  20-50  wall  units.  As  will  be  seen  in  the  following  sections, 
this  is  small  enough  to  resolve  a  significant  portion  of  the  small-scale  structure  and 
events  that  are  the  focus  of  the  investigation,  at  least  outside  of  the  inner  layer. 
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3.3.1  Velocity  statistics 

Figures  3.4  and  3.5  show  the  rms  streamwise  and  wall-normal  velocity  fluctuations, 
respectively.  Outer-layer  scaling  is  used  since  the  disturbed  flow  perturbations  are 
large  scale,  located  outside  of  the  buffer  layer.  Again,  differences  between  the  single 
and  X-wire  runs  are  attributable  to  the  slight  differences  in  edge  velocities. 

In  Figure  3.4,  the  peak  in  the  rms  profile  of  the  equilibrium  boundary  layer  cor¬ 
responds  to  Urms/ur  =  2.42  at  y+  =  18.9,  in  agreement  with  the  criteria  of  Hussain. 
This  compares  with  a  peak  value  of  Urms/ur  ~  2.7  as  compiled  by  Bandyopadhyay 
and  Gad-el-Hak  [3.6]  from  other  investigations  where  l+  <  7.  The  probe  does  indeed 
appear  to  resolve  the  significant  small-scale  motions.  The  disturbed  flow  shows  u„nS 
and  v rms  behavior  similar  to  the  equilibrium  flow  in  the  inner  layer  and  at  the  bound¬ 
ary  layer  edge,  but  also  shows  a  higher  peak  in  turbulence  intensity  at  y/S  ~  0.2. 
This  disturbance  layer  is  the  remnants  of  the  free  shear  layer  that  formed  at  the 
step  and  contacted  the  surface  at  the  point  of  flow  reattachment.  It  moves  away 
from  the  wall  as  the  flow  progresses  downstream,  as  demonstrated  in  the  work  of 
Farabee  [3.7]. 

The  u  and  v  skewness  profiles  are  shown  in  Figures  3.6  and  3.7.  Bandyopadhyay 
and  Gad-el-Hak  [3.6,  p.  62]  point  out  that 

. . .  the  third  and  all  higher  odd-number  moments  retain  the  sign  informa¬ 
tion  and  thus  contain  valuable  statistical  information  related  to  the  co- 
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Figure  3.6:  Profiles  of  the  skewness  factor  of  the  streamwise  velocity  fluctuations 
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Figure  3.7:  Profiles  of  the  skewness  factor  of  the  wall-normal  velocity  fluctuations. 


herent  structures.  The  skewness  of  a  turbulence  quantity  can  be  thought 
of  as  representing  the  flux  of  a  stress  that  is  directly  attributable  to  coher¬ 
ent  structures.  For  example,  u 3  is  the  streamwise  flux  of  the  streamwise 
turbulence  kinetic  energy  u2,  ... 

The  u  skewness  of  the  equilibrium  flow  is  positive  near  the  wall.  This  indicates 
acceleration-dominated  velocity  fluctuations  consistent  with  the  arrival  of  high-speed 
fluid  from  regions  away  from  the  wall  (sweep  events).  In  the  log-law  region,  the 
equilibrium  flow  u  skewness  is  zero  changing  to  negative  farther  away  while  the 
v  skewness  remains  positive.  This  indicates  a  change  to  deceleration-dominated 
(ejection)  events  with  an  outward  transport  of  wall-normal  turbulence  kinetic  energy. 
For  the  disturbed  flow,  the  u  skewness  remains  positive  over  the  log-law  region 
changing  sign  at  approximately  the  same  location  as  the  equilibrium  case.  The 
disturbed  flow  v  skewness,  however,  evolves  from  positive  at  the  edge  of  the  sublayer 
to  negative  over  the  log-law  region  to  positive  again  in  the  outer  layer.  This  shows 
a  net  transport  of  wall-normal  turbulence  kinetic  energy  away  from,  towards,  and 
away  from  the  wall  in  the  three  regions.  The  skewness  behavior  for  both  flows  is 
erratic  in  the  turbulent-nonturbulent  interface  region  at  the  edge  of  the  boundary 
layer. 

Kurtosis  profiles  of  the  u  and  v  velocities  are  shown  in  Figures  3.8  and  3.9.  Unlike 
the  skewness,  it  is  always  positive,  so  no  structural  information  can  be  deduced. 
The  kurtosis  is  seen  to  be  approximately  3  in  the  overlap  region,  the  same  as  a 
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Figure  3.8:  Profiles  of  the  kurtosis  factor  of  the  streamwise  velocity  fluctuations 


Figure  3.9:  Profiles  of  the  kurtosis  factor  of  the  wall-normal  velocity  fluctuations. 
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Figure  3.10:  Reynolds  stress  profiles. 

Gaussian  signal.  The  kurtosis  increases  to  large  values  at  the  boundary  layer  edge, 
and  is  seen  to  be  increasing  as  the  wall  is  approached,  indicating  that  the  turbulence 
is  intermittent  in  both  places. 

The  Reynolds  stress  profiles,  normalized  by  the  square  of  the  edge  velocity,  are 
shown  in  Figure  3.10.  The  disturbed  flow  shows  a  peak  in  Reynolds  stress  at  the 
same  location  as  the  peaks  in  and  v„nS,  again  highlighting  the  added  energy  in 
the  outer-layer  regions.  However,  the  Reynolds  stress  near  the  wall  is  much  reduced 
from  that  of  the  equilibrium  layer,  indicating  that  the  inner  layer  is  either  still 
developing  or  is  contaminated  by  outer-layer  structures.  Figure  3.11  shows  the  same 
profiles  normalized  by  (^rms^Vms).  Similar  observations  can  be  made  from  this  figure. 

3.3.2  Velocity  power  spectra 

Power  spectra  of  the  velocity  components  will  be  considered  next.  Since  the  rms 
can  be  obtained  from  (/  the  spectral  levels  follow  the  trends  of  the  rms 

fluctuations  shown  in  Figures  3.4  and  3.5;  i.e.,  they  increase  with  increasing  y  up 
to  some  point  and  then  decrease.  For  clarity,  the  spectra  inside  and  outside  of  the 
location  of  the  peak  rms  level  will  be  shown  separately,  and  for  selected  y-locations 
only. 

The  u  power  spectra  are  shown  in  Figures  3.12  and  3.13.  From  both  figures,  it 
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Figure  3.11:  Profiles  of  turbulence  stress  correlation  coefficient. 


Figure  3.12:  Streamwise  velocity  spectra  inside  of  the  location  of  peak  equilibrium 
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Figure  3.13:  Streamwise  velocity  spectra  outside  of  the  location  of  peak  equilibrium 

'U'rms' 

is  apparent  that  inside  y/S  «  1,  the  spectra  are  broadbanded,  with  a  high-frequency 
rolloff  that  is  completed  by  10  KHz.  For  y  >  <5,  the  spectra  become  band  limited  with 
a  peak  in  the  neighborhood  of  200  Hz.  The  band-limited  behavior  is  the  result  of 
irrotational  motion  in  the  freestream  induced  by  the  turbulent  boundary  layer.  This 
statement  is  supported  by  the  u-v  cross  spectral  data  of  the  present  investigation  as 
well  as  that  of  Wilczynski  [3.8]. 

Careful  examination  of  Figure  3.12  reveals  that  the  spectral  levels,  though  very 
similar,  are  increasing  slightly  through  y / 5  —  0.022.  Note  that  the  three  lower  sets 
of  curves  have  been  shifted  down  by  10,  20,  and  30  dB,  respectively,  for  clarity.  The 
disturbed  flow  low-  and  high-frequency  components  are  increased  relative  to  the 
equilibrium  flow.  Also,  a  linear  midfrequency  region  is  perceptible  in  the  disturbed 
flow  spectra  for  y/8  >  0.005.  The  extent  of  the  linear  range  grows  with  increasing 
y.  A  line  of  -5/3  slope  is  shown  for  comparison.  Such  scaling  is  also  found  in  other 
wake  flows,  such  as  the  axisymmetric  jet  studied  by  Champagne  [3.9].  Equilibrium 
boundary  layers  also  exhibit  this  behavior  in  the  frequency  range  bridging  the  high 
and  low  wavenumber  spectra.  This  is  known  as  the  “inertial  subrange” ,  the  extent  of 
which  is  dependent  on  the  turbulence  Reynolds  number  [3.4].  It  was  not  apparent  in 
the  equilibrium  data  of  this  investigation  because  the  turbulence  Reynolds  number 
is  not  high  enough. 

Continuing  across  the  boundary  layer  in  Figure  3.13,  the  equilibrium  layer  spec- 
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Figure  3.14:  Selected  wall-normal  velocity  spectra  inside  of  the  location  of  peak 
disturbed  u„„s. 

tral  levels  are  seen  to  peak  at  y/5  =  0.022  to  0.025,  consistent  with  Figure  3.4.  As 
in  Figure  3.12,  the  lower  curve  sets  have  been  shifted  down  by  10,  15,  and  20  dB, 
respectively.  From  y/5  —  0.025  to  0.17,  the  equilibrium  levels  decrease  slowly.  The 
range  of  /~ 5/3  scaling  in  the  disturbed  flow  continues  to  increase,  covering  a  maxi¬ 
mum  frequency  range  of  200<  /  <2000  Hz  at  y/5  =  0.17.  This  coincides  with  the 
peak  in  the  disturbed  flow  turbulence  intensities,  confirming  that  the  disturbance 
layer  is  the  source  of  the  f~5^3  scaling. 

The  v  spectra  of  Figures  3.14  and  3.15  show  behavior  similar  to  the  u  spectra. 
At  the  closest  measurement  locations,  the  equilibrium  and  disturbed  spectra  are 
very  similar,  with  a  wider  low-frequency  flat  section  and  steeper  rolloff  than  the 
corresponding  u  spectra.  Moving  away  from  the  wall,  the  disturbed  v  spectral  levels 
increase  relative  to  the  equilibrium  case  because  of  the  disturbance  layer.  Again, 
the  peak  disturbed  flow  spectral  levels  correspond  in  location  to  the  peak  in  v^. 
However,  /-5//3  scaling  is  not  apparent. 

The  band-limited  character  is  also  more  prevalent,  occurring  closer  to  the  wall 
than  for  the  u  spectra.  The  disturbed  flow  shows  a  peak  at  /  «  150  Hz  at  y/5  =  0.35. 
As  y  increases  further,  the  peak  becomes  more  pronounced  and  moves  to  lower 
frequencies.  This  behavior  was  observed  in  the  step  flow  experiments  of  Farabee  [3.7] 
at  multiple  downstream  locations.  He  suggested  that  the  shift  to  lower  frequencies 
in  the  disturbed  flow  spectra  can  be  explained  by  considering  structures  that  have 
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Figure  3.15:  Selected  wall-normal  velocity  spectra  outside  of  the  location  of  peak 
disturbed  vrms. 


length  scales  equal  to  the  distance  from  the  wall.  Assuming  further  that  these 
structures  convect  with  a  speed  equal  to  the  edge  velocity  results  in 

,  ,  ,  1  Ue 

Jpeak{y )  —  0  (3.8) 

Z7T  y  v  ' 

Frequencies  predicted  by  this  equation  are  indicated  in  Figure  3.15. 

The  equilibrium  flow,  on  the  other  hand,  does  not  show  a  peak  until  y/S  =  0.66, 
where  it  occurs  at  approximately  200  Hz.  This  peak  also  moves  to  lower  frequencies 
with  increasing  y.  By  y/S  «  1.5  (not  shown),  the  peak  is  at  approximately  100  Hz. 
The  movement  of  the  equilibrium  flow  peak  is  linked  to  the  induced  irrotational 
motion  at  the  boundary  layer  edge. 

3.3.3  Wall-pressure  statistics 

Fluctuating  wall  pressure  statistics  are  functions  of  both  flow  Reynolds  number 
and  sensor  resolution.  The  Reynolds  number  determines  the  scales  of  turbulence 
structures  present  in  the  flow.  Sensor  resolution,  generally  expressed  in  wall  units 
as  d+  =  d/lu  =  duT/v  where  d  is  the  effective  diameter,  determines  the  scale  of 
the  turbulence  structures  that  will  be  spatially  averaged  out.  The  wall-pressure 
statistics  were  examined  for  consistency  with  previous  research  and  for  any  effects 
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Table  3.2:  Wall  pressure  statistics,  averaged  over  the  X-wire  surveys. 


Transducer 

rms 

(Pa) 

Skewness 

Kurtosis 

Equilibrium  TBL,  dm{C 

=  32  lv 

Pi 

1.535 

0.0496 

4.2627 

P2 

1.465 

0.0376 

4.4550 

Disturbed  TBL,  = 

:  23/^ 

Pi 

3.903 

0.1426 

4.1281 

P2 

4.123 

0.1861 

4.2201 

Uncertainties 

0.011 

0.039 

0.342 

of  external  disturbances.  The  results  are  presented  in  Table  3.2.  The  values  listed 
are  the  averages  over  the  X-wire  surveys. 

The  rms  levels  for  the  equilibrium  flow  equate  to  3.15th,  and  3.01tw  for  trans¬ 
ducers  pi  and  P2,  respectively.  This  is  in  excellent  agreement  with  other  experi¬ 
mental  results  compiled  by  Farabee  and  Casarella  [3.10]  for  transducer  sizes  less 
than  d+  =  80.  The  disturbed-flow  rms  levels  are  approximately  three  times  those  of 
the  equilibrium  flow,  consistent  with  the  mean  inflow  shown  in  Figure  3.3,  and  the 
increased  turbulence  activity  in  the  vicinity  of  reattachment  [3.11]. 

External  disturbances  that  could  possibly  affect  the  pressure  measurements  in¬ 
clude  facility  noise  and  the  presence  of  the  hot-wire/traverse  assembly.  As  noted 
in  Sections  2.2  and  2.7.2,  the  facility  noise  is  limited  to  frequencies  of  50  Hz  and 
below,  and  is  removed  by  digital  filtering.  Farabee  and  Casarella  [3.10]  note  that 
low-frequency  components  contribute  little  to  the  turbulent  rms  level  and  can  be 
neglected.  The  effects  of  the  hot-wire/traverse  assembly  for  the  equilibrium  flow 
can  be  seen  in  Figure  3.16.  The  hot  wire  influences  transducer  pi  by  increasing 
the  pressure  levels  by  approximately  3  percent  when  near  the  wall.  Transducer  p2, 
situated  upstream,  is  unaffected.  This  effect  is  also  evident  in  the  power  spectra, 
discussed  in  Section  3.3.4.  The  disturbed  flow  case  is  shown  in  Figure  3.17.  For 
this  case,  the  wire/traverse  assembly  affects  both  transducers,  increasing  the  rms 
levels  by  3-5  percent  over  the  average  values.  The  lower  mean  velocity  levels  in  the 
retarded  layer  allow  the  traverse  influences  to  penetrate  further  upstream  than  for 
the  equilibrium  layer.  The  skewness  and  kurtosis  coefficients  for  both  flows  were 
examined  in  a  like  manner  and  found,  within  the  uncertainty  of  the  measurements, 
to  be  unaffected  by  the  traverse. 

The  skewness  and  kurtosis  for  random  data  with  a  Gaussian  probability  distri¬ 
bution  are  zero  and  three,  respectively.  Referring  to  Table  3.2,  the  equilibrium  flow 
wall-pressure  skewness  is  seen  to  be  essentially  zero,  while  the  disturbed  flow  value 
is  slightly  positive.  Both  flows  exhibit  high  values  of  the  kurtosis,  indicating  that 
the  signals  are  characterized  by  spikes  of  large  amplitude.  The  disturbed  flow  ex¬ 
periences  relatively  more  positive,  large  amplitude  fluctuations,  consistent  with  the 
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Figure  3.16:  Effect  of  the  hot  wire  traverse  assembly  on  the  equilibrium  flow  wall 
pressure. 


negative  mean  V  velocity.  The  equilibrium  flow,  on  the  other  hand,  is  symmetric, 
with  nearly  equal  numbers  of  positive  and  negative  large-amplitude  excursions. 

Further  discussion  of  the  wall-pressure  skewness  is  necessary.  Schewe  [3.12]  stud¬ 
ied  the  dependence  of  measured  wall-pressure  statistics  on  transducer  size.  He  con¬ 
cluded  that  as  transducer  size  increases,  spatial  averaging  drives  the  statistics  to 
their  Gaussian  values.  His  results  showed  that  skewness  ->•  0  for  d+  >  «  50  and 
kurtosis  ->  3  for  d+  >«  325  at  Ree  =  1400.  Interpolating  Schewe’s  results,  the  skew¬ 
ness  for  the  present  investigation  equilibrium  flow  should  be  approximately  -0.15 
instead  of  zero. 

This  discrepancy  was  resolved  by  a  series  of  tests  that  focused  on  the  processing 
of  the  pressure  signal,  namely,  the  method  of  high-pass  filtering  to  remove  unwanted 
noise.  The  hot  wire  traverse  assembly  was  removed  from  the  tunnel  to  eliminate  as 
many  external-flow  disturbances  as  possible,  and  wall-pressure  data  were  recorded 
for  a  nominal  16  m/s  equilibrium  TBL.  Three  methods  of  high-pass  filtering  were 
used  to  remove  the  remaining  facility  noise: 

1.  Analog  filter  prior  to  digitization,  elliptic  type  of  order  5,  cutoff  frequency=100  Hz. 

2.  Digital  filter  after  digitization,  Butterworth  type  of  order  5,  single  pass,  cutoff 
frequency=100  Hz. 
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Figure  3.17:  Effect  of  the  hot  wire  traverse  assembly  on  the  disturbed  flow  wall 
pressure. 

3.  Digital  filter  after  digitization,  Butterworth  type  of  order  5,  double  pass,  cutoff 
frequency=100  Hz. 

The  results  are  presented  in  Table  3.3.  The  analog  and  single-pass  digital  filters 
introduce  phase  lags  into  the  data  and  give  equivalent  results  for  the  skewness,  which 
are  consistent  with  the  results  of  Schewe.  The  double-pass  digital  filter,  on  the  other 
hand,  introduces  no  phase  lags  and  produces  results  equivalent  to  Table  3.2.  The 
ways  in  which  phase  lags  influence  the  skewness  are  not  understood,  but  the  effects 
are  clearly  evident.  One  must  therefore  be  careful  when  reviewing  results  in  the 
literature,  as  many  researchers  neglect  to  specify  the  exact  type  of  filtering  performed 


Table  3.3:  Effects  of  high-pass  filtering  on  equilibrium  TBL  wall  pressure  skewness. 


Filter  type 

Skewness 

Analog,  elliptic 

-0.1996 

Digital,  Butterworth, 

-0.1999 

single  pass 

Digital,  Butterworth, 

0.0212 

double  pass 
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Frequency  (Hz) 


Figure  3.18:  Equilibrium  and  disturbed  TBL  wall  pressure  spectra  from  microphone 
Pi-  The  noise-canceled  equilibrium  spectrum  is  from  Farabee  and  Casarella  [3.10]. 


on  the  data.  Within  the  precision  of  the  measurements,  the  rms  and  kurtosis  values 
showed  no  dependence  on  filter  choice.  Peak  event  analysis,  as  detailed  in  Section  3.4, 
is  likewise  unaffected. 

3.3.4  Wall-pressure  power  spectra 

Figure  3.18  presents  the  wall-pressure  power  spectra  for  both  the  equilibrium  and 
disturbed  flows.  For  each  case,  measurements  are  shown  with  the  X-wire  probe 
closest  to  and  furthest  away  from  the  wall.  This  illustrates  again  the  effect  of  the 
traverse  and  hot  wire  probe  on  the  pressure  measurements:  below  1000  Hz,  the 
spectral  levels  with  the  probe  close  to  the  wall  are  slightly  above  those  with  the 
probe  far  away.  Integration  of  these  spectra  over  the  frequency  range  corresponding 
to  the  digitized  data  (100-12800  Hz)  results  in  the  same  percentage  increase  reported 
in  Section  3.3.3. 

The  equilibrium  spectra  show  the  expected  broad-banded  behavior.  Also  in¬ 
cluded  in  Figure  3.18  is  an  equilibrium-flow  spectrum.  This  spectrum  was  obtained  in 
the  same  facility  at  a  flow  speed  Ue  =  15.2  m/s,  using  a  noise  cancellation  technique 
to  eliminate  facility  noise  [3.10].  The  technique  involves  subtracting  the  outputs  of 
two  microphones  located  at  the  same  streamwise  location  but  separated  by  a  span- 
wise  distance.  The  spectrum  is  shown  to  illustrate  the  regions  where  the  current  data 
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are  contaminated  by  unwanted  acoustic  sources.  The  first  region  is  the  frequency 
range  below  40  Hz.  The  noise  source  here  is  the  tunnel’s  centrifugal  blower,  which 
was  identified  by  Farabee  [3.7]  and  Helal,  et  al.  [3.13].  The  second  region  of  con¬ 
tamination  is  in  the  range  150-300  Hz.  This  source,  identified  by  Wilczynski  [3.8], 
is  Strouhal  shedding  from  the  hot-wire  traverse.  The  fact  that  these  do  not  show 
up  in  the  velocity  spectra  confirms  that  they  are  acoustical.  The  blower  noise  was 
removed  by  digital  high-pass  filtering  at  100  Hz  during  data  reduction.  At  this  stage 
of  the  signal  processing,  no  attempt  was  made  to  remove  the  traverse-related  noise, 
and  it’s  presence  will  be  acknowledged.  However,  calculations  made  on  wall-pressure 
data  collected  in  a  clean  wind  tunnel  revealed  no  discernable  effects  of  the  traverse 
on  peak-event  detection  results. 

Compared  to  the  equilibrium  pressure  spectra,  the  disturbed-flow  pressure  spec¬ 
tra  exhibit  excess  low-frequency  components,  as  was  also  seen  in  the  disturbed  flow 
velocity  spectra.  Farabee  [3.7]  showed  that  this  region  of  the  spectrum  scales  on 
“disturbance  variables” ,  given  as 

$p{u,x/h)U{5d)  uh 

(puv)26d  VS'  Ue  (3'9) 

Here,  u>  is  the  frequency  in  rad/s,  the  length  scale,  5d,  is  the  distance  that  the 
disturbance  layer  is  away  from  the  wall,  h  is  the  step  height,  and  the  pressure  scale, 
—puv,  is  the  maximum  value  of  the  turbulence  Reynolds  stress  measured  in  the 
disturbance  layer.  This  scaling  means  that  the  disturbance  layer  is  the  source  of 
the  excess  low-frequency  energy.  The  region  of  /~5/3  behavior,  as  in  the  u  velocity 
spectra,  is  characteristic  of  free  shear  layer  flows. 

The  spectra  for  the  two  flows  show  identical  behavior  in  the  high-frequency  range 
above  3000  Hz.  The  source  of  these  high-frequency  fluctuations  is  the  small-scale 
structures  near  the  wall  [3.10,  3.14],  Since  the  two  flows  show  nearly  identical  near¬ 
wall  turbulence  levels  (Figure  3.4),  their  similarity  in  pressure  spectral  behavior  is 
not  surprising.  However,  a  word  of  warning  is  in  order.  Farabee  had  shown  for 
equilibrium  flows  that  the  high-frequency  data  collapsed  on  inner  variables.  Since 
the  values  of  tu  between  the  equilibrium  and  disturbed  flows  differ  by  a  factor  of 
two,  the  disturbed  flow  will  not  collapse  on  inner  variables.  This  raises  the  question 
of  the  proper  inner  variable  scaling  for  nonequilibrium  flows.  This  is  also  related  to 
the  questions  of  inner/outer  layer  coupling,  since  uT  depends  on  the  outer  layer,  but 
the  near  wall  turbulence  structure  of  urms  and  vrms  appears  to  be  independent. 


3.4  Peak-event  detection 

The  analysis  up  to  this  point  has  been  in  the  form  of  time  averages  and  Fourier  spec¬ 
tra.  Each  of  these  involved  integration  of  the  signals  over  long  times.  Attention  will 
now  be  focused  on  those  sections  of  the  time  series  that  are  believed  to  be  associated 
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Figure  3.19:  Sample  time  record  of  the  equilibrium  flow  fluctuating  wall  pressure. 


with  the  active,  or  turbulence-producing  motions.  This  requires  a  method  of  identi¬ 
fying  or  detecting  turbulence-related  portions  of  a  signal.  Figure  3.19  shows  approxi¬ 
mately  1000  samples  of  the  wall-pressure  signal  from  the  equilibrium  boundary-layer 
flow.  High-amplitude  fluctuations  are  apparent  and  are  intuitively  associated  with 
the  turbulence.  Hence,  the  most  basic  scheme  is  to  declare  an  event  when  the  signal 
exceeds  a  given  threshold  value.  The  threshold  is  expressed  as  a  multiple  of  the 
rms  value.  This  is  known  as  peak-event  detection.  For  example,  peak  wall-pressure 
events  occur  when  \p\  >  Kp^,  where  k  is  the  threshold  level.  In  the  discussion  to 
follow,  positive  and  negative  events  will  be  distinguished  by  the  sign  of  k: 

k>0  =»  positive  event,  p  >  a zprms 

k  <  0  =>  negative  event,  p  <  l  •  J 

The  specific  peak-event  detection  algorithm  was  that  used  by  Karangelen  [3.15] 
and  Wilczynski  [3.8].  The  time  record  was  scanned,  and  four  points  for  each  event 
identified:  the  point  where  the  signal  exceeded  the  given  threshold  (the  trigger 
point),  the  point  of  peak  amplitude,  and  the  two  zero  crossings  bracketing  the  trigger 
point.  Figure  3.20  illustrates  this  for  a  threshold  level  k  =  2,  using  the  two  positive 
events  from  the  center  of  Figure  3.19.  The  indices  to  these  points  were  stored, 
from  which  the  following  event  statistics  were  compiled:  the  number  of  events,  the 
average  event  shape,  the  average  event  duration,  the  distribution  of  the  time  between 
consecutive  events,  and  the  contribution  of  the  events  to  the  total  record  length  and 
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Figure  3.20:  Consecutive  positive  wall-pressure  peak  events  for  threshold  k  =  2. 
to  the  signal’s  fluctuating  power. 

Average  event  shapes  were  computed  by  setting  a  window  around  each  event, 
centered  about  the  peak  magnitude,  and  ensemble  averaging.  The  event  duration, 
At,  was  defined  as  the  time  between  the  zero  crossings.  The  zero  crossings  were 
also  used  to  determine  the  fractional  contribution  of  the  events  to  the  total  record 
length  and  to  the  fluctuating  power.  The  power  fraction  was  computed  relative  to 
the  signal’s  variance  as  opposed  to  the  signal’s  rms  value,  as  reported  by  previous 
investigators.  For  comparison,  the  two  quantities  are  related  by 

rms  fraction  =  (power  fraction)1^2  (3.11) 

The  power  fraction  is  felt  to  be  more  meaningful  and  less  ambiguous,  since  rms 
fractions  often  total  to  values  greater  than  one.  The  time  between  consecutive  events, 
AT,  was  defined  as  the  time  between  trigger  points.  The  analysis  was  performed  on 
both  the  wall  pressure  and  the  Reynolds  stress  time  records. 

3.4.1  Wall-pressure  peak  events 

Table  3.4  summarizes  the  equilibrium  and  disturbed  flow  wall-pressure  peak-event 
detection  statistics.  Results  are  given  separately  for  positive  and  negative  events, 
for  threshold  values  k  =  1-3.  The  equilibrium  flow  results  are  in  close  agreement 
with  the  results  of  Karangelen  [3.15]  and  Wilczynski  [3.8].  The  large  amplitude 
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Table  3.4:  Wall  pressure  peak-event  detection  statistics  for  transducer 


K 

Time 

fraction 

Power 

fraction 

A  t+ 

A  T+ 

Number  of 

events 

Frequency 
per  1000f+ 

Equilibrium 

+ 1  0.330 

0.487 

17.1 

51.9 

5034 

19.3 

-1 

0.343 

0.459 

17.4 

50.7 

5150 

19.7 

+2 

0.124 

0.308 

19.6 

157.7 

1654 

6.3 

-2 

0.114 

0.260 

19.3 

169.6 

1540 

5.9 

+3 

0.029 

0.118 

19.4 

666.8 

391 

1.5 

-3 

0.022 

0.089 

17.1 

765.0 

342 

1.3 

Disturbed 

+1  0.343 

0.494 

18.0 

52.6 

2472 

19.0 

-1 

0.352 

0.454 

19.1 

54.1 

2405 

18.5 

+2 

0.137 

0.326 

19.8 

144.0 

901 

6.9 

-2 

0.105 

0.246 

20.7 

196.6 

659 

5.1 

+3 

0.039 

0.147 

19.7 

503.9 

257 

2.0 

-3 

0.019 

0.075 

18.2 

965.6 

134 

1.0 

events  are  seen  to  contribute  a  majority  to  the  total  power.  Events  exceeding  3prms 
occur  5  percent  of  the  time  but  contribute  almost  21  percent  to  the  total  power. 
To  appreciate  the  significance,  Karangelen  compared  equilibrium  turbulent  wall- 
pressure  data  to  a  random,  Gaussian,  signal.  To  facilitate  the  comparison,  she  used 
only  those  portions  of  the  signals  that  actually  exceeded  the  threshold.  She  showed 
that  for  a  threshold  k  =  3,  the  turbulent  wall  pressure  events  occupied  1  percent  of 
the  time  but  contributed  14  percent  to  the  total  power,  while  the  Gaussian  signal 
events  occupied  0.4  percent  of  the  time  while  contributing  only  4  percent  to  the 
power.  This  is  consistent  with  the  high  kurtosis  values  seen  in  Table  3.2. 

In  Table  3.4,  the  average  event  durations  and  time  between  events  are  given  in 
viscous  time  units,  t+  =  t/tv.  Also,  the  frequency  of  realizations  is  reported  in  the 
form  of  number  of  events  per  1  Cl00t~~ .  The  choice  of  viscous  or  wall  variable  scaling 
reflects  the  belief  that  wall-pressure  events  are  related  to  the  turbulent  burst  fre¬ 
quency  [3.8,  3.15-3.17].  The  average  event  duration  is  seen  to  be  relatively  constant 
with  threshold  value,  ranging  from  17  to  20  viscous  time  units.  As  expected,  the 
average  time  between  events  increases  with  increasing  threshold. 

The  disturbed-flow  event  statistics  are  remarkably  similar  to  the  equilibrium-flow 
results.  The  sum  of  the  positive  and  negative  event  contributions  to  the  total  time 
and  power  for  a  given  threshold  are  equivalent  to  the  equilibrium  flow  case.  However, 
the  disturbed  flow  results  show  relatively  more  positive  than  negative  events,  as 
seen  in  the  power  fraction  and  event  frequency  values.  This  is  consistent  with  the 
positive  skewness  results  shown  in  Table  3.2.  It  is  noteworthy  that  event  durations, 
times  between  events,  and  event  frequencies  for  the  disturbed  flow  reduce  to  the 


52 


NSWCDD/TR-95/222 


Figure  3.21:  Average  positive  wall-pressure  peak-event  shapes,  k  —  2. 

equilibrium  values  when  scaled  on  wall  variables,  particularly  in  light  of  the  drastic 
differences  in  the  number  of  events  observed.  This  implies  that  the  wall-pressure 
peak  events  are  indeed  driven  by  the  near-wall  flow,  and  that  effects  of  the  outer 
layer  disturbances  on  the  near-wall  motion  are  secondary. 

Average  event  shapes  for  a  threshold  level  of  2  are  presented  in  Figures  3.21 
and  3.22.  Shapes  for  the  equilibrium  flow  are  very  similar  to  those  reported  by 
Karangelen  [3.15]  and  Wilczynski  [3.8].  The  positive  events  show  an  overshoot  on 
the  return  from  the  peak  value.  The  negative  events  are  inverted  mirror  images  of 
the  positive  events.  The  disturbed  flow  average  event  shapes  match  the  equilibrium 
case  over  the  central  portions  of  the  events,  but  have  broader  bases.  The  portions 
of  the  pressure  signal  away  from  the  event  also  appear  to  be  more  highly  correlated 
than  in  the  equilibrium  flow,  where  the  ensemble  average  oscillates  about  zero.  These 
observations  seem  to  be  consistent  with  the  increased  lower-frequency  components 
produced  by  the  step  flow.  The  disturbed-flow  event  shapes  are  not  symmetric,  with 
the  rise  in  pressure  to  the  positive  peaks  being  more  gradual  than  the  falloff  beyond 
the  peaks.  Like  the  equilibrium  flow,  the  disturbed-flow  negative  pressure  events  are 
inverted  mirror  images  of  their  positive  counterparts. 

Figures  3.23  and  3.24  show  the  distribution  of  the  time  between  positive  pressure 
events  for  a  threshold  level  of  1.  The  time  axis  has  been  normalized  by  the  mean 
time  between  events  for  that  threshold.  The  equilibrium  flow  shows  the  same  log¬ 
normal-like  distribution  observed  by  Karangelen  and  Wilczynski.  Karangelen  [3.15] 
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Figure  3.22:  Average  negative  wall-pressure  peak-event  shapes,  k  =  2. 

showed  that  a  Gaussian  signal  has  a  log-normal  distribution  for  the  time  between 
events.  As  we  have  seen,  the  equilibrium  wall  pressure  signal  has  a  probability 
distribution  that  is  similar  to  Gaussian,  but  with  an  elevated  value  for  the  kurtosis. 
In  contrast,  the  disturbed-flow  distribution  is  much  more  uniform.  Results  for  the 
negative  events  are  very  similar  and  are  not  shown. 

3.4.2  Reynolds  stress  peak  events 

Peak-event  detection  statistics  for  the  Reynolds  stress  are  presented  in  Table  3.5. 
Results  are  given  for  a  threshold  level  k  =  2  at  four  locations  across  the  boundary 
layer.  Here,  the  threshold  is  relative  to  the  product  UrmsVrms.  Unlike  the  wall  pressure 
events,  wall  variable  scaling  does  not  collapse  the  equilibrium  and  disturbed  results. 
This  is  perhaps  not  surprising  since  the  disturbed  flow  is  dominated  by  large  scale 
eddy  structures  that  span  the  entire  boundary  layer.  Even  more  so  than  for  the  wall 
pressure  events,  the  large  amplitude  uv  events  contribute  a  majority  to  the  signal 
power,  totaling  over  70  percent  for  both  flows.  Turbulence  production  is  clear  for 
both  flows;  the  power  fractions  for  the  negative  events  are  roughly  3-5  times  the 
levels  of  the  positive  events.  The  exception  is  the  disturbed  TBL  near  the  wall, 
where  positive  events  contribute  13  percent  more  than  the  negative  events.  At  this 
location,  the  total  disturbed-flow  Reynolds  stress  is  very  low,  as  seen  in  Figures  3.10 
and  3.11.  This  indicates  again  that  the  disturbed-flow  inner  layer  has  either  not 
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Figure  3.23:  Distribution  of  the  time  between  equilibrium-flow  positive  wall-pressure 
peak  events,  k  —  1. 

reached  a  fully  developed  state  or  is  contaminated  by  the  outer-layer  structures. 

Considering  event  duration,  spacing,  and  frequency,  the  equilibrium  flow  results 
are  roughly  constant  across  the  boundary  layer.  This  corresponds  to  the  relatively 
constant  levels  of  Reynolds  stress  shown  in  Figure  3.10.  The  disturbed-flow  event 
durations  are  even  more  constant  across  the  layer.  Again,  this  is  consistent  with 
a  characteristic  structure  that  spans  the  boundary  layer.  The  disturbed  flow  time 
between  events,  and  hence  the  event  frequency,  however,  are  not  constant.  The 
behavior  of  these  statistics  follows  the  overall  level  of  Reynolds  stress,  which  was 
shown  to  peak  at  y/5  =  0.2.  As  the  Reynolds  stress  increases,  large-amplitude 
negative  events  occur  more  frequently. 

Average  uv  event  shapes  across  the  boundary  layers,  for  a  threshold  level  of  —2, 
are  presented  in  Figures  3.25-3.28.  The  differences  in  event  durations  noted  in 
Table  3.5  are  apparent.  Figures  3.26  and  3.27  are  for  y/5  —  0.08  and  0.25,  cor¬ 
responding  to  the  log-law  region.  Here,  the  event  shapes  and  durations  are  the 
most  similar,  but  wall  variable  scaling  still  cannot  be  said  to  collapse  the  data.  Fig¬ 
ure  3.28  is  for  y/5  =  0.65,  which  is  outside  the  location  of  the  disturbance  layer.  The 
magnitude  of  the  events  has  increased;  a  consequence  of  the  increased  intermittency. 

The  Reynolds  stress  events  were  investigated  further  by  discriminating  the  in¬ 
dividual  signs  of  the  u  and  v  signals;  i.e.  the  quadrants  of  the  u-v  plane.  This  is 
known  as  quadrant-event  detection.  Events  detected  according  to  quadrant  criteria 
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Figure  3.24:  Distribution  of  the  time  between  disturbed-flow  positive  wall-pressure 
peak  events,  k  =  1. 

are  designated  Q  events.  An  event  is  declared  present  as  before,  when  the  uv  prod¬ 
uct  exceeds  a  specified  threshold.  However,  the  limits  of  the  event  are  determined 
by  when  the  signal  “enters”  and  “leaves”  the  quadrant.  It  should  be  noted  that 
because  of  the  discrete  nature  of  the  u  and  v  time  sequences,  a  quadrant  change 
can  occur  with  or  without  a  zero  crossing  of  the  uv  product.  This  means  there  is 
not  necessarily  a  one-to-one  correspondence  between,  for  example,  second  or  fourth 
quadrant  events  and  negative  uv  events. 

The  event  statistics  for  the  turbulence-producing  quadrant  events  across  the 
boundary  layers  are  presented  in  Table  3.6.  Positive  Reynolds  stress  or  negative 
uv  corresponds  to  the  second  (Q2)  and  fourth  (Q4)  quadrants  of  the  u-v  plane.  For 
the  equilibrium  flow,  Q2  events  occur  more  often  and  contribute  more  power  than 
the  Q4  events  at  all  tabulated  locations.  In  contrast,  the  disturbed  flow  exhibits 
a  change  from  Q4  to  Q2  dominance  with  increasing  y.  To  examine  this  further, 
the  event  frequencies  for  all  quadrants  across  the  boundary  layers  are  plotted  in 
Figures  3.29  and  3.30. 

In  Figure  3.29,  Q2  and  Q4  quadrant  events  are  seen  to  dominate  the  equilibrium 
flow  over  most  of  the  boundary  layer.  For  y/6  >  0.7,  Ql,  Q3,  and  Q4  events  occur 
with  roughly  the  same  frequency,  while  Q2  events  remain  dominant.  In  Figure  3.30, 
the  disturbed  flow  is  seen  to  be  dominated  by  Ql  events  at  the  closest  measurement 
station.  This  is  further  evidence  that  the  inner  layer  is  contaminated.  Moving  away 
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Table  3.5:  Reynolds  stress  peak-event  detection  statistics. 


y+  (yS) 

Time 

fraction 

Power 

fraction 

A  t+ 

AT+ 

Number  of 

events 

Frequency 
per  1000t+ 

Equilibrium,  k 
25.9  (0.022) 

=  +2 
0.053 

0.185 

14.0 

261.7 

998 

3.8 

93.1  (0.077) 

0.039 

0.138 

7.5 

191.2 

1364 

5.2 

300.0  (0.250) 

0.041 

0.152 

8.4 

202.0 

1293 

4.9 

781.1  (0.650) 

0.044 

0.146 

9.1 

203.7 

1280 

4.9 

Equilibrium,  k 
25.9  (0.022) 

=  -2 

0.148 

0.576 

14.7 

99.5 

2625 

10.0 

93.1  (0.077) 

0.183 

0.628 

11.7 

63.8 

4095 

15.7 

300.0  (0.250) 

0.191 

0.600 

13.8 

72.3 

3613 

13.8 

781.1  (0.650) 

0.182 

0.677 

15.6 

85.4 

3054 

11.7 

Disturbed,  /t  = 
24.9  (0.020) 

+2 

0.091 

0.488 

5.6 

61.7 

2110 

16.2 

97.2  (0.077) 

0.069 

0.275 

5.7 

82.4 

1581 

12.1 

320.2  (0.255) 

0.040 

0.151 

4.5 

113.1 

1152 

8.8 

833.6  (0.664) 

0.063 

0.182 

5.6 

88.5 

1471 

11.3 

Disturbed,  k  — 
24.9  (0.020) 

-2 

0.103 

0.352 

6.7 

64.9 

2007 

15.4 

97.2  (0.077) 

0.163 

0.523 

8.2 

50.5 

2578 

19.8 

320.2  (0.255) 

0.219 

0.653 

9.2 

41.8 

3112 

23.9 

833.6  (0.664) 

0.156 

0.697 

7.9 

50.6 

2572 

19.7 

from  the  wall,  first  Q2,  then  Q4  events  become  prevalent,  with  Q4  dominating  over 
most  of  the  layer.  The  changeover  point  from  Q2  to  Q4  occurs  at  the  disturbance 
location,  y/5  «  0.2. 

Another  view  of  the  situation  near  the  wall  is  provided  in  Figures  3.31  and  3.32. 
Here,  the  peak  value  of  each  quadrant  event  exceeding  a  threshold  value  k  =  1 
is  plotted  on  the  u-v  plane,  providing  a  measure  of  the  joint  probability  density 
function.  The  dominance  of  the  equilibrium  flow  Q2  events  is  again  evident.  In 
Figure  3.32,  the  dominance  of  the  Ql  and  Q4  events,  both  in  magnitude  as  well  as  in 
number,  is  also  confirmed.  These  two  quadrants  involve  positive  u,  or  accelerating 
perturbations.  This  is  consistent  with  the  larger  values  of  disturbed  flow  u-skewness 
shown  in  Figure  3.6. 


3.5  Discussion 

Overall  flowfield  and  wall-pressure  results  have  been  presented  on  two  turbulent 
boundary-layer  flows:  an  equilibrium  layer  and  a  disturbed  layer  formed  downstream 
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Figure  3.25:  Average  uv  peak-event  shapes  at  y/S  —  0.02,  k  =  —2. 


Figure  3.26:  Average  uv  peak-event  shapes  at  y/S  =  0.08,  k  =  —2. 
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Figure  3.27:  Average  uv  peak-event  shapes  at  y/5  =  0.25,  k  —  —  2 


Figure  3.28:  Average  uv  peak-event  shapes  at  y/5  =  0.65,  k  =  —2 
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Table  3.6:  Reynolds  stress  quadrant-event  detection  statistics,  k  =  -2. 


y+  (yS) 

Time 

fraction 

Power 

fraction 

A  T+ 

Number  of 
events 

Frequency 
per  1000f+ 

Equilibrium,  quadrant  2 

25.9  (0.022) 

0.104 

0.495 

13.94 

133.35 

1957 

7.5 

93.1  (0.077) 

0.101 

0.425 

10.92 

107.61 

2427 

9.3 

300.0  (0.250) 

0.100 

0.367 

12.75 

127.79 

2043 

7.8 

781.1  (0.650) 

0.109 

0.562 

13.38 

122.96 

2122 

8.1 

Equilibrium,  quadrant  4 

25.9  (0.022) 

0.040 

0.081 

15.24 

382.61 

683 

2.6 

93.1  (0.077) 

0.076 

0.201 

11.56 

152.43 

1711 

6.6 

300.0  (0.250) 

0.087 

0.233 

14.36 

164.18 

1589 

6.1 

781.1  (0.650) 

0.069 

0.115 

18.89 

273.78 

951 

3.6 

Disturbed,  quadrant  2 

24.9  (0.020) 

0.048 

0.132 

8.56 

179.23 

727 

5.6 

97.2  (0.077) 

0.072 

0.176 

10.02 

138.18 

938 

7.2 

320.2  (0.255) 

0.123 

0.455 

8.74 

70.81 

1835 

14.1 

833.6  (0.664) 

0.101 

0.624 

6.85 

67.88 

1918 

14.7 

Disturbed,  quadrant  4 

24.9  (0.020) 

0.056 

0.220 

5.67 

101.66 

1280 

9.8 

97.2  (0.077) 

0.091 

0.347 

7.23 

79.44 

1640 

12.6 

320.2  (0.255) 

0.096 

0.198 

9.77 

101.80 

1277 

9.8 

833.6  (0.664) 

0.055 

0.073 

10.90 

199.08 

654 

5.0 

of  reattachment  of  the  flow  over  an  aft-facing  ramp.  The  two  flows  each  had  a 
nominal  external  velocity  of  16  m/s  and  a  boundary-layer  thickness  of  approximately 
1200 The  data  are  of  high  quality.  Comparisons  of  the  statistical  and  spectral 
results  for  the  equilibrium  flow  with  the  archived  literature  indicate  that  the  sensors 
were  of  sufficiently  small  size  to  resolve  the  majority  of  the  small-scale  turbulent 
structures.  The  scaling  for  the  disturbed  flow  was  even  more  favorable  due  to  its 
larger  physical  boundary-layer  thickness  and  lower  skin  friction.  Considering  time- 
series  analyses,  the  sampling  intervals  were  less  than  one  viscous  time  unit. 

The  primary  purpose  of  this  investigation  is  to  shed  light  on  the  question  of  in¬ 
ner/  outer  layer  coupling  in  turbulent  wall-bounded  flows.  The  developing  boundary 
layer  downstream  of  a  ramp  was  chosen  because  it  exhibits  similarities  to  the  base¬ 
line  equilibrium  case  in  the  near-wall  region,  with  the  disturbances  primarily  limited 
to  the  outer  regions.  In  particular,  the  two  flows  show  the  same  classic  log-law 
velocity  profile  for  y+  <  100,  and  equivalent  u™,  and  profiles  for  y/5  <  0.02 
when  scaled  on  outer-flow  variables.  In  addition,  the  high-frequency  behavior  of 
the  u,  v,  and  p  power  spectra  are  similar  both  in  level  and  rolloff.  Also,  the  wall 
pressure  peak-event  detection  statistics  and  the  central,  high-gradient  portion  of  the 
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Figure  3.31:  Peak  values  of  the  equilibrium  flow  Reynolds  stress  quadrant  events  at 
y/S  =  0.02, /c  =  l. 


event  shapes  all  scale  on  wall  variables.  Taken  together,  these  points  indicate  that 
near-wall,  turbulence  production-driven  processes  are  common  to  both  flows. 

The  most  notable  feature  of  the  nonequilibrium  flow  is  the  presence  of  a  distur¬ 
bance  layer  located  at  y/6  =  0.2.  This  disturbance  layer,  directly  visible  as  peaks 
in  the  urms ,  Vj^g^  and  uv  profiles,  is  a  remnant  of  the  free  shear  layer  formed  at  the 
ramp  trailing  edge.  Its  effects  are  felt  throughout  the  boundary  layer  and  extend 
right  down  to  the  wall,  where  it  manifests  itself  as  a  range  of  /~5/3  behavior  in  the 
u  and  p  power  spectra  and  increased  correlation  of  the  tails  of  wall-pressure  peak 
events. 

Other  differences  between  the  two  flows  were  seen  in  the  makeup  of  the  peak-event 
detection  results.  At  the  wall,  the  disturbed  flow  showed  relatively  more  positive 
peak-pressure  events  than  the  equilibrium  flow.  Away  from  the  wall,  different  flow 
structures  were  indicated  by  the  uv  quadrant  analysis.  The  equilibrium  flow  is 
dominated  by  second  quadrant  events  over  the  entire  boundary  layer.  In  contrast, 
the  disturbed  flow  is  characterized  by  first  and  fourth  quadrant  events  near  the  walfl 
changing  to  second  quadrant  dominance  at  the  location  of  the  disturbance  layer. 

While  the  event-detection  results  for  the  disturbed  boundary  layer  are  distorted 
by  the  large-scale,  low-frequency  motions  introduced  into  the  outer  layer,  the  near¬ 
wall  motions  show  a  basic  similarity  to  the  equilibrium  boundary  layer.  Based  on 
these  observations,  the  hypothesis  is  made  that  the  inner  and  outer  layer  portions 
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Figure  3.32:  Peak  value  of  the  disturbed  flow  Reynolds  stress  quadrant  events  at 
y/5  —  0.02,  k  =  1. 

of  turbulent  boundary  layers  are  uncoupled.  Proof  of  this  hypothesis  will  require 
filtering  of  the  disturbance  “noise”  from  the  disturbed  flow  data.  The  remainder 
of  the  research  will  focus  on  separating  the  disturbances  from  the  common  near¬ 
wall  features  using  wavelet  filtering  techniques.  This  will  allow  a  more  thorough 
examination  of  the  extent  of  the  inner-outer-layer  coupling. 
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Chapter  4 

Application  of  wavelet  filtering 


4.1  Introduction 

The  previous  chapter  introduced  the  concept  of  event  detection,  wherein  the  portions 
of  a  signal  that  satisfy  a  given  criteria  are  noted  for  further  analysis.  Identification 
of  the  desired  event  and  determination  of  the  event  statistics  are  hampered  by  the 
presence  of  “noise”  in  the  data.  Noise  in  this  instance  includes  unwanted  signal  com¬ 
ponents  as  well  as  the  more  common  sources  such  as  imperfect  transducer  response 
and  electrical  line  noise.  For  example,  the  wall  pressure  is  obtained  as  the  integral 
of  turbulent  sources  in  the  entire  velocity  field.  If  one  wishes  to  study  the  pressure 
resulting  from  a  particular  flow  region,  the  influences  from  other  regions  need  to  be 
excluded.  This  is  precisely  what  is  required  to  achieve  the  present  goal  of  assessing 
the  relative  effects  of  the  inner  and  outer  layer  flow  structures  on  the  turbulence 
production  mechanism.  This  chapter  examines  a  filtering  scheme  based  on  wavelet 
transforms  to  accomplish  this  selective  exclusion. 

As  reported  by  Eckelmann  [4.1],  Kobashi  and  Ichijo  [4.2]  used  high-  and  low- 
pass  filters  to  separate  pressure  fluctuations  produced  by  the  boundary-layer  bulges 
from  those  originating  within  the  boundary  layer.  Naguib  and  Wark  [4.3]  discrimi¬ 
nated  between  the  influences  of  different  flow  regions  by  using  FIR  band-pass  filters 
matched  to  the  spectral  signature  of  each  region.  For  the  current  goal  of  learning 
more  about  the  sources  of  turbulent  events,  traditional  Fourier-based  analysis  may 
not  be  ideal.  As  noted  by  Tennekes  and  Lumley  [4.4,  p.  259],  and  more  recently  by 
Zubair  [4.5],  eddies  or  organized  structures  are  associated  with  many  Fourier  coeffi¬ 
cients  and  the  phase  relations  among  them.  Therefore,  any  frequency-based  filtering 
will  affect  all  structures  to  some  degree,  those  of  interest  as  well  as  those  considered 
noise.  A  filtering  scheme  based  on  structure  or  scale  instead  of  frequencies  may  be 
more  appropriate. 

This  chapter  will  briefly  review  wavelet  theory,  develop  a  wavelet  filter,  and  as¬ 
sess  the  effectiveness  of  the  filter  for  turbulent  event-detection  applications.  Wavelet 
transforms  are  a  relatively  recent  development  in  applied  mathematics.  The  subject 
has  been  receiving  a  great  deal  of  attention  in  several  different  fields.  A  comprehen- 
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sive  presentation  of  the  theory  can  be  found  in  Daubechies  [4.6].  Applications  to 
fluid  mechanics  are  covered  in  the  review  article  by  Farge  [4.7]. 


4.2  Wavelet  transforms 

Like  the  Fourier  transform,  the  wavelet  transform  is  obtained  from  the  inner  product 
of  the  signal  under  investigation  with  a  set  of  analyzing  functions.  In  the  case  of  the 
Fourier  transform,  the  analyzing  functions  are  sines  and  cosines  of  infinite  extent. 
The  result  is  a  mapping  of  global  frequency  content.  In  the  case  of  the  wavelet  trans¬ 
form,  the  analyzing  functions  (wavelets)  are  dilations  and  translations  of  a  “mother” 
function  which  itself  is  localized  in  the  functional  domain.  Dilation,  or  contraction, 
results  in  decomposition  of  the  signal  into  scales  or  frequencies.  Translation  enables 
the  retention  of  time  information.  The  limited  spatial  support  of  the  wavelets  means 
that  the  behavior  of  the  signal  at  infinity  plays  no  role.  Thus,  the  wavelet  transform 
is  a  means  of  time-scale  decomposition. 

The  continuous  wavelet  transform  of  a  function  f(t )  is  defined  as 

W/(a,  b)  =  |a|“1/2  J  f(t)xp  j  dt  (4-1) 

where  the  function  xp(t)  is  termed  the  “mother”  wavelet.  The  scale  is  controlled 
through  the  parameter  a  while  the  time  localization  through  b.  As  usual,  if  xp  is 
complex,  its  conjugate  is  used  in  forming  the  inner  product.  The  inverse  transform 
exists,  and  is  given  by 


m  =  C;1J  I  W,(a,b)r-b(t (4.2) 
where  the  xpa,b(t)  —  \a\~1/2xp(^)  are  the  wavelets  and  the  constant  is  given  by 

0^  =  21:  J  (4.3) 

For  the  inverse  transform  to  make  sense,  C ^  must  be  finite,  which  leads  to  the 
admissibility  condition  on  xp: 

J  xp(t)dt  =  0  (4.4) 

For  the  present  investigation,  the  functions  to  be  analyzed  are  discrete  time  series. 
Therefore,  the  discrete  wavelet  transform  is  required.  This  is  obtained  by  replacing 
(a,  b )  with  discrete  values  (m,  n),  as  described  by  Daubechies  [4.6].  The  discretization 
of  the  dilation  parameter  is  straightforward:  a  =  a™.  The  discretization  of  the  time 
parameter  must  be  done  to  cover  the  entire  range.  This  means  that  the  shifts 
depend  on  the  wavelet  scale.  Since  a  is  a  measure  of  the  width  of  the  wavelet, 
this  is  accomplished  by  letting  b  =  nb0a™.  The  values  for  m  and  n  range  over  the 
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integers,  and  ao  >  1  and  bo  >  0  are  constants,  which  in  general  depend  on  the 
particular  wavelet,  ip.  A  function  / (t)  is  then  represented  by  its  wavelet  expansion 
with  coefficients  C'm,n  given  by 

Cm,n  =  J  f  (t)ipm'n (t)dt  (4.5) 

In  the  continuous  case,  the  wavelet  transform  is  highly  redundant.  If  the  wavelet 
functions  ipm’n(t)  are  chosen  such  that  they  form  an  orthogonal  basis,  then  the 
redundancies  are  eliminated.  For  this  case,  the  smallest  scale  will  correspond  to  two 
data  points  and  the  shifts  will  be  integer  multiples  of  the  scale  widths:  a0  =  2  and 
bo  =  1.  The  basis  can  then  be  written  as 

~  n2~m)  (4-6) 

This  leads  to  several  useful  properties  of  the  wavelet  transform.  In  particular,  anal¬ 
ogous  to  Parseval’s  theorem  for  Fourier  transforms,  the  total  power  is  conserved: 

j\m?dt  =  Y,Y.\Cm,n\2  (4.7) 

m  n 

Fast  algorithms  have  been  developed  for  the  computation  of  the  orthogonal 
wavelet  transform  and  its  inverse.  The  algorithm  of  Press,  et  al.  [4.8,  4.9]  has  been 
implemented  for  the  present  study.  This  algorithm  will  be  referred  to  as  the  fast 
wavelet  transform  (FWT).  Like  the  fast  Fourier  transform,  the  time  series  under 
analysis  must  have  a  length  that  is  an  integer  power  of  two.  If  the  total  number  of 
samples  is  2N,  then  the  wavelet  expansion  will  also  have  2N  coefficients,  Ce,  with 
0  <  £  <  2n  -  1.  The  index  £  is  related  to  the  scale/time  indices  ( m,n )  by  the 
following  relationships 

_  /  N  +  l,  if  ^  =  0 

m  j  iV  —  [log2  £\ ,  if  £  >  0  (48) 

_  /  1,  if  £  =  0 

n  [  £  -  2^*1  +  1,  if  £  >  0 

where  [zj  is  the  largest  integer  not  exceeding  x.  This  represents  a  dyadic  pattern 
in  the  (m,n)  plane,  as  illustrated  in  Figure  4.1.  The  expansion  is  performed  at  N 
different  scales,  with  each  scale  represented  by  a  row  in  the  grid.  For  each  row  m, 
there  are  2N  m  values  of  the  shift  index,  n.  The  index  i  =  0  corresponds  to  what 
is  known  as  the  scaling  function;  it  is  roughly  analogous  to  the  dc  component  in  a 
Fourier  decomposition.  The  index  £  =  2N  -  1  corresponds  to  the  smallest  scale  in 
the  analysis  located  at  the  end  of  the  record. 

Up  to  this  point,  the  choice  of  wavelet  has  not  been  discussed.  Families  of 
orthogonal  wavelets  with  compact  support  have  been  defined  in  the  literature.  A 
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m  =  N  -  [log2£j 
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Figure  4.1:  Dyadic  arrangement  of  wavelet  coefficients,  N  —  5. 

particular  family  is  described  by  a  set  of  coefficients,  called  wavelet  filter  coefficients, 
which  is  recursively  applied  to  the  data.  Guidelines  for  choosing  a  wavelet  family 
are  given,  for  example,  by  Farge  [4.7].  The  choice  can  also  be  automated,  as  is 
done  in  the  “wavelet-packet”  transforms  utilized  by  Zubair  [4.5].  In  this  method, 
the  family  is  selected  from  a  library  of  basis  functions  according  to  a  minimum 
entropy  criterion.  The  family  so  chosen  provides  the  “best  representation”  of  the 
data  from  the  standpoint  of  data  compression.  Data  compression  is  discussed  further 
in  Section  4.4.1. 

Several  wavelet  families  have  been  studied  for  the  present  analysis  [4.10,  4.11]. 
The  conclusion,  based  on  extensive  computations,  was  that  turbulent  event-detection 
results  are  robust  with  respect  to  the  families  examined.  Therefore,  one  family  of 
wavelets  was  selected  and  used  in  all  subsequent  work,  namely,  the  Daubechies 
extremal  phase  set  with  eight  coefficients  [4.6,  p.  195].  This  family  will  be  referred 
to  as  the  D8  set.  The  shapes  of  the  wavelets  for  four  different  scales  and  locations  are 
illustrated  in  Figure  4.2.  This  curve  was  generated  by  performing  the  inverse  wavelet 
transform  of  a  vector  of  arbitrary  length  1024  ( N  =  11)  consisting  of  all  zeroes, 
except  elements  19,  48,  110,  and  243,  which  were  set  equal  to  one.  These  elements 
represent  unit  vectors  in  the  wavelet  bases.  In  the  dyadic  grid,  they  are  located  in 
rows  6,  5,  4,  and  3,  respectively.  Element  19  is  early  in  the  range  16  <  i  <  31,  row 
m  =  6,  resulting  in  a  wavelet  positioned  early  in  time,  while  element  48  is  in  the 
middle  of  the  range  32  <  £  <  63,  row  m  =  5,  resulting  in  a  wavelet  positioned  near 
the  middle  of  the  sequence,  and  so  forth.  Notice  how  each  successive  wavelet  is  one 
half  as  wide  as  its  predecessor,  with  an  increasing  amplitude  such  that  the  energy  is 
constant. 


69 


NSWCDD/TR-95/222 


Figure  4.2:  Four  different  scales  of  the  Daubechies  extremal  phase  wavelet  with  8 
coefficients. 

4.3  Scale- frequency  relations  and  wavelet  power 

The  fast  wavelet  transform  performs  a  time-scale  decomposition  of  a  signal  in  terms 
of  a  set  of  orthogonal  basis  functions.  Before  constructing  a  wavelet-based  filter,  two 
questions  naturally  arise: 

1.  What  is  the  relationship  between  wavelet  scale  and  conventional  frequency? 

2.  What  is  an  appropriate  wavelet  equivalent  to  the  Fourier  power  spectral  den¬ 
sity? 

The  first  question  can  be  addressed  by  performing  the  Fourier  transform  of  the 
wavelets  at  fixed  scales.  This  was  done  by  forming  a  wavelet  from  a  unit  vector 
as  seen  in  Figure  4.2  and  taking  the  fast  Fourier  transform  (FFT).  The  conversion 
to  Hertz  depends  on  the  sampling  interval,  which  was  equal  to  30.5  ns  for  this 
investigation.  The  results  are  illustrated  in  Figure  4.3  for  the  four  smallest  scales. 
Sequential  wavelet  scales  are  seen  to  contain  overlapping  frequency  bands.  However, 
each  scale  can  be  characterized  by  the  frequency  where  the  FFT  magnitude  is  a 
maximum.  This  is  termed  the  center  frequency,  fc,  and  is  tabulated  in  Table  4.1  for 
1  <  m  <  11.  Notice  how  the  smallest  scale  does  not  have  a  peak  in  the  power  spectral 
density  (PSD),  but  rather  a  wide  plateau  extending  to  the  Nyquist  frequency.  This 
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Figure  4.3:  Power  spectrum  magnitude  of  the  D8  wavelets  at  several  scales  for  a 
sampling  rate  of  32768  Sa/s  (At,  =  30.5  ps ). 

exercise  demonstrates  the  trade-offs  between  localization  in  frequency  of  the  Fourier 
analysis  and  localization  in  scale  of  the  wavelet  analysis. 

Analogous  to  the  Fourier  power  spectral  density  is  the  wavelet  power  scalar 
density,  abbreviated  as  PScD.  This  can  be  constructed  for  a  general  function  / 
by  summing  the  square  magnitude  of  the  wavelet  coefficients  in  each  scale,  m: 

>)  =  Et  (4.10) 

n 

Power  scalar  densities  of  the  equilibrium  flow  u,  v,  and  p  time  records  are  shown 
in  Figure  4.4.  The  majority  of  the  power  for  each  signal  is  located  in  the  scales 
2  <  m  <  8.  The  second  peak  in  the  wall  pressure  at  scale  m  —  7,  centered  at 
179  Hz,  corresponds  to  extraneous  noise  caused  by  the  traverse  system  disturbances 
as  seen  in  Figure  3.18. 

4.4  Wavelet  filtering 

Filtering  based  on  the  orthogonal  wavelet  expansion  can  be  accomplished  by  modifi¬ 
cation,  in  some  rational  way,  of  the  individual  wavelet  coefficients  [4.5,  4.8,  4.9].  For 
example,  setting  a  particular  coefficient  Cm>n  equal  to  zero  would  eliminate  scale  m 
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Table  4.1:  D8  wavelet  scales  and  approximate  center  frequencies  for  a  sampling  rate 
of  32768  Sa/s  (A ts  =  30.5  /rs). _ 


Scale  (row) 
m 

Shift  a  =  2m 
(samples) 

fc 

(Hz) 

1 

2 

16384 

2 

4 

5719 

3 

8 

2858 

4 

16 

1429 

5 

32 

714 

6 

64 

357 

7 

128 

179 

8 

256 

89 

9 

512 

45 

10 

1024 

22 

11 

2048 

11 

Figure  4.4:  D8  power  scalar  densities  of  the  equilibrium  flow  u,  v,  and  p  time  records 
at  y/5  =  0.02. 


from  the  location  n.  Due  to  the  finite  support  of  the  wavelet  basis  functions,  other 
scales  and  locations  remain  unaffected. 
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index 


Figure  4.5:  Normalized  D8  FWT  magnitude  of  a  segment  of  the  equilibrium  flow  u 
velocity  at  y/5  =  0.02. 

4.4.1  Conventional  wavelet  filtering 

A  first-order  wavelet  filter  can  be  implemented  by  setting  to  zero  all  coefficients 
with  magnitude  less  than  some  threshold  [4.5,  4.8-4.10].  The  threshold  can  be  set  in 
terms  of  the  fraction  of  the  total  signal  energy  to  be  retained  in  the  filtered  sequence, 
denoted  by  kf.  First,  the  signal  is  transformed  and  the  resulting  wavelet  coefficients 
arranged  in  decreasing  magnitude  order.  The  cumulative  energy  is  computed  until 
it  equals  the  desired  energy  fraction.  All  coefficients  beyond  this  point  are  then 
set  equal  to  zero.  To  construct  the  filtered  sequence,  the  modified  coefficients  are 
returned  to  their  original  order  and  inverse  wavelet  transformed.  For  a  general 
function  /,  the  wavelet-filtered  sequence,  fwf,  will  have  energy 

E4,=*/£/2  (4-n) 

This  filtering  scheme  will  reject  those  structures  contributing  the  least  to  the 
signal  energy.  The  effects  can  be  predicted  by  looking  at  the  magnitude  of  the 
wavelet  coefficients  of  a  signal  of  interest.  This  is  shown  in  Figure  4.5  for  a  32768- 
point  segment  of  a  typical  u  velocity  trace.  The  normalized  magnitude  of  the  D8 
wavelet  coefficients  are  plotted  versus  index,  i.  Elimination  of  coefficients  would 
generally  begin  with  the  largest  scales  t  <  31,  then  continue  with  the  smallest  scales 
t  >  16384,  etc.  Of  course,  some  coefficients  at  all  scales  will  be  eliminated  based 
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Figure  4.6:  Segment  of  equilibrium  u  velocity  before  and  after  conventional  wavelet 
filtering  with  kf  =  0.99  and  0.90. 


on  their  magnitude.  Since  the  signal  has  a  zero  mean,  very  little  energy  is  in  the 
largest  scales.  Therefore,  most  of  the  eliminated  energy  will  come  from  the  smallest, 
high-frequency  scales.  The  conventional  wavelet  filter  will  therefore  appear  similar 
to  a  Fourier-based,  low-pass  filter. 

The  conventional  wavelet  filter  highlights  the  data  compression  capabilities  of 
the  wavelet  transform.  This  is  illustrated  by  performing  D8  conventional  wavelet 
filtering  on  the  equilibrium  flow  u  velocity  at  location  y/5  =  0.02.  A  filtered  sequence 
with  kf  =  0.99  can  be  reconstructed  using  only  17%  of  the  wavelet  coefficients. 
This  represents  a  6:1  compression  ratio  with  only  a  1%  loss  of  information.  If 
10%  of  the  original  energy  is  filtered  out  (kf  =  0.90),  then  only  5%  of  the  wavelet 
coefficients  are  needed  for  reconstruction.  These  high  compression  ratios  confirm 
that  the  D8  wavelets  are  an  appropriate  basis  for  turbulence  records,  as  conjectured 
at  the  beginning  of  this  chapter.  Figure  4.6  shows  these  results  in  the  time  domain. 
The  low-pass  characteristics  of  the  filter  are  evident  in  the  kf  =  0.90  case.  Small-scale 
(high-frequency)  fluctuations  have  been  filtered  out,  while  larger  scale  structures 
have  been  faithfully  reproduced  with  no  phase  distortions. 
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4.4.2  Weighted  wavelet  filtering 

In  the  preceding  section,  it  was  seen  that  conventional  wavelet  filtering  is  similar 
to  low-pass  filtering  for  typical  turbulence  quantities.  A  more  sophisticated  filter 
can  be  constructed  by  explicitly  taking  scale  into  account  in  the  modification  of 
coefficients.  Several  rational  schemes  for  weighting  the  coefficients  before  selection 
were  investigated  [4.10,  4.11].  For  the  present  application,  the  approach  that  proved 
to  be  appropriate  was  to  first  decide  which  scales  mT  would  be  retained  in  the  filtered 
sequence,  and  then  assign  an  energy  fraction  for  each  retained  scale.  In  terms  of 
conventional  wavelet  filtering, 


kf(m),  if  m  €  mT 
0,  otherwise 


(4.12) 


where  mT  designates  the  discrete  set  of  retained  scales. 

Selecting  which  rows  to  retain  is  analogous  to  selecting  cutoff  frequencies  in 
Fourier  filtering.  This  was  done  based  on  the  physical  features  of  the  flows  and 
numerous  computational  experiments,  with  the  major  consideration  being  the  dis¬ 
turbed  flow  data.  Referring  to  the  wall-pressure  spectra  in  Figure  3.18,  the  dis¬ 
turbed  flow  exhibits  excess  energy  in  frequencies  below  approximately  3000  Hz. 
According  to  Table  4.1  and  Figure  4.3,  this  range  corresponds  to  scales  m  >  4, 
which  will  be  eliminated.  The  smallest  scale  m  =  1  corresponds  primarily  to  fre¬ 
quencies  above  104  Hz,  and  contains  very  little  energy,  as  seen  in  Figure  4.4.  It 
will  also  be  eliminated,  leaving  mr  =  { 2  3  4}  for  the  wall  pressure.  The  energy 
fractions  kf(mr )  =  {0.99  0.99  0.60}  were  selected.  The  value  kf( 4)  =  0.60  was 
chosen  since  scale  m  =  4  is  mostly  below  the  spectral  range  of  interest.  The  values 
kf(  2)  =  kf(  3)  =  0.99  were  arbitrarily  chosen  to  provide  some  level  of  data  compres¬ 
sion. 

The  u  and  v  velocities  were  considered  together.  The  goal  is  to  separate  passive 
disturbances  from  active  turbulence  production;  i.e.,  the  Reynolds  stress  —  uv.  Be¬ 
cause  of  the  increase  in  axis  crossings,  the  uv  product  will  exhibit  higher-frequency 
components  than  either  u  ox  v  separately.  Therefore,  the  scales  for  u  and  v  were 
increased  by  one  relative  to  the  wall  pressure:  mr  =  { 3  4  5}  along  with  kf(mr )  = 
{0.99  0.99  0.99}.  This  set  of  filtering  parameters  was  designated  the  F2  wavelet 
filter,  and  is  summarized  in  Table  4.2. 

The  effects  of  the  F2  wavelet  filter  in  terms  of  the  wavelet  power  (Figure  4.4) 
are  obvious;  the  power  in  each  scale  will  be  multiplied  by  kf(mr).  In  the  frequency 
domain,  the  effect  is  shown  in  Figures  4.7  and  4.8  for  equilibrium  flow  u  and  p 
records.  The  u-velocity  case  is  not  shown,  since  the  filtering  parameters  for  u 
and  v  are  identical.  The  power  spectra  were  computed  using  Welch’s  averaged 
periodogram  method  with  a  window  length  of  1024  points  and  50%  overlap.  At 
first  glance,  the  results  appear  similar  to  a  band-pass  Fourier  filter.  However,  there 
are  several  differences.  The  rolloffs  on  the  low-frequency  side  are  not  as  sharp  as 
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Table  4.2:  The  F2  wavelet,  filter. 


D8  discrete  orthogonal  wavelet  expansion 

Variable 

Scales  retained 

Power  fractions 

mr 

kf(mr ) 

u 

3 

0.99 

4 

0.99 

5 

0.99 

V 

3 

0.99 

4 

0.99 

5 

0.99 

p 

2 

0.99 

3 

0.99 

4 

0.60 

10 

Frequency  (Hz) 


Figure  4.7:  Power  spectra  of  the  equilibrium  flow  u  velocity  before  and  after  F2 
wavelet  filtering. 


typical  Fourier  filters,  especially  for  the  wall  pressure  where  the  value  kf( 4)  =  0.60 
introduces  an  even  smoother  transition.  On  the  high-frequency  side,  a  single  cutoff 
frequency  is  not  evident  in  either  case.  These  features  occur  because  wavelet  scales 
are  not  localized  in  frequency,  as  was  shown  in  Figure  4.3.  Indeed,  these  features 
are  difficult  to  replicate  using  conventional  band-pass  filters. 
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Figure  4.8:  Power  spectra  of  the  equilibrium  flow  wall  pressure  before  and  after  F2 
wavelet  filtering. 


The  effects  of  the  filter  on  segments  of  equilibrium  and  disturbed  flow  u,  v,  and 
p  time  records  are  illustrated  in  Figures  4.9  and  4.10.  The  features  of  the  signals  at 
the  selected  scales  are  left  intact.  For  example,  the  large  amplitude  events  around 
the  time  6  ms  in  Figure  4.9  are  reproduced,  while  others  are  filtered  out.  This 
emphasizes  the  ability  of  the  wavelet  transform  to  highlight  transients  in  a  signal. 
Spectral  filtering,  on  the  other  hand,  would  attenuate  all  events  at  all  locations  in 
time. 


4.5  Comparison  of  results 

To  evaluate  the  F2  wavelet  filter,  the  peak  wall  pressure  and  quadrant  event  detection 
analyses,  pursued  for  the  unfiltered  data  and  discussed  in  the  previous  chapter,  were 
repeated  for  the  wavelet-filtered  equilibrium  and  disturbed  flow  data.  In  addition, 
conditional  sampling  was  performed  for  both  unfiltered  and  F2-filtered  data.  The 
aim  of  the  conditional  sampling  studies  was  to  examine  the  bidirectional  relationship 
between  wall  pressure  and  near-wall  organized  motions  [4.12-4.14],  and  to  assess  the 
ability  of  the  F2  filter  to  preserve/enhance  this  relationship  in  the  equilibrium  and 
disturbed  flow  data. 

Quantities  that  have  been  processed  by  the  F2  filter  will  be  denoted  with  a 
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Figure  4.9:  Effects  of  F2  wavelet  filtering  on  equilibrium  flow  u,  v,  and  pi  at  y/6  — 
0.02. 

subscript  F2.  In  keeping  with  the  length  requirements  of  the  FWT  algorithm,  the 
first  218  =  262144  points  were  used.  This  represents  8  s  of  the  10  s  recorded. 
Thresholds  were  defined  with  respect  to  the  filtered  rms  levels. 

4.5.1  Peak  and  quadrant  event  detection 

Table  4.3  summarizes  the  results  of  peak-event  detection  performed  on  the  equilib¬ 
rium  and  disturbed  wall  pressure  after  F2  filtering.  Comparing  with  Table  3.4,  the 
numbers  of  events  detected  are  increased,  resulting  in  higher  event  frequencies,  in 
spite  of  the  shorter  time  record  examined.  This  is  because  the  rms  level  has  been 
reduced  by  the  filtering,  but  the  amplitude  of  the  detected  events  has  not.  The  aver¬ 
age  event  durations  A t+  have  been  reduced  by  factors  of  3  to  6,  with  corresponding 
reductions  in  the  time  fractions.  This  indicates  that  large-scale  structures  have  been 
filtered  out.  Whereas  wall  variable  scaling  collapsed  the  equilibrium  and  disturbed 
data  before  F2  filtering,  that  is  no  longer  the  case.  Since  identical  filtering  parame¬ 
ters  were  applied  to  both  equilibrium  and  disturbed  data,  large  amplitude  events  of 
the  same  physical  scale  (time)  were  retained.  The  difference  in  the  statistics  for  the 
filtered  data  is  therefore  due  to  the  difference  in  ur  (cf.  Table  3.1).  Average  event 
shapes  will  be  presented  in  the  next  section  with  the  results  of  flowfield  ensemble 
averaging  conditioned  on  the  peak  wall-pressure  events. 
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Time  (s) 


Figure  4.10:  Effects  of  F2  wavelet  filtering  on  disturbed  flow  u,  v,  and  pi  at  y/8  — 
0.02. 

Recall  that  though  the  wall-pressure  event  statistics  collapsed  before  F2  filter¬ 
ing,  the  distributions  of  the  time  between  events  were  of  different  character  (Fig¬ 
ures  3.23  and  3.24).  This  has  been  reversed  with  F2  filtering,  as  shown  in  Figures  4.11 
and  4.12.  The  disturbed  flow  case  now  shows  a  log-normal  distribution,  similar  to 
the  equilibrium  case,  in  place  of  the  more  uniform  distribution  prior  to  F2  filtering. 

The  event  statistics  for  the  F2-filtered  turbulence-producing  quadrant  events 
across  the  boundary  layers  are  presented  in  Table  4.4.  Most  of  the  observations 
concerning  the  wall- pressure  events  also  apply  to  the  quadrant  events.  The  frequency 
of  detected  events  has  increased,  while  the  average  event  durations  and  time  fractions 
have  decreased,  again  implying  that  large-scale  structures  have  been  filtered  out.  The 
quadrant-event  frequencies  across  the  boundary  layers  are  shown  in  Figures  4.13 
and  4.14.  Comparing  with  Figures  3.29  and  3.30,  the  equilibrium  case  is  relatively 
unchanged  by  the  filtering;  Q2  events  still  dominate  everywhere  followed  by  Q4 
events.  For  the  disturbed  case,  Ql  and  Q4  events  still  dominate  in  the  inner  layer, 
and  the  changeover  to  Q2  dominance  still  occurs  at  y/8  &  0.2.  However,  the  filtering 
has  increased  the  frequency  of  Q2  events  relative  to  the  other  quadrants  for  y/8  < 
0.2. 
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Table  4.3:  Wall  pressure  peak-event  statistics  for  transducer  pi  after  F2  wavelet 
filtering. _ 


K 

Time 

fraction 

Power 

fraction 

A ,t+ 

AT+ 

Number  of 

events 

Frequency 
per  1000t+ 

Equilibrium 

+1  0.161 

0.453 

5.5 

34.3 

6111 

29.2 

-1 

0.161 

0.463 

5.5 

33.8 

6191 

29.5 

+2 

0.067 

0.348 

6.2 

92.7 

2261 

10.8 

-2 

0.068 

0.359 

6.1 

89.9 

2330 

11.1 

+3 

0.026 

0.223 

6.3 

243.1 

862 

4.1 

-3 

0.026 

0.235 

6.0 

232.9 

900 

4.3 

Disturbed 

+1  0.106 

0.475 

3.1 

29.1 

3586 

34.3 

-1 

0.108 

0.468 

3.1 

28.8 

3623 

34.7 

+2 

0.059 

0.423 

3.5 

59.5 

1755 

16.8 

-2 

0.059 

0.415 

3.6 

60.5 

1724 

16.5 

+3 

0.031 

0.336 

3.6 

116.6 

895 

8.6 

-3 

0.030 

0.324 

3.7 

122.1 

855 

8.2 
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Figure  4.11:  Distribution  of  the  time  between  equilibrium  flow  pwf  peak  events, 
k  =  1. 
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Table  4.4:  Reynolds  stress  quadrant-event  detection  statistics  after  F2  wavelet  fil- 


tering,  k  =  —  2. 


y+  (y/S) 

Time 

fraction 

Power 

fraction 

A  t+ 

AT+ 

Number  of 

events 

Frequency 
per  1000f+ 

Equilibrium,  quadrant  2 

25.9  (0.022) 

0.073 

0.435 

8.14 

111.24 

1885 

9.0 

93.1  (0.077) 

0.075 

0.414 

7.22 

96.37 

2175 

10.4 

300.0  (0.250) 

0.068 

0.327 

8.08 

119.34 

1754 

8.4 

781.1  (0.650) 

0.067 

0.390 

8.06 

120.59 

1737 

8.3 

Equilibrium,  quadrant  4 

25.9  (0.022) 

0.054 

0.182 

9.38 

174.75 

1199 

5.7 

93.1  (0.077) 

0.063 

0.215 

7.90 

125.23 

1674 

8.0 

300.0  (0.250) 

0.062 

0.252 

8.44 

135.54 

1543 

7.4 

781.1  (0.650) 

0.060 

0.236 

9.01 

149.18 

1403 

6.7 

Disturbed,  quadrant  2 

24.9  (0.020) 

0.037 

0.167 

3.98 

108.35 

965 

9.2 

97.2  (0.077) 

0.049 

0.214 

4.27 

86.91 

1200 

11.5 

320.2  (0.255) 

0.063 

0.309 

3.99 

63.72 

1639 

15.7 

833.6  (0.664) 

0.064 

0.431 

3.85 

60.60 

1724 

16.5 

Disturbed,  quadrant  4 

24.9  (0.020) 

0.037 

0.197 

3.64 

99.11 

1054 

10.1 

97.2  (0.077) 

0.051 

0.249 

3.87 

75.71 

1381 

13.2 

320.2  (0.255) 

0.060 

0.247 

4.24 

71.07 

1467 

14.0 

833.6  (0.664) 

0.055 

0.222 

4.53 

83.13 

1258 

12.0 

4.5.2  Conditional  sampling 

To  assess  the  characteristic  features  of  intermittent  flow  events,  and  their  correlation 
with  wall-pressure  events,  all  flowfield  variables  can  be  ensemble-averaged  over  a 
time  window  based  on  a  detection  criteria  for  the  occurrence  of  that  event.  This 
is  known  as  conditional  sampling.  Yuan  and  Mokhtarzadeh-Dehghan  [4.15]  have 
recently  reviewed  and  compared  many  of  the  methods  found  in  the  literature.  In 
this  study,  conditional  sampling  was  performed  based  on  the  detection  of  peak  wall- 
pressure  events  and  based  on  the  detection  of  quadrant  events.  The  data  will  be 
presented  at  the  location  y/5  —  0.02.  Ensemble  averages,  designated  by  (  ),  were 
computed  for  p,  u,  v,  and  the  quadrant  functions  uv2  and  uv4,  given  by 


uv2(t)  = 


u(t)v(t),  if  u  <  0  and  v  >  0 
0,  otherwise 


(4.13) 


«.<*)  =  {  W()' 


if  u  >  0  and  v  <  0 
otherwise 


(4.14) 
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Figure  4.14:  Disturbed  flow  quadrant-event  frequency  after  F2  wavelet  filtering, 
k  =  2. 

The  averages  were  centered  about  the  points  of  peak  magnitude  in  the  detection 
variable.  In  the  results  to  follow,  the  detection  variable  is  plotted  using  a  thicker 
solid  line. 

Conditional  sampling  on  p 

Averages  conditioned  on  the  detection  of  peak  wall-pressure  events  are  considered 
first.  The  equilibrium  flow  results  are  presented  in  Figures  4.15  and  4.16  for  the 
unfiltered  and  F2-filtered  cases,  respectively.  The  wall-pressure  traces  have  been 
shifted  upwards  for  clarity.  The  most  notable  difference  is  the  shape  of  the  wall- 
pressure  event.  The  filtering  has  exaggerated  the  pre-  and  post-event  dips.  This 
effect  is  responsible  for  the  shorter  average-event  durations.  In  both  cases,  significant 
correlation  is  seen  in  the  flowfield  variables.  The  u  velocity  increases  while  the  v 
velocity  decreases,  moving  from  Q2  to  Q4.  The  trends  are  somewhat  intensified  after 
filtering. 

The  disturbed  flow  results  are  shown  in  Figures  4.17  and  4.18.  In  the  unfiltered 
case,  the  v  velocity  remains  flat  over  the  first  half  of  the  pressure  event,  with  a  sharp 
decrease  after  the  pressure  peak.  Relative  to  the  equilibrium  results,  this  produces 
no  correlation  with  Q2  events  and  exagerated  correlation  with  Q4  motions.  However, 
the  F2  filter  recovers  the  v  and  uv2  correlations  seen  in  the  filtered  equilibrium  flow 
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Figure  4.17:  Ensemble  averages  conditioned  on  peak  wall-pressure  events  with  k  =  2. 
Disturbed  flow  at  y/S  =  0.02  before  filtering. 


Figure  4.18:  Ensemble  averages  conditioned  on  peak  wall-pressure  events  with  k  =  2. 
Disturbed  flow  at  y/S  =  0.02  after  F2  filtering. 
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data.  The  structures  are  narrower  when  normalized  by  viscous  time  units,  for  the 
reasons  noted  previously.  The  filtering  appears  to  capture  the  generic  features  of 
the  near-wall  structure. 

Conditional  sampling  on  uv2 

Conditional  sampling  results  based  on  the  Q2  Reynolds  stress  events  are  shown  in 
Figures  4.19  and  4.20  for  the  equilibrium  flow  and  Figures  4.21  and  4.22  for  the 
disturbed  flow.  The  equilibrium  flow  results  confirm  the  p-uv2  correlation  shown  in 
the  results  based  on  peak  wall  pressure.  That  is,  Q2  motions,  where  the  flow  slows 
and  moves  away  from  the  wall,  coincide  with  local  increases  in  the  wall  pressure. 
Like  the  peak  wall-pressure  results,  F2  filtering  enhances  this  correlation;  the  wall- 
pressure  signature  is  stronger  than  for  the  unfiltered  data.  Also  similar  to  the  peak 
wall-pressure  results,  F2  filtering  produces  pre-  and  post-peak  oscillations  in  the 
detection  variables  u  and  v. 

Looking  at  the  quadrant  detection  results  for  the  disturbed  flow  in  Figures  4.21 
and  4.22,  it  can  be  seen  that  before  filtering,  the  wall-pressure  signature  appears 
primarily  as  a  decrease  prior  to  the  Q2  event,  with  an  increase  back  to  a  zero  level 
at  the  time  of  peak  Reynolds  stress.  The  F2  filter,  however,  recovers  the  behavior 
seen  in  the  equilibrium  results.  Again,  the  ensemble  averaged  structures  are  made 
narrower  by  filtering. 
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4.6  Summary 

The  F2  wavelet  filter  preserves  a  range  of  physical  scales.'  The  scales  retained  are 
analogous  to  the  selection  of  cutoff  frequencies  in  a  Fourier  bandpass  filter.  For 
the  present  application,  the  scales  were  chosen  to  be  representative  of  near- wall 
structures  in  equilibrium  turbulent  boundary  layers. 

Four  sets  of  databases  were  examined:  unfiltered  equilibrium,  unfiltered  dis¬ 
turbed,  filtered  equilibrium,  and  filtered  disturbed.  Two  issues  warrant  considera¬ 
tion.  First,  for  the  equilibrium  flow  data,  what  are  the  effects,  if  any,  of  wavelet 
filtering?  Second,  for  the  disturbed  flow  data,  have  the  generic  features  of  the  near¬ 
wall  flow  structures  been  captured  by  the  filtering  strategy?  The  latter  question  is 
based  on  the  hypothesis  that  the  near- wall  turbulent  structure  is  weakly  affected  by 
the  outer-layer  disturbances,  as  evident  in  the  data  presented  in  Chapter  3. 

The  following  observations  answer  these  questions  with  regard  to  the  performance 
of  the  F2  wavelet  filter: 

•  F2  filtering  retained  the  essential  features  of  the  near-wall  motions  in  the  equi¬ 
librium  flow  and  recovered  these  features  to  a  limited  extent  in  the  disturbed 
flow  results.  By  reducing  overall  energy  levels  and  retaining  only  certain  scales, 
F2  filtering  increases  the  number  and  decreases  the  duration  of  detected  events. 
However,  the  log-normal  distribution  of  the  time  between  peak  wall-pressure 
events  was  preserved  for  the  equilibrium  flow  and  recovered  for  the  disturbed 
flow.  In  the  same  way,  the  relative  contributions  of  the  four  quadrants  to  the 
Reynolds  stress  was  preserved  for  the  equilibrium  flow,  while  the  disturbed 
flow  results  were  improved  over  the  unfiltered  data. 

•  F2  filtering  recovered  the  bidirectional  relationship  between  wall  pressure  and 
Q2  uv  events  in  the  disturbed  flow.  The  correlation  between  wall  pressure  and 
Q2  flow  events,  clearly  seen  in  the  equilibrium  boundary  layer,  is  not  evident 
in  the  original  disturbed  flow  data.  F2  filtering  recovers  this  link.  This  is  an 
important  finding,  which  will  be  expanded  upon  in  the  next  chapter.  The  fact 
that  the  recovered  structures  correspond  to  the  scales  retained  in  the  wavelet 
expansion  does  not  account  for  the  alignment  of  these  structures  in  time.  This 
alignment  is  indicative  of  the  physical  processes  occurring  within  the  turbulent 
boundary  layer. 

In  light  of  these  findings,  the  F2  filter  will  be  employed  for  the  remainder  of  the 
investigation,  and  results  can  be  assumed  to  have  been  obtained  from  F2-filtered 
data  unless  specified  otherwise. 
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Chapter  5 

Detection  of  event  cluster  patterns 


5.1  Introduction 

The  conditional  sampling  results  of  the  previous  chapter  showed  a  correlation  be¬ 
tween  wall-pressure  events  and  flowfield  structure  in  the  form  of  accelerating  u- 
velocity  or  second-quadrant  uv  events.  A  more  careful  examination  of  the  results 
conditioned  on  peak  wall  pressure  in  Figures  4.15-4.18  reveals  that  the  streamwise 
velocity  is  accelerating  from  negative  to  positive  values  while  the  wall-normal  ve¬ 
locity  is  decelerating  from  positive  to  negative.  In  other  words,  a  motion  from  the 
second  to  the  fourth  u-v  quadrant,  or  a  Q2-Q4  event  pair.  Looking  at  the  filtered 
results  conditioned  on  Q2  events  in  Figures  4.20  and  4.22,  Q4  motions  are  evident  on 
either  side  of  the  primary  Q2  event;  i.e.,  a  Q4-Q2-Q4  event  triple.  These  observations 
suggest  that  detection  of  flowfield  structure  based  on  event  groupings,  or  clusters, 
may  provide  a  clearer  picture  than  detection  based  on  single  events.  Therefore,  an 
algorithm  is  required  to  identify  clusters  of  events. 

Another  issue  merits  attention,  namely  event  magnitude.  Event-detection  meth¬ 
ods  identify  those  portions  of  a  signal  that  satisfy  the  detection  criteria,  which  usually 
includes  comparing  some  signal  property  against  a  threshold  value.  For  example, 
the  peak  and  quadrant  event-detection  methods  require  that  the  signal  amplitude 
exceed  a  certain  level.  Another  common  technique,  the  VITA  method,  requires  that 
the  local  variance  of  the  signal  exceed  a  given  threshold.  The  threshold  therefore  dis¬ 
criminates  events  based  on  magnitude  or  strength.  Selection  of  the  threshold  value 
is,  however,  somewhat  arbitrary,  being  based  primarily  on  experience.  Also,  these 
algorithms  ignore  the  fact  that  the  strength  of  an  event  is  not  solely  dependent  on 
the  amplitude  of  the  phenomena  causing  it,  but  also  on  the  relative  spatial  locations. 

The  above  considerations  demonstrate  the  need  for  a  detection  algorithm  that 
can  capture  cluster  patterns  of  events  without  prescribing  an  amplitude  threshold 
level.  The  idea  of  event  clusters  is  not  entirely  new.  Wilczynski  [5.1],  also  using 
peak-event  detection,  observed  Q2-Q4  event  pairs  in  an  equilibrium  TBL,  while 
Guezennec  [5.2]  showed  that  the  VITA  technique,  with  appropriate  alignment,  de¬ 
tected  the  shear  layer  between  Q2-Q4  flow  patterns.  However,  the  development  of 
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detection  algorithms  based  on  the  cluster  concept  is  apparently  without  precedent. 
This  chapter  will  examine  two  methods  to  identify  the  groupings: 

1.  Localized  Variance  Technique. 

2.  Localized  Windowed  Peak  Detection  Technique. 

5.2  Detection  techniques 

5.2.1  Localized  variance  technique 

As  a  first  attempt  at  identifying  event  clusters,  a  variation  of  the  VITA  method 
mentioned  above  was  explored.  The  VITA  method  [5.3,  5.4]  computes  the  local 
variance  of  a  signal,  which  for  an  arbitrary  function  f(t)  is  defined  as 

V{/}  =  V/(t,r.)  =  -  [  f(r)dr  -  f(T)d,T  (5.1) 

where  Tv  is  an  adjustable  short  integration  time  that  can  be  interpreted  as  the  typical 
duration  of  the  event.  An  event  is  detected  when 

V,(t,Tv)  >  Kfi,,  (5.2) 

The  VITA  method  is  almost  always  applied  to  the  streamwise  velocity.  Since 
burst  events  characterized  by  large  Q2  and  Q4  events  are  of  primary  interest,  the 
VITA  approach  was  applied  to  the  uv(t)  =  u(t)v(t)  time  record.  Since  it  is  also 
desired  to  retain  the  pattern  of  these  burst  events  and  to  correlate  them  with  the 
wall-pressure  events,  contiguous  V  time  records  for  uv(t)  and  p(t)  were  computed. 
Thus,  the  local  variance  records  Vp(t)  and  Vuv(t)  were  generated  from  the  original 
time  records.  The  application  of  the  VITA  approach  to  uv  and  p  in  this  way  is 
referred  to  as  the  localized  variance  (LVAR)  technique. 

The  VITA  method  locates  regions  of  large  time  derivative  whose  durations  are 
approximately  equal  to  the  short-time  averaging  period  Tv  [5.4],  In  this  sense,  VITA 
is  reminiscent  of  a  wavelet  filter  with  fixed  scale.  For  this  study,  the  value  Tv  = 
12  samples  was  chosen  to  reflect  the  scales  retained  in  the  F2  wavelet  filter  (cf. 
Tables  4.1  and  4.2).  In  nondimensional  terms,  this  equates  to  T+  =  9.6  for  the 
equilibrium  flow,  which  is  in  the  range  of  values  used  by  other  researchers  [5. 5-5. 7]. 
For  the  disturbed  flow,  T+  =  4.8. 

The  correspondence  of  the  local  variance  with  Tv  =  12  samples  with  the  F2 
wavelet  filter  is  illustrated  in  Figures  5.1  and  5.2.  These  figures  show  segments 
of  original  and  F2-filtered  wall-pressure  signals  and  their  corresponding  V  functions 
for  the  equilibrium  and  disturbed  boundary  layers,  respectively.  The  event  clusters 
identified  by  the  LVAR  time  window  in  the  unfiltered  data  correspond  to  the  flow 
structures  retained  in  the  F2-filtered  data. 
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Figure  5.1:  Equilibrium  wall  pressure  and  localized  variance  with  T+  =  9.6.  (a),  (b) 
Prior  to  F2  filtering,  (c),  (d)  After  F2  filtering. 


Similar  results  were  obtained  for  the  uv  time  records.  Thus,  the  LVAR  technique 
appears  to  be  a  simple  alternative  to  the  wavelet  filter.  However,  locating  the  indi¬ 
vidual  clusters  would  require  querying  the  V  function,  thus  introducing  a  threshold 
level  similar  to  the  VITA  method.  Therefore,  V  functions  offer  no  advantage  over 
other  threshold  techniques  for  conditional  sampling.  Also,  since  the  LVAR  function 
is  positive-definite,  discrimination  of  positive  and  negative  events  must  be  done  by 
returning  to  the  original  time  records.  However,  this  property  means  that  correlation 
computations  will  be  enhanced  if  conducted  on  V  functions  instead  of  the  signals 
themselves.  The  use  of  V  functions  in  correlations  will  be  discussed  in  Section  5.3. 


5.2.2  Localized  windowed  peak-detection  technique 

A  new  algorithm,  called  Localized  Windowed  Peak  Detection  (LWPD),  was  devel¬ 
oped  that  captures  the  cluster  patterns  of  peak  events  while  being  independent  of 
a  magnitude  threshold.  The  method  is  based  on  the  conjecture  that  the  events  of 
interest  are  the  largest  peaks  in  the  event  clusters. 

The  algorithm  is  described  as  follows  for  a  sequence  s  =  {si . . .  s^}  sampled  at 
discrete  times  t  =  {t\ . . .  t#}: 

1.  Find  the  local  maxima,  P,  of  s  at  times  tn  =  {tni . . .  tnp}.  Define  the  event 
start/end  as  the  local  minima  immediately  preceding/following  each  peak. 
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Figure  5.2:  Disturbed  wall  pressure  and  localized  variance  with  T+  =  4.8.  (a),  (b) 
Prior  to  F2  filtering,  (c),  (d)  After  F2  filtering. 

2.  Sequentially  compare  each  Pn  with  its  neighbors  and  compute  the  times  be¬ 
tween  peaks  A t+  —  tni+1  —  tni  and  At_  =  tni  —  tn%_1 .  Retain  peaks  from  P  to 
form  a  new  subset  of  peaks,  P',  if  one  of  the  following  conditions  is  met: 

(a)  A  peak  is  itself  a  local  maximum  of  P.  Let  the  start  point  be  equal  to  the 
end  point  of  the  preceeding  (now  eliminated)  peak  and  let  the  end  point 
be  equal  to  the  start  point  of  the  following  (now  eliminated)  peak. 

(b)  A  peak  is  isolated  on  only  one  side  such  that  A t+  >  Tw  or  A >  Tw, 
and  it  is  a  local  maximum  within  the  window  Tw.  If  it  is  isolated  from  the 
next  peak  (A t+  >  Tw ),  let  the  start  point  be  equal  to  the  end  point  of  the 
preceeding  (now  eliminated)  peak.  If  it  is  isolated  from  the  preceeding 
peak  (A >  Tw),  let  the  end  point  be  equal  to  the  start  point  of  the 
following  (now  eliminated)  peak. 

(c)  A  peak  is  isolated  on  both  sides  such  that  At+  >  Tw  and  Af_  >  Tw.  The 
start/end  points  remain  unchanged. 

3.  If  P'  =  P,  i.e.  if  no  peaks  have  been  eliminated,  stop.  Otherwise,  let  P  —  P' 
and  iterate  from  step  2. 

A  complete  listing  of  the  LWPD  code  is  contained  in  appendix  C.  The  final  output 
contains  the  location  of  the  largest  peak  within  each  cluster  and  the  duration  of 
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Figure  5.3:  Localized  windowed  peak  detection  with  Tw  =  50  samples. 


each  individual  cluster.  The  time  window,  Tw,  is  the  only  adjustable  parameter 
of  the  algorithm.  It  specifies  the  minimum  time  between  successive  peaks  in  a 
sequence.  Without  this  parameter,  the  “peak  of  the  peaks”  operation,  step  2a,  would 
eventually  eliminate  all  points  except  the  maximum  peak  in  s.  This  algorithm  is 
more  computationally  intense  than  the  traditional  peak  level  detection  algorithm, 
but  surprisingly  robust  for  a  wide  range  of  random-like  time  records. 

The  LWPD  process  is  illustrated  in  Figure  5.3  for  a  segment  of  u-velocity  from 
the  equilibrium  TBL  at  y/S  —  0.02.  A  time  window  equal  to  50  samples  or  40t+  was 
used.  Circles  mark  the  peaks  retained  after  each  step.  The  step  functions  show  the 
evolution  of  the  cluster  indicator  functions  based  on  the  cluster  beginning  and  end 
points.  The  initial  pass  retains  only  “peaks  of  the  peaks,”  with  the  window  criteria 
coming  into  play  in  the  subsequent  passes.  Three  event  clusters  are  indicated  in  the 
segment  shown. 

The  final  cluster  indicator  functions  merit  further  attention.  The  durations  of 
the  first  two  clusters  appear  to  be  roughly  correct,  if  not  in  optimum  alignment. 
The  indicator  function  for  the  third  event,  while  centered  on  the  peak,  appears  too 
short.  It  is  concluded  from  this  limited  sampling  of  the  data  that  the  definition  of 
event  start/ end  points  needs  further  testing  and  possibly  some  refinements.  There¬ 
fore,  event  duration  results  from  the  LWPD  algorithm  should  be  viewed  with  some 
caution. 

These  results  can  be  contrasted  with  traditional  peak-event  detection  with  a 
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Table  5.1:  LWPD  statistics  for  positive  wall  pressure  peaks  with  Tw  =  75  samples . 


Case 

T+ 

w 

A  t+ 

AT+ 

Number  of 
clusters 

Frequency 
per  s 

Frequency 
per  1000t+ 

ETBL 

60 

47.3 

112.0 

2334 

234.1 

8.9 

ETBL,  F2 

60 

47.7 

112.4 

1866 

233.3 

8.9 

DTBL 

30 

22.5 

56.0 

2327 

233.4 

17.9 

DTBL,  F2 

30 

24.2 

57.1 

1831 

228.9 

17.5 

magnitude  threshold  k  =  2.  The  traditional  method  would  rigidly  exclude  the  peaks 
at  points  86  and  164,  which  have  magnitudes  of  approximately  1.5m,™,,.  Instead,  it 
would  identify  only  the  peak  at  sample  220.  In  addition,  if  there  were  a  zero  crossing 
between  the  peaks  at  samples  220  and  248,  both  would  be  identified  by  traditional 
peak-event  detection.  Ensemble  averaging  over  a  time  interval  greater  than  about  20 
samples  would  then  include  these  peaks  twice  in  different  relative  temporal  locations. 
By  aligning  ensemble  averages  on  the  largest  peak  in  each  event  cluster,  this  type  of 
duplication  can  be  minimized. 

The  single  parameter  Tw  determines  statistical  properties  such  as  number  of 
events,  average  duration,  and  average  time  between  cluster  patterns.  It  is  important 
to  note  that  the  number  of  cluster  patterns  detected  by  the  LWPD  technique  is 
mainly  governed  by  Tw,  with  weak  dependence  on  the  distinct  features  of  the  time 
record.  This  can  be  understood  if  one  considers  a  signal  that  has  large  amplitude 
clusters  separated  by  periods  of  relatively  low  amplitude  activity.  If  Tw  is  less  than 
the  time  between  large  amplitude  clusters,  the  LWPD  algorithm  will  still  locate  the 
largest  peak  in  this  interval.  It  is  just  that  such  a  cluster  will  have  a  smaller  impact  on 
correlations  or  ensemble  averages  based  on  the  detected  clusters.  This  characteristic 
of  the  LWPD  technique  was  exploited  by  defining  the  same  time  window  value  Tw 
for  both  uv  and  p  time  records.  This  insured  conformity  in  the  number  of  events  for 
both  records.  The  particular  choice  of  Tw  was  based  on  maximizing  the  correlation 
between  concurrent  events,  which  is  discussed  in  Section  5.3.  The  bounds  for  this 
number,  established  from  extensive  computations  on  the  equilibrium  time  records, 
are  40  <  T+  <  120  (50  <  Tw  <  150  data  samples). 

The  value  of  Tw  =  75  samples,  was  chosen  for  the  equilibrium  data  (T+  =  60) 
based  on  numerous  test  cases  and  applied  also  to  the  disturbed  flow  data  (T+  =  30) . 
The  resulting  cluster  statistics  for  the  positive  wall  pressure  peaks  are  displayed  in 
Table  5.1.  Using  the  same  physical  value  for  Tw  for  both  flows  results  in  equivalent 
statistics  in  physical  terms,  as  evidenced  by  the  number  and  frequency  of  physical 
clusters  detected.  Recall  that  the  F2-filtered  time  records  are  8-s  long  compared 
with  10  s  for  the  unfiltered  data.  In  nondimensional  terms,  the  statistics  directly 
reflect  the  factor  of  two  difference  in  viscous  time  scales.  Comparing  with  Table  3.4, 
the  equilibrium  LWPD  statistics  are  equivalent  to  peak-event  detection  results  for 
a  threshold  between  1  and  2,  while  the  disturbed  flow  statistics  are  equivalent  to 
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Figure  5.4:  LWPD  cluster  detection  on  equilibrium  pF2  and  uvF2  with  T+  =  60. 
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a  threshold  equal  to  1.  In  both  cases,  the  LWPD  durations  are  longer,  reflecting 
clusters  rather  than  individual  events. 

Using  this  window  size,  results  for  F2-filtered  equilibrium  and  disturbed  wall 
pressure  and  Reynolds  stress  are  shown  in  Figures  5.4  and  5.5,  respectively.  The 
indicator  functions  marking  the  cluster  beginning  and  end  points  are  shown  shifted 
upwards  for  clarity.  The  clusters  were  obtained  by  applying  the  LWPD  algorithm 
to  the  positive  peaks  for  the  pressure  and  to  the  negative  peaks  for  the  Reynolds 
stress.  The  time  segments  are  the  same  as  those  in  Figures  5.1  and  5.2  illustrating 
the  LVAR  technique.  Comparing  with  these  figures,  one  sees  that  there  is  good 
correspondence  between  LWPD  clusters  and  increased  local  variance.  Cluster  time 
records  will  be  denoted  by  a  subscript  cl,  and  are  defined  as  the  (F2-filtered)  time 
records  multiplied  by  the  LWPD  indicator  function: 

fd{t)  =  fFs(t )  •  ^lwpd(^)  (5.3) 


where 


^lwpd(^)  = 


1,  if  cluster  start  <  t  <  cluster  end 
0,  otherwise 


(5.4) 


5.3  Correlation  of  clusters 

Correlations  were  performed  on  three  time  records:  original  data,  V  functions,  and 
LWPD  clusters  [5.8].  These  correlations  were  computed  between  the  two  spatially 
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Figure  5.5:  LWPD  cluster  detection  on  disturbed  pF2  and  uvF2  with  T+  =  30. 


separated  wall  transducers,  pi  and  P2,  between  the  Reynolds  stress  and  the  wall 
pressure  at  the  same  streamwise  location  uv  and  px,  and  between  the  Reynolds  stress 
and  the  upstream  wall  pressure  uv  and  p2.  The  uv-p  correlations  were  computed  for 
all  y-locations. 

5.3.1  Estimation  of  convection  velocity 

The  experimental  arrangement  included  two  wall-pressure  transducers  separated  by 
a  streamwise  distance  Ax  =  1.27  cm.  This  allowed  the  convection  velocity  of  the 
organized  motions  to  be  computed  from 

Uc  =  AxJt j max  R  (5.5) 

where  r|max  R  is  the  time  delay,  which  maximizes  the  cross  correlation  function 
between  px  and  p2,  RPl,P2(T)-  The  correlation  was  formed  by  dividing  the  total 
time  records  into  windows  of  equal  size,  performing  the  correlation  between  pairs  of 
windows,  and  then  ensemble  averaging.  A  window  size  of  2048  points  (62.5  ms)  was 
used.  Similar  calculations  were  made  using  the  localized  variance  records  VPl  ( t )  and 
VP2(t)  and  the  LWPD  cluster  records  P\cl{t)  and  P2c,(t)-  For  completeness,  both  the 
unfiltered  and  F2-filtered  cases  were  examined. 

Special  attention  must  be  given  to  computing  correlations  between  the  localized 
variance  records,  since  these  functions  are  positive-definite.  Such  functions,  even 
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Table  5.2:  Convection  velocity  of  organized  structures  for  equilibrium  data 


time 

record 

u 7 

Uc/ue 

Unfiltered  Filtered 

Unfiltered  Filtered 

Original 

17.4  16.5 

(37)  (39) 

0.70  0.66 

LVAR 

(clusters) 

16.5  16.9 

(39)  (38) 

0.66  0.68 

LWPD 

(clusters) 

16.5  16.9 

(39)  (38) 

0.66  0.68 

(  )  time  delay  in  sampled  points. 


U+  =  16.9  ±  0.5 
Uc/Ue  =  0.68  ±  0.02 


Table  5.3:  Convection  velocity  of  organized  structures  for  disturbed  data. 


time 

record 

u: 

uc/ue 

Unfiltered  Filtered 

Unfiltered  Filtered 

Original 

18.5  19.8 

(50)  (47) 

0.52  0.55 

LVAR 

(clusters) 

20.2  24.4 

(46)  (38) 

0.56  0.68 

LWPD 

(clusters) 

20.6  22.1 

(39)  (42) 

0.58  0.62 

(  )  time  delay  in  sampled  points. 

U+  =  22.1  ±  2.0 
Uc/Ue  =  0.62  ±  0.06 


uncorrelated  random  noise,  will  produce  correlation  functions  with  a  characteristic 
sawtooth  shape,  the  peak  of  which  is  at  zero  time  delay.  Therefore,  when  performing 
correlations  on  localized  variance  records,  the  results  are  reported  in  dB  referenced 
to  the  base  sawtooth  shape.  This  will  be  referred  to  as  a  normalized  correlation, 
NC. 

The  results  for  the  convection  velocity  computed  for  the  equilibrium  and  dis¬ 
turbed  flows  are  shown  in  Tables  5.2  and  5.3,  respectively.  The  equilibrium  cal¬ 
culations  obtained  from  the  three  distinct  time  records  are  surprisingly  consistent. 
The  nominal  result  Uc/Ue  =  0.68  compares  with  a  value  of  0.67  obtained  by  Thomas 
and  Bull  from  similar  measurements  in  an  equilibrium  TBL  [5.9].  For  the  disturbed 
flow,  the  results  cover  a  wider  range,  particularly  between  the  unfiltered  and  filtered 
cases.  However,  the  filtered  LVAR  and  LWPD  time  records  provided  results  that 
were  most  consistent  with  the  equilibrium  results,  indicating  the  validity  of  these 
detection  schemes. 
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Time  delay  (s)  x  1  q-3 

Figure  5.6:  Normalized  correlation  between  VPl  and  Vuv  across  the  equilibrium  TBL. 
The  different  ^-locations  are  shifted  upwards  for  clarity. 

5.3.2  Correlations  between  turbulence  and  wall  pressure 

To  illuminate  the  shape  and  extent  of  the  coherent  motions,  correlations  were  com¬ 
puted  between  the  wall  pressure  and  the  Reynolds  stress  across  the  boundary  layers. 
These  calculations  were  performed  using  the  original,  V,  and  LWPD  cluster  time 
records.  Only  the  results  on  V  functions  are  presented.  Similar  results  were  ob¬ 
tained  for  the  original  and  LWPD  time  records. 

The  normalized  correlations  between  the  wall  pressure  and  the  Reynolds  stress 
were  made  at  each  y-location  across  the  boundary  layers.  The  results  are  shown  in 
Figure  5.6  for  the  equilibrium  flow  and  Figure  5.7  for  the  disturbed  flow.  These 
results  are  for  transducer  pi  at  the  same  streamwise  location  as  the  X-wire.  For 
clarity,  each  y-location  is  shown  shifted  upwards  by  1  dB  relative  to  the  previous 
station.  The  wall  pressure  and  Reynolds  stress  are  correlated  out  to  a  distance  y/5  « 
0.2.  Beyond  this  point,  significant  peaks  in  the  correlations  cannot  be  unequivocally 
defined.  Notice  that,  for  both  flows,  the  time  delay  for  maximum  correlation  is  zero 
near  the  wall,  shifting  to  increasingly  negative  values  across  the  TBL.  A  quantitative 
comparison  of  the  correlations  between  the  equilibrium  and  disturbed  flows  can  be 
made  from  Figure  5.8,  which  plots  results  at  two  specific  locations.  For  both  flows, 
the  negative  time  delay  for  maximum  correlation  at  y/5  =  0.14  indicates  that  the 
organized  motion  burst  mechanism  extends  to  the  log-law  region  and  has  a  forward 
angle  of  inclination. 
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Time  delay  (s)  x  1 Q-3 

Figure  5.7:  Normalized  correlation  between  VPl  and  VUXJ  across  the  disturbed  TBL. 
The  different  y-locations  are  shifted  upwards  for  clarity. 

The  time  delays  for  maximum  correlation  versus  position  in  the  boundary  layer 
are  presented  in  Figure  5.9.  Results  are  also  included  for  the  correlations  of  Reynolds 
stress  with  the  upstream  pressure  transducer  p2.  These  data  are  more  scattered  than 
for  pi  because  of  the  larger  distance  between  transducer  p2  and  the  X-wire  probe. 
Some  points  from  both  plots  have  been  omitted  due  to  the  previously  mentioned  dif¬ 
ficulties  in  identifying  a  dominant  peak  outside  of  y/5  =  0.2.  Even  so,  the  trends  are 
well-defined,  and  indicate  roughly  the  same  shape  for  both  flows.  A  representative 
inclination  angle  can  be  inferred  from  the  data  in  Figure  5.9(b)  by  noting  the  y- 
location  where  the  time  delay  is  zero.  For  the  equilibrium  flow,  this  point  occurs  for 
a  y/S  between  0.12  and  0.14,  which  equates  to  an  angle  between  15.5  and  18.3  deg. 
This  is  consistent  with  the  results  of  Snarski  [5.10]  and  Thomas  and  Bull  [5.9],  who 
found  the  dominant  structures  to  be  inclined  at  18  deg.  For  the  disturbed  flow,  the 
point  of  zero  time  delay  is  at  y/5  «  0.17,  which  roughly  equates  to  an  inclination 
angle  of  30.4  deg. 

5.4  Conditional  sampling  on  event  clusters 

The  features  of  the  cluster  patterns  were  examined  from  conditional  sampling  based 
on  the  locations  obtained  from  the  LWPD  algorithm.  The  ensemble  averages  were 
centered  around  the  maximum  peak  in  each  detected  cluster,  using  the  original  (F2- 
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Figure  5.8:  Normalized  correlation  between  VPl  and  Vuv  at  y/6  =  0.02  and  0.14  in 
the  equilibrium  and  disturbed  flows. 

filtered)  time  records.  As  in  Section  4.5.2,  the  detection  variable  is  plotted  using  a 
thicker  solid  line.  In  the  present  section,  results  at  one  y-location  will  be  compared  to 
those  obtained  in  Chapter  4  using  the  traditional  event-detection  methods  based  on 
an  amplitude  threshold  level.  Results  spanning  the  boundary  layers  will  be  discussed 
in  Chapter  6. 

5.4.1  Quadrant  cluster  detection 

The  LWPD  algorithm  was  applied  to  the  negative  of  the  Q2  Reynolds  stress  time 
records,  uv2(t),  to  determine  the  cluster  locations  based  on  negative  uv2{t)  peaks. 
Figures  5.10  and  5.11  display  the  conditional  sampling  results  based  on  these  loca¬ 
tions  for  the  equilibrium  and  disturbed  flows,  respectively.  The  results  for  both  flows 
are  almost  identical  with  averages  based  on  peak  quadrant  event-detection,  shown 
previously  in  Figures  4.20  and  4.22.  The  primary  difference  is  that  the  present  re¬ 
sults  are  of  a  reduced  magnitude.  The  LWPD  algorithm  is  able  to  identify  significant 
flow  structure  without  specifying  a  magnitude  threshold.  However,  the  Q4-Q2-Q4 
event  pattern  was  not  enhanced  by  conditioning  on  clusters. 
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Figure  5.9:  Time  delays  for  maximum  correlation  between  wall  pressure  and 
Reynolds  stress  from  Figures  5. 6-5.8.  (a)  pi  vs.  uv.  (b)  p2  vs.  uv. 

5.4.2  Peak  wall-pressure  cluster  detection 

Since  one  objective  of  the  research  is  to  confirm  a  bidirectional  relationship  between 
wall-pressure  and  turbulent  structures,  conditional  sampling  was  performed  trig¬ 
gered  on  positive  wall  pressure  clusters.  The  results  are  displayed  in  Figures  5.12 
and  5.13  for  the  equilibrium  and  disturbed  flows,  respectively.  Again,  the  results 
are  almost  identical  to  the  conditional  sampling  based  on  traditional  peak-pressure 
event  detection  shown  in  Figures  4.16  and  4.18.  Unlike  the  conditional  sampling 
based  on  Q2  clusters,  triggering  on  LWPD  wall  pressure  clusters  has  not  sacrificed 
the  strengths  of  the  correlations.  The  most  important  result  is  that  conditional  sam¬ 
pling  on  pressure  indicates  that  the  Q2  events  are  aligned  with  the  positive  pressure 
peaks.  This  is  most  evident  in  the  equilibrium  case.  The  uv2  peak  in  the  disturbed 
case  is  less  pronounced.  The  degradation  is  probably  due  to  an  increase,  over  the 
equilibrium  case,  in  the  variation  of  convection  velocities,  or  phase,  between  individ¬ 
ual  realizations.  This  effect  is  known  as  phase  jitter  [5.11].  Even  so,  the  following  uvA 
excursion  is  equally  significant  in  both  flows.  One  difference  between  the  equilibrium 
and  disturbed  flows  is  the  relative  strengths  of  the  u  and  v  motions.  The  equilibrium 
flow  shows  a  strong  positive  u  acceleration  with  a  relatively  weaker  v  deceleration. 
In  contrast,  the  v  motion  is  much  stronger  in  the  disturbed  flow  case,  consistent  with 
the  increased  vertical  motions  presented  in  Chapter  3.  By  comparing  the  results  in 
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Figure  5.10:  Ensemble  averages  conditioned  on  equilibrium  flow  Q2  LWPD  clusters 
at  y/5  —  0.02. 


Figures  5.10  through  5.13,  it  appears  that  a  bidirectional  relationship  exists  between 
positive  pressure  peaks  and  Q2  events  for  both  equilibrium  and  disturbed  flows. 

5.5  Summary 

Two  new  methods  have  been  developed  to  identify  and  study  clusters  of  events: 
LVAR  and  LWPD.  The  positive-definite  nature  of  the  localized  variance  function 
makes  it  ideal  for  correlation  applications,  while  the  LWPD  algorithm,  with  its 
precise  locations  of  clusters  in  time,  is  more  suited  to  conditional  sampling.  These 
methods  differ  from  traditional  event-detection  methods  in  that  they  do  not  rely 
on  the  selection  of  an  amplitude  threshold  level.  Instead,  both  methods  require 
the  specification  of  a  time  parameter.  For  the  LVAR  technique,  this  parameter 
corresponds  to  a  characteristic  time  scale,  similar  to  a  wavelet  filter  at  fixed  scale. 
For  LWPD  detection,  this  parameter  specifies  the  minimum  time  between  successive 
clusters.  The  algorithm  insures  that  the  most  significant  clusters  are  identified. 

Correlations  using  V  functions  were  successful  in  deducing  the  convection  velocity 
of  the  wall-pressure  signature  of  the  organized  motions  and  the  spatial  extent  of  the 
organized  motions.  The  results  were  consistent  with  the  findings  obtained  by  other 
researchers. 
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Figure  5.11:  Ensemble  averages  conditioned  on  disturbed  flow  Q2  LWPD  clusters  at 
yfS  =  0.02. 


Conditional  sampling  based  on  wall  pressure  and  Q2  LWPD  clusters  reproduced 
the  results  based  on  traditional  peak  and  quadrant  event-detection  techniques  pre¬ 
sented  in  Chapter  4.  The  results  conditioned  on  Q2  clusters  showed  an  alignment 
of  the  positive  pressure  peak  with  the  Q2  event,  with  Q4  events  on  either  side.  The 
results  conditioned  on  positive  wall  pressure  clusters  showed  a  correlation  with  Q2- 
Q4  event  pairs.  The  bidirectional  relationship  between  positive  wall-pressure  peaks 
and  accelerating  u-velocity  or  second-quadrant  uv  motions  was  confirmed.  While 
the  results  may  not  be  a  significant  improvement  over  traditional  methods,  in  the 
sense  of  increased  correlation,  the  fact  that  LWPD  detection  produced  equivalent 
results  suggests  that  the  idea  of  event  clusters  is  valid.  The  problem  of  phase  jitter 
between  the  conditionally  sampled  events  remains  to  be  addressed. 
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Figure  5.12:  Ensemble  averages  conditioned  on  equilibrium  flow  wall  pressure  LWPD 
clusters  at  y/8  =  0.02. 


Figure  5.13:  Ensemble  averages  conditioned  on  disturbed  flow  wall  pressure  LWPD 
clusters  at  y/5  =  0.02. 
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Chapter  6 

Interpretations  and  conclusions 


The  results  shown  in  Chapters  4  and  5  confirm  that  wall-pressure  peak  events  are 
footprints  of  near-wall  motions  for  both  equilibrium  and  disturbed  flows.  Specifi¬ 
cally,  conditional  sampling  results  near  the  wall  revealed  a  Q2-Q4,  or  ejection-sweep 
sequence  bracketing  the  positive  pressure  peak.  From  correlation  results  across  the 
boundary  layers,  it  was  inferred  that  the  structures  are  inclined  downstream.  This 
establishes  the  framework  to  address  the  two  central  questions  of  the  research: 

1.  What  are  the  physical  features  of  these  near- wall  flow  structures? 

2.  Are  these  near- wall  structures  coupled  to  and  affected  by  the  outer-layer  dy¬ 
namics? 

The  first  question  will  be  addressed  by  examining  the  ensemble-averaged  flow  field 
conditioned  on  wall-pressure  cluster  events  for  the  equilibrium  TBL.  The  second 
question  will  be  discussed  by  comparing  the  equilibrium  and  disturbed  flows  for 
common  near- wall  structures.  Implicit  to  both  questions  is  the  issue  of  filtering, 
which  was  necessary  for  the  extraction  of  the  near- wall  motions. 

6.1  Physical  features  of  near- wall  flow  structures 

To  explore  the  spatial  extent  of  the  observed  organized  motions,  the  ensemble  average 
of  u  and  v,  conditioned  on  positive  wall-pressure  LWPD  clusters,  was  constructed  at 
each  y  measurement  station.  If  one  assumes  that  the  structures  remain  intact  over  a 
short  distance  and  that  they  convect  with  the  local  mean  streamwise  velocity,  then 
the  time  domain  can  be  mapped  into  space: 

xu,v  —  —U(y)  •  t  (6.1) 

Likewise,  the  wall-pressure  distribution  can  be  inferred  from  the  ensemble-averaged 
time  history  using  the  convection  velocities  deduced  in  Section  5.3.1: 

xp  =  -Uc  ■  t  (6.2) 
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These  assumptions  allow  the  construction  of  a  composite  picture  of  the  near-wall 
flowfield  structure  associated  with  positive  wall-pressure  peaks.  This  approach  was 
validated  by  Kim  [6.1],  who  computed  ensemble  averages  based  on  the  VITA  tech¬ 
nique  using  data  obtained  from  a  large-eddy  simulation  of  turbulent  channel  flow. 
The  results  are  presented  in  Figures  6.1  and  6.2  for  the  equilibrium  flow  before  and 
after  F2  filtering,  respectively. 

In  these  figures,  the  wall- pressure  distributions  associated  with  each  y  location  are 
plotted  to  emphasize  the  consistency  of  the  trigger  condition.  The  flow  field  displays 
organized  motions  that  extend  out  to  a  distance  y/5  =  0.2,  beyond  which  correlated 
structure  is  not  evident.  For  clarity,  the  second  and  fourth  y  measurement  stations 
have  been  omitted,  and  the  magnitude  of  the  velocity  vectors  has  been  arbitrarily 
scaled.  The  scaling  factor  for  the  post-filtered  case  (Figure  6.2)  was  twice  that  of 
the  pre-filter  case,  reflecting  the  lower  energy  content. 

The  most  important  feature,  visible  in  both  the  unfiltered  and  the  filtered  data, 
is  that  the  positive  pressure  peaks  are  aligned  with  ejections  followed  by  strong 
sweeps  from  upstream.  This  finding  was  shown  earlier  in  Chapter  5  in  Figure  5.12. 
Moreover,  conditional  sampling  based  on  Q2  clusters  confirmed  the  bidirectionality 
of  this  result  (cf.  Figure  5.10). 

First,  consider  the  case  before  wavelet  filtering  (Figure  6.1).  The  Q2-Q4  interface 
displayed  in  the  composite  picture  forms  a  shear  layer  that  is  inclined  downstream 
with  an  inclination  of  approximately  45  deg.  Overall,  the  flow  field  is  consistent  with 
the  schematic  diagram  of  the  near-wall  motion  proposed  by  Astolfi  and  Forestier 
(Figure  1.3).  The  extent  of  the  wall-pressure  footprint,  approximately  160 v/uT  or 
0.145,  also  agrees  with  their  finding  of  110  to  130 u/uT.  Note  that  the  pressure  minima 
on  either  side  of  the  peak  do  not  appear  in  the  unfiltered  data  but  do  appear  after 
filtering. 

Next,  consider  the  effects  of  filtering  (Figure  6.2).  Application  of  the  F2  filter 
changes  the  flowfield  picture  in  two  ways.  First,  the  ejection/sweep  motions  are 
steepened,  and  the  dividing  shear  front  is  sharpened.  Second,  significant  ancillary 
structures  are  revealed  in  the  form  of  three  vortical  motions  of  alternating  rotation. 

These  structures  bracket  the  shear  front  and  merit  discussion.  Proceeding  down¬ 
stream,  the  first  structure  appears  as  ejection  motions  that  arch  around  downstream 
to  be  entrained  in  a  sweep  motion.  This  vortex  motion  is  not  confined  to  one  partic¬ 
ular  vertical  location,  but  spans  the  region  y/5  <  0.2.  The  second  structure  consists 
of  a  vortex  core  located  above  the  upstream  pressure  minimum  at  y/5  —  0.15.  This 
is  the  strongest  of  the  ancillary  structures  and  appears  to  be  driving  the  main  sweep 
motion.  The  third  structure  is  a  vortex  core  just  downstream  of  the  pressure  peak 
at  y/5  =  0.10.  This  vortex  appears  to  be  inducing  the  ejections  near  the  wall. 

The  flow  field  pictured  in  Figure  6.2  is  consistent  with  the  models  of  turbulent 
boundary-layer  structure  discussed  in  Chapter  1.  Specifically,  it  is  suggested  that  the 
two  structures  with  negative  vorticity  are  the  heads  of  inclined  horseshoe  vortices. 
However,  the  strong  vortex  in  between  is  puzzling,  since  its  sign  is  opposite  to  that 
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Figure  6.1:  Ensemble-averaged  equilibrium  perturbation  flow  field  conditioned  on  wall-pressure  LWPD  clusters  prior  to 
F2  filtering  and  in  a  reference  frame  moving  with  the  mean  streamwise  velocity. 
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Figure  6.2:  Ensemble-averaged  equilibrium  perturbation  flow  field  conditioned  on  wall  pressure  LWPD  clusters  after  F2 
filtering  and  in  a  reference  frame  moving  with  the  mean  streamwise  velocity. 
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Figure  6.3:  The  mechanics  of  streak  breakup  after  Kline,  et  al.  1967  [6.2]. 


of  a  head  or  neck.  More  likely,  the  sweep  and  inclined  shear  layer  appear  so  strong 
and  highly  correlated  precisely  because  of  their  relative  alignment  with  the  ensemble 
trigger  event,  namely,  the  positive  pressure  peak.  In  other  words,  the  fundamental 
flow  structure  consists  of  two  horseshoe  vortices  of  negative  vorticity,  one  ahead  and 
above  the  other,  that  induce  the  ejection/sweep  pattern  and  a  concomitant  wall- 
pressure  peak.  This  situation  was  proposed  by  Kline  et  al.,  who  put  forth  the  model 
of  streak  breakup  shown  in  Figure  6.3.  The  two  vortex  shapes  and  their  induced 
flowfield  agree  with  Figure  6.2.  However,  the  breakup  of  the  structure  because  of  a 
secondary  instability  has  been  superseded  (cf.  Section  1.2.1). 


6.2  Effects  of  outer-layer  disturbances 

To  assess  the  universality  of  the  near-wall  motions  and  their  wall-pressure  signatures 
observed  for  the  equilibrium  flow,  a  disturbed  TBL  was  studied.  The  data  presented 
in  Chapter  3  showed  that  the  disturbed  TBL  differed  from  the  equilibrium  TBL 
by  the  addition  of  large-scale,  low-frequency  disturbances  in  the  outer  regions.  The 
inner  layer  appeared  to  have  been  established  quickly  downstream  after  reattachment 
while  the  outer  layer  decayed  more  slowly.  The  hypothesis  was  made  that  the  inner- 
and  outer-layer  portions  of  turbulent  boundary  layers  are  uncoupled. 

The  composite  flow  pictures  for  the  disturbed  flow  case  are  shown  in  Figures  6.4 
and  6.5  for  the  unfiltered  and  filtered  data,  respectively.  Because  of  the  different 
time  scales  of  the  flows,  twice  as  many  data  points  are  available  for  the  disturbed 
flow.  For  clarity  of  comparison,  only  every  other  data  point  has  been  plotted.  The 
scaling  factors  for  the  vector  magnitudes  are  identical  to  those  used  in  Figures  6.1 


113 


250 


NSWCDD/TR-95/222 


BJ|ap/A 


suiid/d 


114 


Figure  6.4:  Ensemble-averaged  disturbed  perturbation  flow  field  conditioned  on  wall-pressure  LWPD  clusters  prior  to 
F2  filtering  and  in  a  reference  frame  moving  with  the  mean  streamwise  velocity. 
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and  6.2,  respectively.  This  allows  a  direct  comparison  of  the  structures  between  the 
equilibrium  and  disturbed  flows. 

As  in  the  previous  section,  consider  first  the  case  before  wavelet  filtering.  Fig¬ 
ure  6.4  shows  significant  differences  from  the  equilibrium  case  (Figure  6.1).  The 
most  significant  difference  is  that  a  shear  front  is  not  apparent.  A  large-scale,  grad¬ 
ual  turning  of  the  flow  dominates  even  near  the  wall.  The  outward  motion  is  spread 
out  downstream  of  the  wall-pressure  peak.  At  the  wall,  the  pressure  signature  is 
also  extended  with  a  less  pronounced  peak.  The  near-wall  motions  are  clearly  con¬ 
taminated  by  the  outer-layer  disturbances. 

Application  of  the  wavelet  filter  in  Figure  6.5  shows  a  more  complex  picture 
than  for  the  equilibrium  flow,  but  significant  structure  is  again  revealed.  The  ejec¬ 
tion/sweep  interface  is  sharpened  and  is  aligned  with  the  pressure  peak.  Further¬ 
more,  the  same  sequence  of  vortical  motions  revealed  in  the  equilibrium  flow  has 
become  evident.  Arching  ejection  motions  span  the  vertical  range  and  feed  into  the 
primary  sweep  motion.  However,  in  this  case,  the  motions  near  the  wall  form  a 
closed  vortex  over  the  upstream  pressure  minimum  at  y/8  m  0.03.  A  vortex  of  posi¬ 
tive  vorticity  is  again  seen  above  the  primary  sweep  motion,  but  the  negative  vortex 
downstream  of  the  shear  front  is  not  as  highly  correlated  as  in  the  equilibrium  case. 
Overall,  the  ejection/sweep  motion  is  more  compact  than  for  the  equilibrium  flow, 
spanning  approximately  0.15,  but  the  angle  of  inclination  is  the  same. 

The  observations  made  for  the  filtered  data  support  the  hypothesis  that  the  near¬ 
wall  motions  associated  with  wall-pressure  peak  events  are  largely  unaffected  by  the 
outer-layer  disturbances. 


6.3  Accomplishments 

A  summary  is  presented  of  the  major  achievements  of  the  research  investigation. 

•  High-resolution,  high-quality  turbulence  data  for  two  flow  configurations  have 
been  obtained.  These  archived  databases  are  available  to  the  research  com¬ 
munity  for  validation  of  computational  and  analytical  modeling  of  turbulent 
flows. 

•  Extensive  software  was  developed  for  both  wavelet  filtering  and  signal  process¬ 
ing  of  the  raw  data.  The  documented  software  was  written  in  modular  form 
(MATLAB  and  C)  and  has  been  made  available  to  the  naval  laboratories  and 
scientific  community. 

•  The  research  questions  were  addressed  by  performing  analyses  of  the  raw  data 
for  both  the  equilibrium  and  disturbed  flows.  The  following  deductions  were 
made  on  the  causal  relationship  between  near-wall  organized  motions  and  wall- 
pressure  events: 
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-  The  results  strongly  confirm  that  patterns  of  large-amplitude  wall-pressure 
events  are  footprints  of  the  active  (turbulence-producing)  motions.  These 
pressure  signatures  contain  both  positive  and  negative  peak  events. 

-  Wavelet  filtering  techniques  appear  to  be  successful  in  extracting  the  wall- 
pressure  signatures  of  the  near-wall  motions.  This  was  most  apparent  in 
the  case  of  the  disturbed  flow. 

-  Clusters  of  burst  events  and  wall-pressure  events  are  highly  correlated. 
Distinct  patterns  of  burst  events  (pairing  of  Q2  and  Q4  events)  are  ob¬ 
served  from  conditional  sampling  based  on  locations  of  clusters  of  wall- 
pressure  events.  The  results  show  two  vortical  structures,  one  ahead 
and  above  the  other,  that  induce  the  ejection/sweep  pattern.  These  mo¬ 
tions  are  detected  by  the  alignment  of  Q2  events  with  positive  peak  wall- 
pressure  events. 

-  The  near-wall  motions  appear  to  be  weakly  coupled  to  the  outer-layer 
dynamics.  This  was  tested  by  comparing  the  results  between  equilibrium 
and  disturbed  flows.  The  same  generic  features  of  the  wall  structures 
were  observed  for  both  flows. 

-  The  organized  motions  observed  by  wall-pressure  detection  are  inclined 
downstream  and  extend  to  the  log-law  region. 

-  The  physical  features  of  the  near- wall  motions  displayed  in  flowfield  map¬ 
pings  are  consistent  with  the  models  proposed  by  several  investigators. 


6.4  Future  work 

Many  of  the  findings  reported  in  this  investigation  require  further  validation.  In  ad¬ 
dition,  improvements  in  both  the  equipment  and  spatial  extent  of  the  measurements 
are  suggested.  Therefore,  recommendations  for  future  work  are  made  in  three  areas: 
testing,  measurements,  and  signal  processing. 


Testing 

Wall-pressure  data  indicate  that  the  traverse  system  for  the  hot  wire  causes  some 
contamination  in  the  midfrequency  range  due  to  vortex  shedding  from  the  mounting 
struts.  Filtering  techniques  were  used  to  minimize  this  noise,  but  it  is  suggested 
that  a  new  traverse  system  be  installed  to  eliminate  this  problem. 

Measurements 

It  is  strongly  recommended  that  hot-wire  measurements  be  made  at  several  stream- 
wise  and  spanwise  locations  across  the  boundary  layer.  Correlation  and  conditional 
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sampling  applied  to  these  data  will  provide  a  more  detailed  composite  picture  of 
the  flowfield  and  validate  the  proposed  physical  models.  In  addition,  it  is  proposed 
that  an  array  of  pressure  transducers  be  mounted  along  the  surface.  This  will  allow 
detailed  information  on  the  time  scales,  decay  rates,  and  convection  properties  of 
the  organized  motions.  This  type  of  information  is  required  in  the  design  of  adaptive 
techniques  for  turbulent  control. 


Signal  processing 

A  major  problem  in  the  analysis  of  turbulence  data  by  conditional  sampling  methods 
is  the  problem  of  phase  jitter  between  the  bursts  and  wall-pressure  events  that  are 
being  correlated.  When  the  data  are  ensemble  averaged,  significant  distortion  of  the 
true  correlation  can  appear.  Most  investigators  incorporated  a  phase-shift  correction 
to  each  event  based  on  the  localized  correlation.  Since  the  present  data  had  high 
resolution,  this  correction  was  not  needed.  However,  the  improvements  in  the  quality 
of  the  conditionally  sampled  results  presented  in  Chapters  4  and  5  will  be  apparent. 
An  algorithm  should  be  developed  that  implements  phase-jitter  corrections. 
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Appendix  A 

Hot-wire  and  wall-pressure  surveys 


This  appendix  documents  the  laboratory  conditions  and  the  primary  anemometer 
measurements  in  engineering  units.  Tabulations  are  given  for  the  four  surveys  from 
which  results  have  been  presented,  namely,  one  single-  and  one  X-wire  survey  for 
the  equilibrium  and  disturbed  TBL  flows.  Also,  results  are  presented  for  a  series  of 
wall-pressure  measurements  at  various  locations  downstream  of  the  disturbed  flow 
ramp.  This  was  done  to  ascertain  the  location  of  the  point  of  flow  reattachment. 
Figure  A.l  shows  the  wall-pressure  power  spectra  for  these  measurements.  The 
maximum  low-frequency  pressure  levels  occur  at  a  location  x/h  —  5.4,  while  the 
maximum  high-frequency  levels  occur  at  x/h  =  7.8.  Integrating  these  spectra  from 
50  Hz  to  the  noise  floor  at  12500  Hz  gives  the  rms  distribution  shown  in  Figure  A. 2. 
The  maximum  rms  levels  are  seen  to  be  at  x/h  ss  7.  From  these  observations,  it 
is  seen  that  the  primary  measurement  location  of  transducer  px  at  x/h  =  8.6  is  at 
least  one  step  height  downstream  of  reattachment. 
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Equilibrium  TBL  single-wire  survey 

Test  date:  06-Aug-94 

Reduction  date:  13-Aug-94 

Probe:  TSI  1261A-T1.5,  #914074 

Atmospheric  pressure:  101.6  KPa  (30.01  inHg) 

Static  temperature:  295.7  K  (532.2  R) 

Freestream  velocity,  Ue:  14.79  m/s  (  48.52  ft/s) 

Shear  velocity,  uT :  0.60  m/s  (  1.95  ft/s) 

Boundary-layer  thickness,  6:  2.72  cm  (  1.072  in) 


Survey 

station 

V 

(cm) 

y/s 

y+ 

U 

(m/s) 

u/ue 

U+ 

'U'rms 

(m/s) 

r^rms/ rr 

01 

0.01 

0.003 

3.5 

2.33 

0.158 

3.92 

0.70 

1.169 

02 

0.01 

0.004 

4.5 

3.36 

0.227 

5.63 

1.02 

1.720 

03 

0.01 

0.005 

5.5 

3.77 

0.255 

6.33 

1.13 

1.900 

04 

0.02 

0.007 

7.9 

4.48 

0.303 

7.53 

1.27 

2.132 

05 

0.02 

0.009 

9.4 

5.01 

0.339 

8.41 

1.35 

2.272 

06 

0.03 

0.011 

11.4 

5.32 

0.360 

8.93 

1.38 

2.317 

07 

0.04 

0.014 

15.4 

6.12 

0.414 

10.28 

1.43 

2.402 

08 

0.05 

0.018 

18.9 

6.58 

0.445 

11.05 

1.44 

2.417 

09 

0.06 

0.022 

23.3 

7.12 

0.481 

11.95 

1.42 

2.391 

10 

0.08 

0.029 

30.8 

7.69 

0.520 

12.90 

1.38 

2.325 

11 

0.11 

0.039 

41.2 

8.26 

0.559 

13.87 

1.31 

2.203 

12 

0.12 

0.045 

48.2 

8.55 

0.578 

14.36 

1.28 

2.144 

13 

0.19 

0.069 

73.5 

9.16 

0.619 

15.38 

1.18 

1.979 

14 

0.22 

0.083 

87.9 

9.43 

0.638 

15.84 

1.15 

1.938 

15 

0.26 

0.097 

102.8 

9.64 

0.652 

16.19 

1.13 

1.903 

16 

0.33 

0.120 

127.6 

10.10 

0.683 

16.95 

1.11 

1.857 

17 

0.42 

0.152 

162.4 

10.49 

0.709 

17.61 

1.08 

1.819 

18 

0.52 

0.190 

202.1 

10.98 

0.742 

18.43 

1.06 

1.774 

19 

0.74 

0.271 

288.5 

11.64 

0.787 

19.53 

1.01 

1.694 

20 

1.00 

0.367 

390.8 

12.23 

0.827 

20.54 

0.95 

1.589 

21 

1.30 

0.476 

506.5 

12.85 

0.869 

21.57 

0.86 

1.440 

22 

1.59 

0.583 

620.8 

13.35 

0.903 

22.41 

0.76 

1.281 

23 

1.92 

0.704 

749.4 

13.88 

0.939 

23.31 

0.64 

1.073 

24 

2.30 

0.843 

897.9 

14.33 

0.969 

24.06 

0.49 

0.823 

25 

2.84 

1.042 

1109.4 

14.72 

0.996 

24.72 

0.25 

0.412 

26 

3.82 

1.402 

1493.3 

14.79 

1.000 

24.83 

0.11 

0.183 
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Equilibrium  TBL  X-wire  survey 


Test  date: 

Reduction  date: 

Probe: 

Atmospheric  pressure: 

Static  temperature: 
Freestream  velocity,  Ue : 
Shear  velocity,  uT: 
Boundary-layer  thickness,  5: 


13-Aug-94,  run  1 
05-Sep-94 

TSI  1249A-10,  #51100 
101.6  KPa  (29.99  inHg) 
300.3  K  (540.5  R) 

16.13  m/s  (  52.93  ft/s) 
0.64  m/s  (  2.11  ft/s) 
2.95  cm  (  1.161  in) 


Survey 

station 

y 

(cm) 

v/s 

y+ 

U 

(m/s) 

^ rms 

(m/s) 

ms/ 

^rms 

(m/s) 

^ rms / 

01 

0.06 

0.022 

26.2 

8.31 

1.45 

2.258 

0.85 

1.318 

02 

0.07 

0.025 

30.9 

8.57 

1.44 

2.235 

0.83 

1.284 

03 

0.08 

0.028 

34.6 

8.76 

1.41 

2.197 

0.80 

1.237 

04 

0.10 

0.034 

41.9 

9.12 

1.37 

2.131 

0.74 

1.158 

05 

0.12 

0.041 

49.7 

9.36 

1.33 

2.074 

0.71 

1.104 

06 

0.16 

0.053 

64.9 

9.75 

1.28 

1.984 

0.67 

1.043 

07 

0.19 

0.065 

79.6 

10.05 

1.25 

1.935 

0.66 

1.023 

08 

0.23 

0.077 

94.2 

10.28 

1.22 

1.893 

0.65 

1.004 

09 

0.28 

0.093 

113.6 

10.54 

1.20 

1.862 

0.64 

1.000 

10 

0.35 

0.119 

145.0 

10.91 

1.17 

1.820 

0.64 

0.996 

11 

0.42 

0.142 

173.3 

11.34 

1.14 

1.778 

0.64 

0.990 

12 

0.52 

0.177 

214.7 

11.77 

1.13 

1.750 

0.64 

0.996 

13 

0.74 

0.250 

303.7 

12.57 

1.06 

1.643 

0.64 

0.989 

14 

1.00 

0.340 

413.6 

13.22 

0.99 

1.532 

0.61 

0.949 

15 

1.31 

0.444 

540.3 

13.93 

0.89 

1.382 

0.56 

0.877 

16 

1.61 

0.545 

663.4 

14.52 

0.79 

1.224 

0.51 

0.798 

17 

1.92 

0.650 

790.6 

15.02 

0.67 

1.037 

0.45 

0.692 

18 

2.30 

0.779 

947.6 

15.55 

0.50 

0.774 

0.35 

0.543 

19 

2.73 

0.926 

1125.7 

15.91 

0.30 

0.460 

0.24 

0.380 

20 

3.42 

1.158 

1408.4 

16.09 

0.13 

0.198 

0.12 

0.186 

21 

4.51 

1.530 

1860.2 

16.12 

0.10 

0.157 

0.07 

0.115 

22 

6.42 

2.175 

2645.0 

16.14 

0.09 

0.146 

0.07 

0.101 
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Disturbed  TBL  single-wire  survey 

Test  date:  24-Oct-94 

Reduction  date:  27-Oct-94 

Probe:  TSI  1261A-T1.5,  #914074 

Atmospheric  pressure:  101.1  KPa  (29.86  inHg) 

Static  temperature:  299.0  K  (538.2  R) 

Freestream  velocity,  Ue:  13.74  m/s  (  45.08  ft/s) 

Shear  velocity,  uT:  0.32  m/s  (  1.05  ft/s) 

Boundary-layer  thickness,  5:  4.05  cm  (  1.593  in) 


Survey 

station 

y 

(cm) 

y/s 

y+ 

U 

(m/s) 

u/ue 

U+ 

'U’rms 

(m/s) 

01 

0.02 

0.004 

3.1 

1.97 

0.143 

6.16 

1.02 

3.174 

02 

0.02 

0.005 

4.4 

2.52 

0.183 

7.87 

1.19 

3.724 

03 

0.03 

0.007 

5.5 

2.75 

0.200 

8.58 

1.25 

3.919 

04 

0.04 

0.009 

7.5 

3.00 

0.218 

9.38 

1.31 

4.103 

05 

0.05 

0.011 

9.3 

3.21 

0.234 

10.04 

1.36 

4.240 

06 

0.05 

0.013 

11.2 

3.40 

0.248 

10.63 

1.39 

4.337 

07 

0.07 

0.016 

13.5 

3.58 

0.260 

11.18 

1.41 

4.418 

08 

0.09 

0.021 

17.4 

3.80 

0.276 

11.87 

1.43 

4.484 

09 

0.11 

0.028 

22.8 

4.02 

0.292 

12.55 

1.46 

4.548 

10 

0.13 

0.032 

26.7 

4.15 

0.302 

12.96 

1.46 

4.578 

11 

0.20 

0.049 

40.2 

4.44 

0.323 

13.89 

1.50 

4.688 

12 

0.23 

0.057 

47.3 

4.57 

0.332 

14.27 

1.53 

4.768 

13 

0.27 

0.067 

55.6 

4.69 

0.341 

14.66 

1.55 

4.839 

14 

0.33 

0.082 

68.0 

4.87 

0.355 

15.24 

1.58 

4.954 

15 

0.47 

0.116 

96.3 

5.50 

0.400 

17.18 

1.76 

5.492 

16 

0.71 

0.176 

145.9 

6.63 

0.482 

20.71 

1.95 

6.097 

17 

0.90 

0.223 

184.9 

7.58 

0.552 

23.69 

2.00 

6.263 

18 

1.30 

0.322 

266.1 

9.55 

0.695 

29.84 

1.73 

5.398 

19 

1.59 

0.394 

325.9 

10.74 

0.781 

33.56 

1.32 

4.141 

20 

1.94 

0.480 

397.5 

11.69 

0.851 

36.53 

0.92 

2.865 

21 

2.30 

0.569 

471.0 

12.35 

0.899 

38.59 

0.72 

2.240 

22 

2.88 

0.711 

587.9 

13.04 

0.949 

40.76 

0.47 

1.482 

23 

3.83 

0.945 

782.1 

13.58 

0.988 

42.45 

0.14 

0.424 

24 

5.10 

1.259 

1041.7 

13.71 

0.997 

42.84 

0.08 

0.237 

25 

6.17 

1.526 

1262.4 

13.74 

1.000 

42.95 

0.07 

0.214 
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Disturbed  TBL  X-wire  survey 


Test  date: 

Reduction  date: 

Probe: 

Atmospheric  pressure: 

Static  temperature: 
Freestream  velocity,  Ue: 
Shear  velocity,  uT: 
Boundary-layer  thickness,  6: 


08-Oct-94,  run  2 
27-Oct-94 

TSI  1249A-10,  #51100 
101.6  KPa  (30.00  inHg) 
296.1  K  (533.0  R) 

16.18  m/s  (  53.09  ft/s) 
0.45  m/s  (  1.47  ft/s) 
4.30  cm  (  1.691  in) 


Survey 

station 

y 

(cm) 

y/s 

y+ 

U 

(m/s) 

^ 77725 

(m/s) 

^rros/ 

Vrms 

(m/s) 

Vrms/ U- 

01 

0.07 

0.015 

19.3 

5.53 

1.61 

3.599 

0.78 

1.782 

02 

0.09 

0.020 

24.9 

5.71 

1.63 

3.654 

0.83 

1.854 

03 

0.11 

0.026 

32.6 

5.87 

1.67 

3.738 

0.88 

1.970 

04 

0.13 

0.030 

38.2 

5.94 

1.69 

3.773 

0.91 

2.043 

05 

0.20 

0.046 

57.1 

6.18 

1.72 

3.856 

1.04 

2.319 

06 

0.23 

0.054 

67.5 

6.30 

1.75 

3.924 

1.10 

2.453 

07 

0.27 

0.063 

79.4 

6.45 

1.79 

3.993 

1.16 

2.600 

08 

0.33 

0.077 

97.2 

6.66 

1.83 

4.103 

1.25 

2.806 

09 

0.43 

0.100 

125.0 

7.11 

1.96 

4.382 

1.38 

3.081 

10 

0.52 

0.122 

152.9 

7.55 

2.04 

4.567 

1.45 

3.245 

11 

0.60 

0.140 

175.1 

7.91 

2.09 

4.669 

1.47 

3.293 

12 

0.74 

0.173 

217.4 

8.67 

2.17 

4.845 

1.49 

3.322 

13 

0.88 

0.205 

257.5 

9.46 

2.17 

4.845 

1.46 

3.262 

14 

1.10 

0.255 

320.2 

10.71 

2.04 

4.552 

1.35 

3.007 

15 

1.30 

0.303 

380.3 

11.75 

1.80 

4.022 

1.19 

2.650 

16 

1.59 

0.371 

465.6 

13.02 

1.32 

2.960 

0.94 

2.104 

17 

1.92 

0.448 

561.7 

13.92 

0.94 

2.100 

0.70 

1.570 

18 

2.30 

0.536 

672.3 

14.64 

0.71 

1.588 

0.53 

1.178 

19 

2.85 

0.664 

833.6 

15.39 

0.49 

1.090 

0.36 

0.812 

20 

3.77 

0.877 

1100.0 

15.96 

0.17 

0.377 

0.17 

0.372 

21 

4.90 

1.141 

1432.1 

16.09 

0.11 

0.246 

0.09 

0.192 

22 

5.93 

1.380 

1731.5 

16.18 

0.10 

0.230 

0.07 

0.165 

23 

7.44 

1.732 

2174.1 

16.26 

0.10 

0.221 

0.07 

0.147 
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Appendix  B 

Measurement  uncertainty 


Measurement  uncertainties  were  estimated  using  the  principles  set  forth  in  [B.l,  B.2], 
In  general,  the  uncertainty  in  a  measurement  was  composed  of  a  combination  of  fixed 
error  or  bias,  B ,  and  random  error  or  precision,  P.  The  root-sum-square  model  was 
used  to  estimate  the  uncertainties  at  the  95  percent  confidence  level: 


Urss 


±[B2  +  P2j'/2 

±[b2  +  (<95s)2]1/2 


(B.l) 

(B.2) 


where  Urss  is  the  uncertainty,  S  is  the  sample  standard  deviation,  and  tg5  is  the 
95th  percentile  point  for  the  two-tailed  Students  “t”  distribution  (95  percent  confi¬ 
dence  interval).  For  sample  sizes  greater  than  30,  t95  is  considered  equal  to  2.  Bias 
and  precision  errors  were  propagated  through  to  calculated  results  individually,  then 
combined  into  overall  uncertainties.  Traceability  of  working  standards  to  the  Na¬ 
tional  Institute  of  Standards  and  Technology  (NIST)  was  established  through  the 
Navy  Metrology  and  Calibration  (METCAL)  Program  and  through  manufacturer- 
supplied  calibrations. 

The  estimated  uncertainties  are  presented  in  Table  B.l.  In  general,  two  contrib¬ 
utors  to  the  bias  were  considered:  the  bias  of  the  working  standard  against  which 
the  instrument  was  calibrated,  and  the  bias  of  the  calibration  as  determined  through 
the  actual  measurement  system.  The  total  bias  was  then  estimated  as  the  root-sum- 
square  of  the  two  contributors.  Correlated  biases  were  handled  as  outlined  in  [B.2]. 
This  procedure  proved  important  for  the  velocity  measurements,  where  the  single 
thermocouple  measuring  the  tunnel  static  temperature  was  also  used  to  determine 
the  calibration  reference  temperature.  The  correlated  biases  cancel,  allowing  an  in¬ 
strument  with  relatively  large  bias  to  be  used.  The  precision  was  estimated  from 
in-place  readings  using  the  entire  measurement  system  as  P  =  tg5S.  The  combi¬ 
nation  of  uncertainty  components  by  root-sum-square  has  the  effect  that  any  single 
component  less  than  about  10  percent  of  the  maximum  contributor  can  be  neglected. 
Negligible  contributors  are  indicated  by  a  dash  in  Table  B.l. 
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Table  B.l:  Instrumentation  specifications  and  estimated  uncertainties. 


Quantity 

(units) 

Type 

B 

P 

Urss 

Atmospheric 
pressure  (mmHg) 

Welch  Scientific 
barometer 

2.5e-3 

2.5e-3 

Pitot/static 
difference  (mmHg) 

Edwards  High  Vacuum 
570D-10T-2C2-V1X 

5.5e-3 

4.0e-3 

6.8e-3 

Atmospheric 
temperature  (deg  C) 

mercury-in-glass 

thermometer 

1.0 

— 

1.0 

Tunnel  temperature,  T 
(deg  C) 

Omega  DP  11  type  K 
thermocouple 

2.2 

0.2 

2.2 

Fluctuating  wall 
pressure,  p  (Pa) 

Briiel  &  Kjaer  4138 

— 

5.0e-3 

5.0e-3 

u  velocity 
(m/s) 

DISA/TSI  single  wire 
DISA/TSI  X-wire 

0.025 

0.030 

0.010 

0.004 

0.027 

0.030 

v  velocity 
(m/s) 

DISA/TSI  X-wire 

0.030 

0.004 

0.030 

Probe  location,  y 
(m) 

telescope  with 

Ono  Sokki  DG-357 

2.8e-5 

2.8e-5 

4.0e-5 

Power  spectra 
(dB) 

HP  3562A  dynamic 
signal  analyzer 

0.2 

0.03 

0.2 
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Appendix  C 

Localized  windowed  peak  detection  algorithm 


This  appendix  contains  the  source  code  for  the  Localized  Windowed  Peak  Detection 
(LWPD)  algorithm.  The  code  is  written  in  MATLAB.  A  few  notes  are  provided  for 
those  not  familiar  with  the  notation.  Comments  are  preceeded  by  the  °/.  character. 
The  first  contiguous  block  of  comment  lines  are  displayed  by  the  help  facility.  The 
notation  a  (m ,  n)  denotes  a  matrix  a  with  m  rows  and  n  columns.  The  :  character  is 
used  to  designate  a  range,  as  in  ni=l:m. 


function  table=lwpd(data,tw) 

'/.LWPD  Localized  Windowed  Peak  Detection  (LWPD)  Algorithm. 
‘/,  table=lwpd[data,tw) 

*/, Input  Parameters : 

‘/.  data  — time  record 

‘/.  tw  — window  size  in  samples 

‘/.Output  Parameters: 

*/.  table  — Table  for  LWPD  events.  It  contains: 

*/.  column  1  — trigger  index 
’/,  column  2  — start  index 
’/.  column  3  — end  index 
’/,  column  4  — peak  index 
% 

*/,  Note:  This  function  calls  PEVENT  at  each  recursion  level. 

‘/,  By  P.  Penafiel 
*/.  April,  1995 

°/,  Last  revision:  April,  1995 

[zO , flag] =pevent (data) ; 
while  flag 

[zO,flag]=pevent(zO,tw) ; 
end 

table=[zO( : ,2) ,z0( : ,3:4) ,z0( :  ,1)]  ; 
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function  [z,flag]=pevent(y,xO) 

7PEVENT  One  pass  in  the  recursion  of  the  LWPD 

7. 


algorithm. 


’/.Input  Parameters: 

%  y 

’/.  xO 

^Output  Parameters : 
'/.  z 
7. 


7. 

7. 


7. 


7. 


7.  flag 
7. 


time  record 

window  size  in  samples 

Table  for  LWPD  events 

level.  It  contains: 

column  1 

column  2 

column  3 

column  4 

0  if  end  of  recursion, 
1  otherwise, 


;  the  present  recursion 

trigger  index 
start  index 
end  index 
peak  index 


7.  By  P.  Penafiel 
7.  Mar/ 1995 

7.7.  Comments  beginning  with  7.7.  added  by  M.  Kammeyer,  18  Aug  95 


[m,n]=size(y) ; 

z=zeros  (m,4) ;  7.  z(:,3)  start;  z(:,4)  end 

f lag=0 ; 

nj=0; 

7.7.  If  first  pass  through  LWPD,  find  the  local  peaks  and  write  them 
/./  to  the  keeper  matrix,  z.  Start  and  end  of  each  peak  event  is  one 

7.7.  point  shy  of  the  local  minimum  on  either  side  of  the  peaks, 
if  nargin==l 
if  n  ~=  1 

error(’input  argument  should  be  a  column  vector’); 
end 

for  ni=2:m-l  7.7.  for  time  indices  2  thru  m-1 

if  (y(ni)>y (ni+1) )  &  (y(ni)>y(ni-l))  7.7.  if  y(ni)  is  a  local  peak 
nj=nj+l; 

z(nj,l)=ni;  %7.  write  the  time  and  value 

z(nj ,2)=y(ni) ;  7.7.  to  the  keeper  matrix,  z 

nk=ni ; 

while (y(nk)>y(nk+l)) 
nk=nk+l ; 

if  (nk==m)  nk=nk+l;  break;  end 
end 

z(nj ,4)=nk-l; 
nk=ni ; 

while (y(nk)>y (nk-1)) 
nk=nk-l ; 

if  (nk==l)  nk=nk-l;  break;  end 
end 

z(nj ,3)=nk+l ; 
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end 
end 
if  nj 
flag=l; 
end 


7.7. 


If  in  the  recursion  loop  of  LWPD 
else if  nargin==2 

for  ni=l:m  7.7.  Loop  over  time  indices 

if  ni==l  7.7.  If  first  peak 


if  (y(ni+l,l)-y(ni,l))<xO 
if  (y(ni,2)>y(ni+l,2)) 
nj=nj+l; 

z(nj , l:3)=y(ni,l:3) ; 
z(nj ,4)=y(ni+l ,3) ; 
flag=l; 
end 
else 

nj=nj+l ; 

z(nj , 1 :4)=y (ni ,1:4); 
end 


7.7. 

7.7. 

7.7. 


if  first  peak  is  not  isolated 
but  it  is  a  local  peak 
then  keep  it 


7.7. 


if  first  peak  is  isolated 
7.7.  keep  it 


elseif  ni==m  7.7.  If  last  peak 
if  (y(ni,l)-y(ni-l,l))<xO  7.7. 

if  (y(ni ,2)>y(ni-l,2) )  7.7. 

nj=nj  +  l;  7.7. 

z(nj , 1 :2)=y(ni ,1:2) ; 
z(nj ,3)=y(ni-l ,4) ; 
z(nj ,4)=y(ni,4) ; 
flag=l; 
end 

else  7.7. 

nj=nj+l;  77. 

z(nj ,l:4)=y(ni,l:4) ; 
end 


if  last  peak  is  not  isolated 
but  it  is  a  local  peak 
then  keep  it 


if  last  peak  is  isolated 
keep  it 


else  7.7.  all  interior  peaks 

if  (y(ni,l)-y(ni-l,l))<xO 
if  (y(ni+l,l)-y(ni,l))<xO 


7.7. 

7.7. 


if  this  peak  is  not  isolated 
on  either  side 


if  (y(ni,2)>y(ni-l,2))  &  (y(ni,2)>y(ni+l,2))  7.7.  and  it  is  a 


nj=nj+l; 
z  (nj , 1 : 2) =y (ni ,1:2) ; 
z(nj ,3)=y(ni-l,4) ; 
z(nj ,4)=y(ni+l ,3) ; 
flag=l; 
end 
else 

if  (y(ni,2)>y(ni-l,2)) 
nj=nj+l; 

z (nj ,1:2) =y (ni ,1:2) ; 


7.7.  local  peak 

7.7.  then  keep  it 


7.7.  if  this  peak  is  isolated 

7.7.  from  the  NEXT  ONLY, 

7.7.  AND  is  greater 

7.7.  than  the  previous  peak,  then 
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z(nj ,3)=y(ni-l,4) ; 
z(nj ,4)=y(ni,4) ; 
flag=l; 
end 
end 
else 

if  (y(ni+l,l)-y(ni,l))<xO 
if  (y (ni ,2)>y (ni+1 ,2) ) 
nj=nj+l; 

z(nj , 1 :3)=y (ni ,1:3) ; 
z(nj ,4)=y(ni+l ,3) ; 
f lag=l ; 
end 
else 

nj=nj+l ; 

z(nj , l:4)=y(ni,l:4) ; 
end 
end 

end 

end 

else 

error (’Maximum  2  input  arguments 
end 

z=z(l :nj , :) ; 


VL  keep  it 


if  this  peak  is  isolated 
7X  from  the  PREVIOUS  ONLY, 
7.7.  AND  is  greater 
7.7.  than  the  next  peak,  then 
7.7.  keep  it 


yX  if  this  peak  is  isolated 
7.'/.  on  both  sides,  keep  it 


are  accepted’); 
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